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FOREWORD 


The  Army  Research  Office  is  pleased  to  sponsor  the  workshop,  "Dynamic 
Response  of  Composite  Structures"  to  be  held  in  New  Orleans,  August  30-  September  1, 
1993.  The  purpose  of  the  workshop  is  to: 

(i)  assess  the  state  of  the  art  in  the  modeling  and  analysis  of  modem  composite  structures 

that  are  relevant  to  the  Army's  need  to  improve  the  design-and  operation  of  land 
vehicles,  helicopters,  weapon  systems,  etc.  _ 

(ii)  identify  new  opportunities  for  relevant  research. 


The  special  areas  of  interest  are  crash  worthiness  of  composite  structures,  composite 
helicopter  rotor  blade  dynamics,  stability  and  control,  smart  composite  structures,  etc. 
Other  area  of  composite  structural  research  are  welcomed  and  the  authors  are  encouraged  to 
show  the  relevance  of  their  work  to  the  U.S.  Army. 

The  workshop  contains  six  keynote  and  distinguished  lectures  and  34  lectures  for  a 
total  of  40  paper  presentations.  The  regular  sessions  include  (i)  Vibration  of  Composite 
Structures  (ii)  Smart  Structures  (iii)  Impact  of  Composite  Structures  (iv)  Damage  and 
Failure  (v)  Design  of  Composite  Structures  (vi)  Structural  Mechanics.  We  are  pleased  with 
the  many  outstanding  presentations  at  this  ARO  workshop  and  would  like  to  acknowledge 
the  financial  and  moral  support  of  the  Army  Research  Office.  Thanks  are  due  to  the 
contributions  from  the  authors,  the  reviewers  and  the  support  and  encouragement  of  Dr. 
John  Crisp,  Dean  of  Engineering,  University  of  New  Orleans. 


David  Hui 

ARO  workshop  chairman 
University  of  New  Orleans 
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One-dimensional  phase  transition  problems  in  shape- 
memoiy  alloys  and  certain  polymers  modeled  as  nonlinear 
elastic  solids  have  been  extensively  studied,  e.g.  see  Ericksen 
(1975),  and  the  solution  corresponding  to  a  phase  transition 
between  two  stable  states  has  been  referred  to  as  necking, 
Hutchinson  and  Neale  (1983).  Physicists  (e.g.  Jacobs  (1985)) 
working  in  the  area  of  ferroelasticity  have  examined  phase 
transformations  in  two-dimensional  problems.  In  two- 
dimensional  solidification  problems,  it  has  been  found  (Davis, 
1990)  that  the  straight  interface  between  the  solid  and  liquid 
phases  will  become  unstable  and  deform  into  a  cellular  one 
under  certain  conditions,  and  may  even  develop  into  a  more 
complicated  dendritic  pattern. 

Here  we  consider  a  fiber-reinforced  flat  sheet  of  unit 
thickness  with  fibers  aligned  along  the  x-direction,  and  the  end 
surfaces  of  the  sheet  subjected  to  uniform  surface  tractions  in 
the  x-direction.  We  assume  that  the  sheet  undergoes 
infinitesimal  deformations  so  that  linear  kinematic  relations 
apply,  and  the  fibers  are  densely  packed  so  that  lateral 
deformations  of  the  sheet  are  negligible.  Thus  every  point  of 
the  sheet  undergoes  a  time-dependent  displacement  u  in  the 
x-direction,  and  u  is  in  general  also  a  function  of  x  and  y. 
Whereas  both  the  fibers  and  the  matrix  material  are  taken  to 
be  elastic,  the  stress-strain  curve  for  the  former  is  taken  to  be 
nonlinear  with  two  stable  branches,  cf.  Fig.  1,  and  that  for  the 
latter  linear.  The  material  of  the  fibers  is  presumed  to  exhibit 
phase  transformations.  Also,  fibers  are  assumed  to  cany  all  of 
the  tensile  load  and  the  matrix  material  supports  the  shear 
stress.  Thus 


°xy  ~  2G«xy  '  G  ^  “  G%  ' 

Here  <r„  and  are  the  nonzero  components  of  the  Cauchy 
stress  tensor,  eu  and  are  nonvanishing  components  of  the 
infinitesimal  strain  tensor,  W  is  the  strain-energy  density  for 
the  fiber  material  normalized  to  vanish  for  null  value  of  «D, 
and  G  is  the  shear  modulus  of  the  matrix  material.  The  axial 
displacement  u  is  governed  by 

Punt  '  f/(*«)u-«  +  Guw  <2> 

We  assume  that  p  >  0  and  P(«n)  >  0-  Here  P  denotes  the 
derivative  of  f  with  respect  to  its  argument.  We  note  that  eqn. 
(2)  differs  from  that  in  a  solidification  problem  wherein  the 
governing  equation  is  parabolic.  The  second  term  on  the 
right-hand  side  of  eqn.  (2)  describes  the  coupling  effect 
between  adjacent  fibers. 


o=f(e) 


Fig.  1.  The  stress-strain  curve  for  the  shape-memory  fibers. 


Since  fiber  material  is  presumed  to  exhibit  phase 
transformations,  it  is  conceivable  that  there  is  a  surface  of 
discontinuity,  or  a  singular  surface,  in  the  sheet  in  the  sense 
that  fibers  are  in  different  phases  on  the  two  sides  of  this 
surface.  Across  this  surface,  the  displacements  must  be 
continuous,  i.e., 

I«1  -  0  (3) 

where  (  u  ]  »  u  +  -  u"  equals  the  difference  in  the  values  of  u 
on  the  positive  and  negative  sides  of  the  surface;  the  positive 
side  of  the  surface  being  the  one  on  which  the  outward  normal 
points  in  the  direction  of  propagation  of  the  surface.  If  this 
surface  propagates  with  speed  V  in  the  x-direction  then  the 
Rankin  e-Hugo  niot  condition  (see  Hutchinson  and  Neale,  1983) 
and  the  balance  of  total  energy  require  that 


-pV  K1  •  if(u,x)l,  (4) 

-V  IW  ♦  1  fi( u„)2l  -  lf(u,x)u„l.  (S) 

We  consider  a  necking  solution  that  is  independent  of  y, 
and  with 

z  »  x  -  Vt  (6) 

can  be  expressed  as 

u(x,t)  -  U(z).  (7) 


According  to  Mullins  and  Sekerka’s  method  (e.g.  see 
Davis,  1990),  we  consider  a  perturbation  of  the  interface 
geometry  in  the  reference  frame  {z  «  x- Vt,y,  t}  moving  with 
the  speed  V  (cf.  Fig.  2).  Thus,  let  the  interface  geometry  be 
perturbed  to  that  given  by 

z-  z  *  (y,t)  «  6 coskye"1,  (8) 

where  S  is  an  infinitesimal  number. 

We  consider  a  few  simple  cases  for  which  V  *  0  and 
examine  the  interface  stability  of  a  propagating  neck.  For  the 
following  three  cases,  the  interface  z  «  0  is  found  to  be  stable 
in  the  sense  of  Mullins  and  Sekerka. 

(i)  k  »  0.  This  corresponds  to  a  one-dimensional  problem. 

(ii)  There  is  no  surface  tension,  and  the  two  modulii 
corresponding  to  the  stretches  on  either  side  of  the  neck 
are  equal  to  each  other. 
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unr«eck«ri 
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Fig.  2.  A  schematic  sketch  of  a  propagating  neck  in  a 
uni  axially  reinforced  sheet. 

(iii)  The  matrix  material  is  very  weak  so  that  the  shear 
modulus  G  can  be  approximated  as  zero. 

(iv)  The  fibers  are  extremely  light  but  are  quite  strong. 
However,  the  stationary  interface  is  always  unstable. 
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Delamination  failure  is  a  major  considera¬ 
tion  in  the  design  of  laminated  composite  struc¬ 
tures.  The  impact  of  foreign  objects  is  one  of  the 
common  causes  of  delamination  in  composite 
structures.  However,  because  impact-induced 
delamination  appears  inside  the  laminates,  it 
can  hardly  be  detected  by  visual  inspection.  There¬ 
fore,  locating  the  embedded  damage  is  critically 
important  for  maintaining  and  repairing  the  struc¬ 
tures  to  prevent  a  catastrophic  failure  in  service. 

In  this  presentation,  a  damage  diagnostic 
method  is  proposed  to  detect  the  location  and 
size  of  the  embedded  damage  in  composite  struc¬ 
tures,  using  measured  dynamic  responses  from 
distributed  piezoelectrics  built  into  the  struc¬ 
tures.  The  method  consists  of  three  major  parts: 
a  dynamic  structural  analyzer,  a  response  com¬ 
parator,  and  a  signal  generator  as  shown  in  Fig¬ 
ure  1.  The  system  would  function  as  follows: 


Damage  Detection  Method 


Figure  1  A  schematic  of  the  proposed  dam¬ 
age  detection  method. 


first,  upon  impact,  the  comparator  compares  the 
measured  responses  from  the  piezoelectrics  with 
the  calculated  responses  from  the  analyzer  to 
estimate  the  impact  energy  and  the  impact  lo¬ 
cation.  Second,  appropriate  diagnostic  signals 
will  be  selected  by  the  generator  and  then  dis¬ 
patched  from  designated  piezoelectric  actuators 
to  selected  piezoelectric  sensors.  The  measured 
diagnostic  signals  from  the  sensors  will  then  be 
compared  in  the  comparator  with  the  calcula¬ 
tions  from  the  analyzer  to  determine  the  damage 
location  and  size. 

Composite  plates  with  and  without  an  im¬ 
planted  delamination  were  fabricated  with  four 
pairs  of  piezoceramics  patched  on  the  surfaces. 
One  pair  of  piezoceramics  was  used  as  actua¬ 
tors,  and  the  other  pairs  were  treated  as  sensors. 
The  measured  dynamic  responses  between  the 
damaged  and  undamaged  plates  were  compared. 
Figure  2  presents  the  comparison  of  the  dynamic 
responses  measured  from  the  designated  sensors. 
Considerable  differences  in  dynamic  responses 
between  the  damaged  and  undamaged  plates  were 
detected,  suggesting  the  presence  of  damage  in 
the  plates. 

At  present,  a  method  has  been  developed 
for  detecting  the  location  and  size  of  a  delamina¬ 
tion  in  composite  beams.  The  predictions  com¬ 
pared  with  experimental  results  are  presented 
in  Figure  3.  The  technique  successfully  detected 
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OCLAMNATtO* 


SCNSOK 


the  presence  of  dei&minations  and  predicted  with 
reasonable  accuracy  the  locations  and  sizes  of 
the  damage  in  the  beams. 

A  method  has  also  been  proposed  to  esti¬ 
mate  impact  location  and  the  impact  force  his¬ 
tory  for  a  beam  subjected  to  a  transient  dy¬ 
namic  load.  Figure  4  shows  that  the  predicted 
impact  locations  based  on  two  sensor  measure¬ 
ments  agreed  with  the  data  very  well. 

It  is  anticipated  that  a  health  monitoring 
system  such  as  the  one  proposed  could  be  devel¬ 
oped.  However,  considerable  work  both  on  the 
development  of  the  methods  and  on  the  imple¬ 
mentation  of  the  technique  is  still  required. 


Figure  2  Comparison  of  sensor  measure¬ 
ments  of  composite  plates  with  and  with¬ 
out  a  delamination.  (Upper)  Piezoceramic 
measurements  from  the  sensor  at  upper  right 
corner.  (Below)  Piezoceramic  measurements 
from  the  sensor  at  lower  left  corner. 
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Figure  3  Comparison  of  predicted  delami¬ 
nation  locations  and  sizes  in  composite  beams 
with  actual  measurements. 


Figure  4  Comparison  of  impact  location  I  •• 
tween  the  predictions  and  the  actual  m<  > 


surements. 


Dynamic  Stresses  in  Composites  at  Low  Temperatures 
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EXTENDED  ABSTRACT 


Behavior  of  polymer  composites  at  low  temperature  is  a 
neglected  area  of  research,  and  the  dynamic  behavior  at  low 
temperatures  is  more  so.  ft  is  nc*  because  that  the  thermo-physical 
problems  of  research  at  low  temperatures  are  less  interesting  than 
the  high  temperature  problems  but  because  that  the  phenomenal 
development  in  technology  of  composites  in  the  past  had  been 
driven  by  the  needs  of  the  aerospace  industries  which  historically 
have  seen  more  problems  with  composites  at  high  temperatures  than 
at  low  temperatures.  Of  comae,  it  is  tine  that  the  cryogenic  range  of 
low  temperature  effects  on  composites  have  been  studied  in  the  past 
primarily  far  specific  applications,  but  modi  wide  spread 
applications  of  composites  through  the  newly  developed  land 
vehicles  and  water  crafts  and  new  applications  of  composites  for 
civil  engineering  structures  demand  a  reassessment  of  the  focus  of 
composites  research  on  the  dynamic  effects  at  low  temperatures. 

It  is  weD  known  that  the  lowering  of  temperatures  modifies 
the  behavior  of  composites.  Hie  two  basic  changes  that  happen  in 
flic  Mm™  ^ fwmpfuhiff  mtm  fh*  Aibjm  ill 

dflftic  modulus  of  the  individual  oompoocots  md  tocoodly  in  die 
changes  m  stress  hocwsc  of  outmatched  tfaennal  expansion 
coefficients  of  those  components*  Soperoupoeed  on  these  effects 
wiE  be  the  effects  of  higb-atrain-rate.  Figure  1  extracted  from 
Hartwig  (1979)  shows  the  dramatic  inflnmnr  of  both  the  low 
temperature  and  strain  rate  on  the  behavior  of  the  epoxy  resin 
matrix  of  composites,  which  would  effect  the  behavior  of  the 
composites.  The  figure  shows  that  low  temperature  and  high  strain 
rate  increases  the  Young’s  modulus  value,  E  of  the  polymer  matrix. 
Extending  the  Rosen's  (1965)  model  of  composites  failure  to  tensile 
loading  condition,  Dutta  (1992)  has  shown  that  the  maximum 
critical  fiber  stress  gfa(D,  that  would  be  indneed  in  the  composites 
is  given  fay 


indicative  of  the  new  increased  modulus  of  the  matrix  under  the 

combined  effects  of  low  temperature  and  high  strain  rate.  It  is 
obvious  from  equation  (1)  that  cra(T)  will  rapidly  increase  to  die 
value  of  the  fiber  ultimate  strength  (Figure  2)  if  the  values  of 
(E.(T0)  +  K(T0  -1))  in  the  numerator  increases  as  a  result  of 
temperature  decrease,  or  attain  rate  increase  or  both.  However,  no 
direct  experimental  results  are  currently  available  to  identify  the 
individual  effects  of  these  two  factors.  Currently  investigations  are 
underway  at  GRREL  to  determine  the  effects  of  tbeae  two  factors 
separately. 

One  of  the  exanqdes  of  these  two  combined  effects  on  the 
composite  laminate  behavior  is  seen  in  the  results  of  the  work 
performed  by  Dutta.  Hni  and  Altamiiano  (1991).  Figure  3  shows 
the  results  of  the  energy  absorption  tests  by  a  30-ply  AS4/3502 
quasi-iaotropic  laminated  composite  when  indented  by  a 
hemispherical  bit  under  the  influence  ofan  incident  stress  pulse  of 
only  180  microsecond  duration,  bis  seen  that  in  general,  at  higher 
velocities  of  impact,  energy  absorption  increases  both  at  room 
temperature  (2I°Q  and  at  low  temperature.  But  at  lower  velocities 
where  the  primary  mode  of  energy  absorption,  Ec,  is  by  hertzian 
contact  deformation,  as  outlined  by  Creszcznk  (1982) 


where 


Ee  “  (PVJK4/5K3*/16)m[(K1  +Kj)/(R/2)] 
K1«(l-ui2)/<JrEi) 

Kj-a-O/OfE.) 


(2) 


(3) 

(4) 


owCO-2 


VfEf(E„(Tn)  +  K(TC-T)) 
3(1-Vf) 


I03 


(1) 


where  Ef  -  fiber  stress,  Vf  -  fiber  volume  fraction,  E.  -  matrix 
■nnHniiif  To  -  low  temperature,  and  K  is  an  empirical  constant 


E,  and  u,  are  the  Young's  modulus  and  Poisson's  ratio  of  the  steel 
indent  or,  respectively,  Ez  is  the  laminate  modulus  in 
direction,  and  ur  is  the  laminate's  Poisson's  ratio.  The  value  Ee 

dearty  decreases  with  the  increased  values  of  Ez  as  a  result  of 
temperature  and  strain  rate. 


ssni* 


Similarly,  the  strain  rate  dependent  behavior  of  composites 
in  low  temperature  environment  is  also  important  for  predicting  the 
crashworthiness  of  composites  used  in  vehicles  or  guard  rails.  It  is 
necessary  that  serious  research  effect  is  undertaken  to  clearly 
understand  the  combined  effort  of  low  temperature  and  high  strain  <_ 

rate  on  composite  behavior.  w 

D 

to 
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fiber  critical  stress. 
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Abstract 

Thin-walled  structures  of  interest  to  the 
U.  S.  Army,  such  as  rifle  barrels,  rocket  casings, 
helicopter  blades  and  containment  vessels  are 
often  constructed  of  layers  of  anisotropic, 
filament  or  fiber-reinforced  materials  which  must 
be  designed  to  remain  elastic.  A  proper 
assessment  of  end  or  edge  effects  in  such 
structures  is  of  fundamental  importance  to  U.  S. 
Army  technology.  The  extent  to  which  local 
stresses,  such  as  those  produced  by  fasteners  and 
at  joints,  can  penetrate  girders,  beams,  plates  and 
shells  must  be  understood  by  the  designer.  Thus 
a  distinction  must  be  made  between  global 
structural  elements  (where  Strength  of  Materials 
or  other  approximate  theories  may  be  used)  and 
local  elements  which  require  more  detailed  (and 
more  costly)  analyses  based  on  exact  elasticity. 
The  neglect  of  end  effects  is  usually  justified  by 
appeals  to  some  form  of  Saint-Venant’s  principle 
and  years  of  experience  with  homogeneous 
isotropic  elastic  structures  has  served  to  establish 
this  standard  procedure.  Saint-Venant’s  principle 
also  is  the  fundamental  basis  for  static  mechanical 
tests  of  material  properties.  Thus  property 
measurements  are  made  in  a  suitable  gage  section 
where  uniform  stress  and  strain  states  are  induced 
and  local  effects  due  to  clamping  of  the  specimen 
are  neglected  on  invoking  Saint-Venant’s 
principle.  Such  traditional  applications  of  Saint- 
Venant’s  principle  require  major  modifications 
when  strongly  anisotropic  and  composite 
materials  are  of  concern.  See  e.g  list  of 
references  and  the  references  cited  therein.  For 
such  materials,  local  stress  effects  persist  over 
distances  far  greater  than  is  typical  for  isotropic 
materials.  In  this  lecture,  we  provide  a  survey  of 
situations  in  linear  elasticity  where  anisotropy 
and  material  inhomogeneity  induce  such  extended 
Saint-Venant  end  zones.  Both  static  and  dynamic 
problems  are  described,  although  in  the  latter  case 
relatively  few  results  are  known.  The  implication 


for  the  analysis  and  design  of  structures  using 
advanced  composite  materials  are  addressed. 
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ABSTRACT 

The  replacement  of  steel  with  advanced 
composite  materials  represents  one  of  the  most 
significant  advances  under  consideration  in  the 
field  of  Civil  Engineering.  Since  composite 
materials  are  unproven  in  their  substitution  for 
steel  in  concrete  structures  there  is  considerable 
motivation  to  instrument  test  structures 
incorporating  these  materials  in  order  to  gain  a 
through  understanding  of  how  well  composite 
materials  serve  in  this  new  challenging  role. 
We  have  developed  a  fiber  optic  Bragg  grating 
strain  gauge  that  is  well  suited  for  embedding 
within  composite  structures  and  has  many  other 
attractive  attributes.  These  include  the  ability  to 
follow  the  local  strain  with  a  very  high 
frequency  response  while  at  the  same  time 
tracking  creep.  Lastly,  these  sensors  are 
immune  to  electrical  interference. 

Many  civil  engineering  structures  like:  bridges, 
dams,  tunnels....  could  benefit  from  a 
structurally  integrated  fiber  optic  sensing 
system  that  could  monitor  the  state  of  the 
structure  throughout  its  working  life.  Such  an 
integrated  sensing  system  could  also  monitor 
the  various  structural  components  during 
construction,  possibly  leading  to  improved 
quality  control.  In  certain  instances  this 
resident  sensing  system  could  also  provide 
valuable  information  on  the  usage  of  the 
structure.  For  example,  in  the  case  of  a  bridge 
the  sensing  system  could  be  used  to  generate 
information  on  the  traffic  itself. 

Two  types  of  carbon  prestressing  tendons  are 
being  tested  in  a  new  two-span  concrete  road 
bridge  that  is  under  construction  in  Calgary.  It 
was  originally  planned  that  UTIAS  would 


insirumcni  six  of  these  girders  3ragg 

grating  laser  sensors,  two  with  \  rope, 
two  with  Leadline  and  two  with  steel. 

In  the  intracore  Bragg  grating  fiber  optic  sensor 
the  centre  (Bragg)  wavelength  of  the  reflected 
signal  is  linearly  dependent  upon  the  product  of 
the  scale  length  of  its  periodic  index  variation 
and  the  mean  core  index  of  refraction.  Changes 
in  strain  or  temperature,  to  which  the  optical 
fiber  is  subjected,  will  shift  this  Bragg 
wavelength  leading  to  a  spectrally  encoded 
optical  measurement  that  can  he  determined  by 
the  simple  wavelength  dependent  ratiometric 
technique  that  we  haye  recently  developed.  We 
have  recently  demonstrated  that  the  sensing 
Bragg  grating  can  be  used  to  control  the 
wavelength  of  a  laser,  and  for  the  Calgary 
bridge  project  we  have  developed  a  4-channel 
Bragg  grating  fiber  laser  sensor  demodulation 
system.  This  system  involves  4  -  independent 
Bragg  grating  tuned  fiber  lasers  that  are 
multiplexed  in  order  to  be  pumped  by  one 
semiconductor  laser.  Each  fiber  laser  is  tuned 
by  a  connectorized,  remote  and  structurally 
embedded  Bragg  grating  that  serves  as  a 
sensor.  The  measurand  controlled  wavelength 
of  each  fiber  laser  is  determined  by  a  passive 
ratiometric  wavelength  demodulation  system, 
pioneered  in  the  UTIAS  Fiber  Optic  Smart 
Structures  Laboratory.  The  resulting  4-channel 
Bragg  grating  fiber  laser  sensor  demodulation 
system  with  a  laptop  computer  is  rugged, 
compact  and  transportable  to  construction  sites 
where  it  has  been  used  to  undertake 
measurements  on  large  concrete  girders  as  used 
for  civil  engineering  applications.  A  schematic 
of  the  system  is  presented  as  figure  L 


To  test  Che  performance  of  Ihc  Bragg  grating  the  surface  of  the  prestressing  Tokyo  ropes,  ai 

laser  strain  sensor,  Bragg  gratings  have  been  roughly  the  same  location  along  the  tendon, 

adhered  to  the  surface  of  cantilever  beams  and  Figure  2  shows  the  the  strain,  as  measured  by 

the  wavelength  of  the  fiber  laser  measured  as  a  the  Bragg  grating  laser  sensor,  plotted  against 

function  of  the  applied  strain.  Good  linearity  is  the  applied  load  to  the  girder, 

observed  provided  the  strain  is  limited  to  about 

3000  pstrain.  Once  this  is  exceeded  stress  In  the  case  of  the  two-span  concrete  road 

induced  birefringence  in  the  optical  fiber  leads  bridge  in  Calgary  20  Bragg  gratings  were 

to  the  appearance  of  a  second  peak  in  the  laser  embedded  into  two  concrete  girders  prestressed 

spectrum.  This  undesirable  affect  can  be  with  steel,  two  prestressed  with  Tokyo  rope 

avoided  if  polarization  maintaining  optical  fiber  and  one  girder  prestressed  with  Leadline.  In 

is  used.  ~  each  girder  there  is  at  least  one  Bragg  grating 

temperature  sensor  and  two  Bragg  grating 
strain  sensors.  At  the  time  of  writing  the  bridge 
At  the  University  of  Manitoba  a  series  of  static  is  under  construction  so  we  hope  to  provide 

and  fatigue  loading  tests  were  performed  on  more  details  on  the  use  of  the  Bragg  grating 

concrete  girders  that  were  prestressed  with  the  fiber  laser  sensing  system  in  monitoring  the 

graphite  Tokyo  Rope  and  Leadline.  Resistive  strains  experienced  by  the  steel  and  graphite 

foil  strain  gauges  was  adhered  to  the  surface  of  prestressing  tendons  at  the  Workshop  itself, 

each  of  the  prestressing  cables  at  roughly  their 
midpoint.  Bragg  gratings  were  adhered  only  to 


Figure  2. 
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Design  of  advanced  aircraft,  ground-vehicle,  and 
naval  structures  with  multi-layered  organic-matrix  com¬ 
posite  materials  requires  easy-to-use  computational 
tools  which  provide  accurate  predictions  of  structural 
response  and  integrity  under  mechanical,  thermal,  and 
dynamic  loads.  Unlike  metals,  composites  commonly 
possess  significantly  lower  strength  and  stiffness  in  the 
direction  transverse  to  the  reinforcing  fiber  as  compared 
to  the  fiber  direction;  hence,  transverse  shear  and  nor¬ 
mal  strain  and  stress  effects  can  significantly  contribute 
to  the  response  and  failure  of  such  structures.  Conse¬ 
quently,  computational  tools  that  are  employed  in  the 
industrial  environment  must  accurately  predict  these  so- 
called  "secondary”  stains  and  stresses.  Conventional 
shear-defonnable  theories  that  are  prevalent  in  com¬ 
mercial  codes  generally  do  not  provide  accurate 
through  the  thickness  deformations  and  stresses. 
Further,  nearly  all  higher-order  theories  in  the  literature 
are  impractical  for  general  computational  purposes 
because  they  do  not  result  in  finite  elements  that  are 
compatible  with  standard  analysis  software  nor  with 
standard  engineering  boundary  conditions. 

The  paper  first  examines  recent  advances  in  the  area 
of  clastostatic  and  dynamic  analysis  of  thick  laminated 
composite  plates  and  shells  with  the  specific  focus  on 
identifying  suitable  theories  for  general-purpose  finite 
element  software  and  industrial  applications.  Advan¬ 
tages  and  shortcomings  of  first-  and  higher-order  theo¬ 
ries  are  discussed  both  from  analytical  and  computa¬ 
tional  viewpoints. 

A  special  type  of  Higher-Order  Theory  (HOT)  for 
plate  and  shell  analysis  which  is  compatible  with  exist¬ 
ing  structural  analysis  software  products  and  standard 
engineering  boundary  conditions  is  closely  examined.  l* 
3  The  theory  is  based  on  independent  tfarough-the- 
thickness  expansions  of  displacements,  transverse  shear 
strains,  and  transverse  normal  stress.  Unlike  other 
higher-order  theories,  it  uses  a  least-squares  constraint 
on  the  compatibility  of  deformations  and  transverse 


strains  which  also  permits  exact  satisfaction  of  upper 
and  lower  surface  stress  equilibrium  conditions.  The 
virtual  work  principle  (and,  in  dynamics,  Hamilton's 
variational  principle)  is  applied  to  derive  the 
appropriate  equations  of  equilibrium/motion  and  engi¬ 
neering  boundary  conditions.  As  analytical  results  indi¬ 
cate,  the  theory  is  capable  of  accurately  predicting  dis¬ 
placements,  natural  frequencies,  mode  shapes,  and 
three-dimensional  strains  and  stresses.  It  yields  im¬ 
proved  predictions  over  the  first-order  theory,  particu¬ 
larly  for  thick  laminates  and  higher  vibration  frequen¬ 
cies.  Moreover,  thickness-stretch  modes  and  natural 
frequencies,  which  are  critical  for  delamination  initia¬ 
tion  and  which  are  suppressed  in  the  first-order  and 
many  higher-order  theories,  are  accurately  predicted. 

Figure  1  shows  a  typical  frequency  spectra  plot  for  a 
simply-supported  orthotropic  square  plate  where  the 
HOT  frequencies  for  the  lowest  six  modes  correspond¬ 
ing  to  the  wave  numbers  (umiMM)  am  compared  with 
those  of  three-dimensional  elasticity  theory  (designated 
as  EXACT). 


Fig.  1.  Frequency  spectra  for  simply-supported, 
orthotropic  square  plate  for  six  lowest-frequency  thick¬ 
ness  modes  for  wave  numbers  (m,nWl,l). 


As  can  be  seat,  results  agree  very  favorably  with  the 
available  exact  solutions  even  in  the  very  thick  regime. 
The  corresponding  mode  shapes  for  a  thick  plate  are 
shown  in  Fig.  2.  Though  not  shown,  these  modes  are 
also  closely  correlated  with  the  exact  ones. 


Antisymmetric  Modes  Symmetric  Modes 
Ha) 


Fig.  2.  Six  lowest-frequency  mode  shapes  correspond¬ 
ing  to  wave  numbers  (nuiMl.l)  for  simply-supported, 
orthotropic  thick  square  {date  (2h/L*0.3). 

From  the  perspective  of  finite  element  approxima¬ 
tions  suitable  for  general  purpose  finite  element  soft¬ 
ware,  this  higher-order  theory  is  shown  to  rely  on 
simple  C°-cootinuous  displacement  assumptions  that 
permit  efficient  and  robust  finite  element  formulations. 
The  appropriate  approximations  for  elements  of  this 
type  are  the  same  as  those  for  conventional  first-order 
elements  that  have  been  explored  extensively  by  die 
author4  and  others.  This  computational  capability 
enables  large-scale  analyses  of  complex  multi-layered 
composite  structures,  providing  improved  response  and 
failure  predictions  using  standard  engineering  software. 
This  is  achieved  without  the  additional  computational 
cost  that  is  usually  associated  with  the  use  of  higher- 
order  theories. 

A  conforming  three-node  facet  shell  element  and  a 
companion  two-node  beam  element  have  been  formu¬ 
lated  and  implemented  in  the  general-purpose  code 
COMET  (Computational  MEchanics  Testbed)  so  that 
any  structural  configuration  can  be  considered.  The 
shell  element  has  also  been  used  in  the  form  of  a  user- 
supplied  capability  in  the  commercial  code  ABAQUS. 

Additional,  ‘auxiliary’  computational  techniques 
which  serve  to  improve  recovery  of  three-dimensional 
displacements  and  mode  shapes,  strains,  and  stresses 
are  also  discussed,  and  some  of  their  key  features  are 


emphasized.  These  techniques  are  generally  applied  at 
the  post-processing  stage  of  the  finite  element  analysts 
and  address  the  following  areas: 

(1)  Improved  ply-by-ply  recovery  of  displacements 
and  inplane  strains  and  stresses.3 

(2)  Recovery  of  transverse  stresses/strains  via  inte¬ 
gration  of  three-dimensional  equations  of  equilibrium. 

(3)  Recovery  of  physically  smooth  strain  and  stress 
quantities  and  their  partial  derivatives  by  way  of  the 
Smoothing  Element  Analysis  (SEA).6  The  latter  quanti¬ 
ties  are  required  in  order  to  perform  integration  of 
three-dimensional  equations  of  equilibrium  in  (2). 

(4)  Estimation  of  residual-equilibrium  discretization 
errors  in  the  finite  element  analysis  in  order  to  guide  in 
adaptive  mesh  refinement  procedures.  Such  physically 
meaningful  error  estimators — which  require  partial 
derivatives  of  stress  resultants  in  plate  and  shell  analy¬ 
ses — are  readily  computed  by  way  of  SEA. 

Finally,  remarks  are  offered  on  the  potential  use  of 
the  aforementioned  computational  capabilities  in  the  ar¬ 
eas  of  global-local  finite  element  modeling,  nonlinear 
dynamic  response,  impact  damage,  and  progressive 
failure  in  thick  composite  plate  and  shell  stnictrures. 
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sions  in  the  original  theory  that  give  rise  to  the  problemat¬ 
ical  C1  continuity  requirement  In  order  to  ensure  the 
continuity  of  transverse  shearing  stresses,  a  constraint  con¬ 
dition  is  then  added  to  the  statement  of  the  displacement 
field  and  its  associated  functional,  so  the  new  form  of  die 
theory  can  be  shown  to  have  the  same  properties  as  the 
original  theory  of  DiSciuva.  By  imposing  such  a  constraint 
via  the  penalty  function  concept,  the  transverse  shear 
stresses  are  made  continuous  without  introducing  addi¬ 
tional  degrees  of  freedom.  Furthermore.  C°  continuity  of 
all  degrees  of  freedom  is  maintained,  allowing  for  a  sim¬ 
ple,  efficient,  and  accurate  finite  element  model  of  the  gov¬ 
erning  equations. 

The  finite  element  model  based  on  the  above  laminate 
theory  takes  advantage  of  the  interdependent  interpolation 
concept  introduced  by  Tessler  and  Hughes  (Comp.  Meth. 
AppL  Mech.  Eng.,  Vol.  39,  pp.  311-335,  1983).  As  a  result, 
the  linear  part  of  the  constraint  is  satisfied  identically,  the 
approximation  order  of  the  element  is  effectively 
increased  without  adding  more  degrees  of  freedom,  and 
the  possibility  of  shear  locking  is  eliminated.  The  pro¬ 
posed  discrete-layer  finite  element  model  is  thus  much 
more  efficient,  convenient  and  widely  applicable  than  any 
one  developed  to  date. 


Introduction 

The  currently  available  layerwise  (or  discrete  layer) 
theories  provide  significant  improvements  over  the  equiv¬ 
alent  single-layer  (or  smeared)  theories  for  analysis  of  thin 
and  moderately  thick  laminated  structures.  For  practical 
applications,  however,  all  existing  discrete-layer  theories 
give  rise  to  finite  element  models  that  are  either  too  com¬ 
putationally  expensive  or  poorly  suited  for  general  pur¬ 
pose  analysis.  In  order  to  take  advantage  of  this  powerful 
class  of  laminate  theories,  a  new  generalized  discrete-layer 
theory  and  an  associated  computational  model  for  the 
analysis  of  moderately  thick  laminated  beams  has  been 
developed1.  The  theory  assumes  a  piecewise  linear  distri¬ 
bution  of  the  inplane  displacements  through  the  thickness 
of  the  laminate,  and  satisfies  the  stress  continuity  condi¬ 
tions  at  the  interface  between  each  layer.  In  this  way,  the 
number  of  degrees  of  freedom  in  the  theory  is  made  inde¬ 
pendent  of  the  number  of  layers  in  the  laminate. 

Unlike  available  theories,  the  current  theory  is  spe¬ 
cially  formulated  to  be  ideally  suited  for  solution  by  the 
finite  element  method.  A  f'O-noded  beam  element  has 
been  developed  that  combines  the  penalty  method  with  an 
interdependent  interpolation  scheme,  resulting  in  a  simple 
and  efficient  C°  element  that  is  rank  sufficient,  accurate 
for  both  thick  and  thin  beams,  and  has  no  apparent  defect 
The  accuracy  of  the  computational  model  has  been  dem¬ 
onstrated  for  static  and  vibration  analysis  of  thin  and  thick 
laminated  beams,  with  and  without  delaminations  and/or 
ply  damage.  It  has  been  shown  that  layerwise  variations  of 
the  inplane  displacements  must  be  taken  into  account 
when  modelling  either  the  global  or  local  response  of  dam¬ 
aged  laminates. 

Approach 

In  the  current  theory,  the  continuity  requirements  on 
the  displacement  variables  of  the  discrete-layer  laminate 
theory  of  DiSciuva  (J.  Sound  and  Vib.,  Vol.  105,  pp.  425- 
442,  1986)  are  reduced  by  momentarily  relaxing  the  con¬ 
dition  of  continuous  transverse  shear  stresses.  Two  new 
degrees  of  freedom  are  introduced  to  replace  the  expres- 


1  Extension  of  the  theory  and  computational 
model  to  plates  and  shells  is  currently  underway. 


Numerical  Results 

A  simply-supported  (0/d/90/0)  beam  has  been  ana¬ 
lyzed,  where  d  represents  a  full-length  delamination,  mod¬ 
elled  as  a  very  thin  compliant  layer.  The  normalized  center 
deflection  due  to  a  sinusoidal  load  is  plotted  versus  the 
span-to- thickness  ratio  in  Figure  1.  The  center  deflections 
are  normalized  by  those  of  the  exact  solution  obtained  by 
Pagano  (J.  Comp.  Mat.,  Vol.  3,  pp.  398-411.  1969).  For 
beams  with  an  aspect  ratio  less  than  500,  the  first-order 
shear  deformation  theory  (FST)  and  the  third-order  theory 
of  Lo,  Christenson,  and  Wu  (LCW)  do  not  accurately  pre¬ 
dict  the  global  deflection  behavior,  while  the  cunent  dis¬ 
crete-layer  model  (ZZ)  is  nearly  exact  for  all  thickness 
ratios  greater  than  8. 

The  reason  for  the  poor  predictions  of  the  equivalent 
single  layer  theories  (FST  and  LCW)  is  illustrated  in  Fig¬ 
ure  2,  where  the  inplane  displacement  at  the  end  of  the 
beam  is  plotted  versus  the  normalized  thickness  coordinate 


W/W  (exact) 


for  the  case  L/ h  =  10.  It  can  be  seen  in  Figure  2  that  a 
piecewise  linear  variation  of  displacement  through  the 
thickness  of  a  laminate  is  necessary  to  adequately  model 
the  kinematics  of  a  laminate  with  a  delamination  or  any 
type  of  damage  resulting  in  adjacent  layers  with  greatly 
different  shear  stiffnesses.  The  predictions  of  the  current 
model  are  indistinguishable  from  those  of  the  exact  solu¬ 
tion  of  Pagano.  Predictions  of  the  through-thickness  varia¬ 
tion  of  inplane  stress  components  (not  shown)  are  equally 
as  accurate  using  the  current  model. 

The  layerwise  kinematics  also  affect  the  dynamic 
response  predictions,  as  seen  in  Figure  3.  Here,  the  ratio  of 
the  first  four  natural  frequencies  of  the  delaminated  beam 
predicted  by  FST  and  ZZ  are  plotted  versus  the  span-to- 
thickness  ratio.  For  moderately  thick  beams,  the  FST  pre¬ 
dicts  the  first  four  natural  frequencies  of  the  beam  to  be 
approximately  two  times  greater  than  those  predicted  by 
ZZ.  The  difference  in  the  predictions  increases  as  the 
mode  number  increases,  and  decreases  as  the  beam 
becomes  thinner. 

These  results  have  implications  for  all  post-damage 
analysis  as  well  as  analysis  of  laminated  smart  structures. 
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Figure  2. 

Axial  displacement  versus  normalized  thickness 
coordinate  at  the  end  of  a  simply  supported 
(0/d/90/0)  beam  subjected  to  a  sinusoidal  load. 


Figure  1. 

Normalized  center  deflection  versus 
span-to-thickness  ratio  of  a  simply  supported 
(0/d/90/0)  beam  subjected  to  a  sinusoidal  load. 


Figure  3. 

Ratio  of  flexural  frequencies  predicted  by  FST  and 
ZZ  versus  span-to-thickness  ratio  for  a  simply 
supported  (0/d/90/0)  beam. 
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A  study  od  the  working  principles  and  structural  he  torque-plate  fin  design  of  Barrett  (19921—  )  used 

arrangements  of  active  composite  torque-plate  missile  fins  is  directionally  attached  piezoelectric  (DAP)  elements  bonded  to 

presented.  The  design  evolution  from  active  fin  twist  to  torque-  “  isotropic  substrate  antisymmetrically.  The  plate  was 

plate  design  shows  that  early  designs  proposed  by  optimized  to  provide  the  maximum  twist  from  commercially 

experimenters  like  Crawley.  Lazarus  and  Warkemin  (1989)  available  stock.  However,  another  design  which  uses  symmetric 

provided  a  valuable  starting  point  for  a  host  of  new  designs.  elements  can  provide  significantly  higher  twist  deflections.  An 

This  early  design  of  adaptive  lifting  surface  was  composed  of  a  examination  of  the  plate  performance  using  laminated  plate 

bending-twist  coupled  graphite-epoxy  plate  which  had  lead  ‘ta«y  clearly  illustrates  the  potential  of  this  new  arrangement 

zhvonate  titanate  (PZT)  piezoelectric  sheets  laminated  on  either  Using  the  analysis  techniques  outlined  by  Barrett 

side  of  the  plate.  As  the  PZT  sheets  were  alternatively  (19921-23).  the  properties  of  DAP  lamina  can  be  estimated, 

contracted  and  extended  to  induce  a  bending  deformation,  the  Fundamentally,  they  are  modeled  as  having  different  stiffnesses 

plate  also  twisted.  The  plate  was  exposed  to  air  loads.  in  ***  longitudinal  and  lateral  directions,  but  with  equal 

Accordingly,  small  active  twist  rie format^  resulted  in  larger  actuation  strains.  For  a  symmetric  torque-plate  with  an  isotropic 

deflections  through  active  aeroservoelasticity.  substrate  and  a  main  spar  that  stiffens  the  plate  in  bending,  the 

This  first  successful  design  lead  to  the  development  of  the  laminate  twist  can  be  solved  for  as  shown  in  equation  1 
torque-plate  fin  configuration.  One  of  die  design  compromises 
that  was  necessary  for  the  bending-twist  coupled  plate  design  Kn 
was  to  make  it  thick  enough  so  that  bending  loads  could  be 


_pni(Bi»  +  B»).-p»i(Bii  +  l>ii). 
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carried.  Unfortunately,  as  the  thickness  was  increased,  the 

deflections  decreased.  So  a  new  design  was  conceived.  This  If  ^  P,ale  is  further  constrained  so  that  neither 

torque-plate  fin  design  alleviates  this  thickness-deflection  longitudinal  nor  transverse  betiding  is  allowed,  then  equation  1 

design  compromise.  Fig.  1  shows  the  evolution  to  the  torque  simplifies  and  can  be  broken  down  into  its  components  as 

plate  design  as  initially  proposed  by  Barren  (1992u-3).  This  shown  in  equation  2. 

K  =  "141.= 

design  is  strengthened  in  bending  by  a  spar  at  the  quarter-chord  12  0 


and  provides  much  more  control  than  the  twisting  {date  design 
as  the  entire  aerodynamic  shell  undergoes  a  pitch  change.  If  the 
spar  is  moved  aft,  then  aenxcrvoelastic  properties  may  also  be 


Active  Beading-Twist  Coupled  Plate  Stiffened  Torque  Plate  Torque  Plate  Fin 

•low  drag  •  high  twist  rate  •  high  pitch  deflections 

•  weak  in  beadiag  *  strong  in  bending  •  low  drag 

•  flatter  and  divergence  proee  •higbdrag  •  strong  in  beading 


Fig.  1  Design  evolution  from  active  plate  to  torque-plate  fm  configuration 


Similar  equations  can  be  derived  for  other  forms  of 
ortbotropic  actuator  elements  like  piezoelectric  fiber  composites 
(PFC).  PFCs  are  similar  to  DAP  elements  analytically  in  that 
they  have  differing  stiffnesses  longitudinally  and  laterally. 
However,  their  active  strains  are  also  unequal.  For  analysis,  a 
nondimensional  twist  parameter  is  defined  in  equation  3  so  that 
the  different  types  of  materials  may  be  compared  on  equal 
thickness  elements,  where  OR  is  the  ortfaotropy  ratio  of  the 
actuator  element,  AR  is  the  activation  ratio.  TR  is  the  thickness 
ratio  of  the  actuator  to  the  substrate,  and  SR  is  the  stiffness  ratio 
of  the  actuator  to  die  substrate. 
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Figure  2  shows  that  TP  values  approaching  1.5  are 
regularly  attainable  using  well  proven  assembly  techniques  and 
two-axis  bending  constraint  with  a  symmetric  layup. 
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Fig.  2  Twist  Parameter  for  a  symmetric  DAP  laminate 
with  an  uncoupled  substrate 

If  PFCs  are  used  to  induce  twist  in  a  similar  manner,  then 
intetdigitated  electrodes  must  be  used  to  achieve  TP  levels 
approaching  that  of  DAP  elements.  The  interdigitation  of  the 
electrodes  will  induce  actuation  strains  of  opposite  signs  in  the 
actuator  plies.  As  can  be  seen  in  Fig.  3,  only  at  very  high  fiber 
volume  fractions  and  very  large  negative  strain  ortho  tropics  can 
the  PFC  laminate  surpass  the  performance  of  the  DAP  laminate. 
However,  it  should  also  be  noted  that  if  the  modulus  stiffness 
ratio,  MSR,  is  increased  to  1,  then  the  !«»»««♦»  will  resemble  a 


sheet  of  PZT  in  stiffness  and  accordingly,  the  twist  induced  by 
the  PFC  composite  will  far  surpass  that  of  the  DAP  laminate 
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Fig.  3  Twist  parameter  for  a  symmetric  PFC  laminate 
with  an  uncoupled  substrate 

A  comparison  of  the  stiffnesses  of  the  two  materials  shows 
similar  trends.  (The  stiffness  is  important  in  maintaining  a  high 
flutter  and  divergence  speed. ) 

A  final  examination  of  the  performance  gains  that  are 
possible  through  the  new  structural  arrangement  is  also 
conducted.  This  study  demonstrates  that  the  amounts  of  twist 
deflection  obtained  by  Barren  (19921-2-3)  can  be  nearly  doubled 
and  the  torsional  stiffness  may  also  be  increased.  This  new 
discovery  indicates  that  aeroservoelastic  coupling  may  no 
longer  be  needed  for  induction  of  high  pitch  deflections  in 
missile  fins. 
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ABSTRACT 

Shap*  aaaory  fibers,  made  of  a  Nlckel-Tltaniua  alloy 
(NITINOL),  aro  embedded  Ins  Id*  a  now  class  of  SMART 
coaposlto  boaas  In  order  to  control  the  dynaale 
characteristics  and  the  daapod  response  of  those  beams 
«d»en  subjected  to  external  excitations.  The  NITINOL 
fibers  are  tuned  by  adjusting  their  initial  tension  and 
operating  teaperature  to  achieve  an  optiaal  balance 
between  the  the  real  aof  toning  of  the  coaposlto  aatrlx, 
the  stiffening  effect  iaparted  by  the  activated  fibers 
•nd  the  enhanced  daaping  of  the  aatrlx  as  It  Is  heated 
towards  Its  glass  transition  region. 

A  finite  eleaent  aodel  Is  developed  to  aodsl  the 
dynaalcs  of  daapod  NITINOL-relnforeed  coaposlto  booae. 
The  aodel  is  utilized  to  eeapute  the  natural 
frequencies,  the  nodal  loss  factors  and  the  frequency 
response  functions  of  this  class  of  SMART  boaas.  The 
frequency  responses  of  the  NITINOL-relnforeed  boaas  are 
conpared  with  those  of  the  un reinforced  boaas  In  ardor 
to  eaphaslze  the  laportance  of  the  NITINOL  reinforceaent 
and  Its  optiaal  tuning  in  significantly  attenuating  the 
vibration  ef  these  boaas. 

1.  INTRODUCTION 

Considerable  attention  has  boon  dovotad  recently  to 
th*  utl fixation  of  the.  Shape  Meaory  Nlckal-TXtaaiiai 
alloy  (NZTZMOL)  In  developing  SMART  eoaposltes  that  are 
capable  of  adapting  Intelligently  to  external 
disturbances  (  Baz  et  al.  1990  and  1991).  The  eaphaals 
has.  however,  been  placed  on  stiffening  the  NITINOL 
eoaposltes  through  proper  activation  of  the  shape  aaaory 
fibers.  Such  an  eaphasls  Is  particularly  laportant  whan 
the  static  characteristics,  as  the  critical  >»— »<  l-g 
loads,  aro  to  bo  controlled.  But.  when  the  dynaale 
response  Is  of  concern,  the  stiffening  aedtanlsa  alone 
be  cones  Ineffective  la  attenuating  vibrations  result!^ 
froa  broad  hand  excitations  unless  It  is  augmented  with 
a  control  lad  energy  dissipation  mechanism,  Fortunately. 
NITINOL  eoaposltes  have  a  unique  source  of  ■*— c'rg  which 
Is  built-in  the  coaposlto  aatrlx  Itself.  No  atteapt 
has.  however,  been  node  to  exploit  the  excellent 
features  of  this  built-in  daaping  source. 

It  is  therefore  the  purpose  of  this  paper  to 
simultaneously  utilize  both  the  stiffening  and  the 
energy  dissipation  aochanisas  to  achieve  optiaal 
vibration  attenuation  ever  brood  frequency  spectrun. 
The  nature  of  Interaction  between  the  two  aechanlaaa 
will  be  Investigated  and  the  conditions  that  ensure 
optiaal  balance  between  then  will  be  datarained. 

2.  DYNAMICS  OF  DAMPED  NtTZMOL-BRNnHCZD  COMPOSITE  BEAKS 
The  dynaale  characteristics  of  a  NITINOL-relnforeed 
heaa  eleaent.  shown  la  Figure  Cl).  Is  governed  by.  the 
following  equation  of  notion: 

CM.)  til  ♦  (K.J  t«J  -  IF).  (1)  ' 

where  (g)  and  (1)  denote  the  nodal  daflactlon  and 
acceleration  vectors.  The  stiffness  aatrlx  IK.]  of  the 
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Fig.  (1)  -  NITINOL-relnforeed  been  eleaent. 


eleaent  are  given  by: 

t  t 

(K.)  -  I.  J  IDITC0I  dx  -  P,  J  CCITIC)  dx.  (2) 
o  o 

where  E^I.  la  the  flexural  rigidity  of  the  heaa  eleaent. 
<  la  Its  length,  end  P.  la  the  net  axial  fores  acting  on 
the  eleaent.  Aleo.  the  eatrlcee  [Cl  and  (0)  are 
determined  froa  the  transverse  displaces ent  function  w 
aa  fellows: 

w  •  (A)tdl.  dw/dx  e  [C)(«]  and  d*u/dx*»  (01(d).  (3) 

The  elaaanta  M.(l.  J)  of  the  boss  aatrlx  (H.J  of  the 
booa  eleaent  are  obtained  using  tbs  consistent  aass 
formulation  aa  follows: 

t 

H.C1.J)  -  p.  A,  J  [A,]  (Ajl  dx  (4) 


where  (A,  1  and  (A.)  are  tbs  1th  and  Jth  elaaanta  of  the 
vector  (A]  given  by  equation  (3).  Aloe  In  aquation  (4). 
p.  and  A.  denote  the  density  and  cross  sectional  area  of 
the  coaposlto  booa.  In  aquation  (2).  the  force  Pa  •  \ 
♦  P*  -  Tt  where  r.  Pt  and  Tt  denote  the  axial 
ae  chan  ire  l  loads,  Jthe  sxlal  thermal  lead  and  the  total 
tension  developed  by  the  NITINOL  fibers  respectively. 
The  thermal  lead  Pt  •  «  H  ^  A,,  la  generated  by  the 
teaperature  difference  AB  caused  by  the  activation  aid 
de-activation  of  the  NITINOL  fibers.  Also,  a  and  C. 
donate  the  thermal  expansion  coefficient  of  the 
composite  and  Its  modulus  of  elasticity. 

It  Is  laportant  hero  to  note  that  the  activation 
and  deactivation  of  the  NITINOL  fibers  vary,  on  one 
hand,  the  tension  Tt  and  control  the  stiffness  ef  the 
eoapmlte  and  coopensste  for  the  softening  of  the  matrix 
due  to  the  associated  boating.  On  the  other  hand,  it 
controls  the  doaplag  characteristics  of  tbs  eooposlte. 
This  Is  attributed  to  the  fact  that  the  modulus  of 
elasticity  ^  of  the  composite  bean  is  s  complex  modulus 
(  t*  •  V  (1  *1  x)J  whose  storage  modulus  V  and  leas 
factor  e  aro  controlled  by  tbo  teaperature  and  frequency 
ef  tbo  coaposlto. 


3.  PERFORMANCE  OF  DAMPED  NIT  I  MOL -REINFORCED  BEAMS 

Material*  The  characteristics  of  NITINOL-rainforced 
beaas  are  detsrainad  for  a  coaposite  beaa  ude  of 
randoaly  oriented  glass  fibers  eabedded  In  a  low  cure 
temperature  polyester  resin.  The  beaa  Is  30  ca  long. 
2.  S  ca  wide  and  0.  1S6  ca  thick  mounted  In 
claaped-claaped  arrangement.  The  temperature  and 
frequency  dependence  of  the  complex  modulus  of  the  beam, 
is  shown  In  Figure  (2).  Such  reduced  frequency  nomograph 
Is  obtained  experimentally.  The  displayed 
characteristics  show  typical  behavior  of  structural 
composites  where  the  storage  aodulus  E'  decreases  as  the 
operating  temperature  Is  Increased.  It  shows  also  that 
the  damping  characteristics  of  composites  Increase 
significantly  as  their  operating  temperature  is 
Increased  towards  the  glass  transition  region. 

In  NITINOL  composites.  It  Is  possible  to  compensate 
for  the  degradation  of  the  storage  aodulus  of  the 
composite  with  temperature  by  tuning  the  NITINOL  fibers 
properly.  The  fibers  can  be  tailored  to  produce  their 
phase  recovery  forces  to  counterbalance  the  softening  of 
the  coaposite  matrix  with  Increased  temperature.  Such  an 
Important  feature  when  combined  with  the  enhanced 
damping  at  high  temperature  can  be  extremely  effective 
In  attenuating  structural  vibrations. 

In  the  present  study,  four  NITINOL  55  fibers  that 
are  0.S5  mm  In  diameter,  are  eabedded  inside  the  beaa 
through  vulcanized  rubber  sleeves  that  have  outer 
diameter  of  0. 95  mm. 

The  characteristics  shown  In  Figure  (2)  is  utilized 
along  with  the  finite  element  model  to  compute  the 
frequency  response  of  NITINOL  composites  at  different 
temperatures  and  Initial  tensions. 
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Fig.  (2)  -  Reduced  frequency  nomograph  for  the 
fiberglass/resln  coaposite  beaa 


Vibration  control  The  optimal  tuning  of  the  Initial 
tension  and  operating  temperature  of  the  NITINOL  fibers 
is  determined  In  order  to  minimize  the  maximum  amplitude 
of  vibration  of  the  coaposite  bean  when  subjected  to 
external  excitations.  Such  amplitude  of  vibration  is 
measured  by  the  maximum  frequency  response  function 
(FRF)  which  Is  determined  from: 

«•„  -  «a*  T  (d„(l)d11(J)l  /  -  »*  ♦  t  s.  wj  (5) 
!,),«  ■■ 


which  Is  given  by  equation  (l).  In  equation  (S),  u 
denote  the  excitation  frequency  of  a  force  located  at  J. 

Figure  (3)  shown  the  effect  of  the  temperature  and 
Initial  tension  on  the  aaxiaum  FRF  Indicating  that  there 
Is  an  optimal  operating  temperature  for  each  value  of 
the  Initial  tension.  At  that  temperature,  the  maximum 
FRF  attains  Its  minimum  value.  Figure  (4)  shows  the 
FRF,  at  beaa  aid-span  due  to  excitation  at  the  same 
location,  as  a  function  of  the  excitation  frequency  for 
optimally  tuned  and  activated  fibers  at  Initial  tensions 
of  13.34  and  31.14  N/flber  respectively. 
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Fig.  (3)  -  Contours  of  lso-FRF. 


Fig.  (4)  -  Effect  of  tension  *  frequency  on  FRF. 


4.  CONCLUSIONS 

The  dynamic  characteristics  of  damped 
NITIHOL-reinf oread  coaposite  beaas  have  been  presented. 
The  fundamental  Issues  governing  the  behavior  of  this 
new  class  of  SMART  composites  have  been  Introduced. 
Particular  emphasis  Is  placed  on  the  utilization  of  both 
the  stiffening  and  the  energy  dissipation  mechanisms, 
which  are  inherent  to  NITINOL  composite,  to  achieve 
optimal  vibration  attenuation  over  broad  frequency 
spectrum.  The  nature  of  interaction  between  the  two 
mechanisms  Is  Investigated  and  the  conditions  that 
ensure  optimal  balance  between  them  Is  determined. 
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where  ♦„(!),  w„  and  %  denote  the  nod#  shape  at  location 
1  for  the  nth  noda,  the  natural  frequency  of  the  nth 
mode  and  the  loas  factor  at  the  nth  noda  respectively.  •» 
The  three  parameters  dn(l),  «„  and  a„  era  obtained  from 
the  solution  of  the  eigenvalue  problem  of  the 
homogeneous  equation  of  notion  of  the  composite  beam 
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The  classical  theory  of  composite  beam 
vibrations  fails  to  predict  a  correct  dynamical 
response  if  high  frequencies  (oiWG/p  -  1,  co  - 
frequency,  h  -  beam  thickness,  G  -  shear 
modulus,  p  -  mass  density)  and  short  waves  ( h/1 
■  1,  1  -  wave  length)  are  involved.  We  are 
developing  an  improved  theory  of  composite 
beam  vibrations  which  captures  die  major 
phenomena  unaccessible  to  classical  theory. 

Probably  the  first  attempt  to  extend  the 
classical  one-dimensional  beam  theory  was 
made  by  Rayleigh,  taking  into  account  the 
inertia  of  cross  sectional  motion  to  describe 
the  beam  behavior  more  properly  under 
impact  load.  Then,  Timoshenko  generalized 
the  classical  one-dimensional  beam  theory 
by  introducing  shear  deformation  to  describe 
die  stress  state  of  short  waves.  Although 
corrections  by  Rayleigh  and  Timoshenko  are 
physically  meaningful,  quantitative  results 
were  usually  poor. 

A  more  fundamental  base  for 
construction  of  improved  one-dimensional 
theory  was  created  by  Mindlin.  He  found 
the  key  reasoning  to  improve  classical 
theory:  the  theory  should  predict  branches  of 
the  dispersion  curve.  The  more  branches  are 
predicted,  die  better  the  accuracy  of  the 
theory.  Classical  theory  describes  only  low- 
frequency  branches.  Mindlin  suggested  a 
method  to  incorporate  high-frequency 
branches,  but  his  approach  contained  some 
free  parameters  which  should  be  chosen 
from  fitting  the  dispersion  curves  of  ID 
theory  and  the  exact  dispersion  curves. 
Unfortunately,  the  exact  dispersion  curves 
are  known  only  for  simplest  cases,  such  as 
an  isotropic  homogeneous  plate  or  an 
isotropic  beam  with  circular  cross-section. 
A  method  of  derivation  of  ID  theory  which 
automatically  guarantees  coincidence  of 


dispersion  curves  of  ID  theory  and  exact 
dispersion  curves  has  been  proposed  by  V. 
Berdichevsky  about  15  years  ago.  Later,  it 
was  applied  to  developing  an  improved  ID 
theory  of  isotropic  homogeneous  (or  layered) 
plate,  shells  and  beams  by  V.  Berdichevsky, 
Le  Khank  Chao,  M.  Riazantzeva  and  S. 
Kvashnina.  The  method  has  proven  to  be 
very  efficient  In  the  study  in  progress  now, 
we  are  going  to  apply  this  method  to 
developing  a  ID  dynamical  theory  of 
anisotropic  inhomogeneous  beam.  In  this 
talk,  some  preliminary  results  are  presented. 

The  key  point  is  to  find  the  cross- 
sectional  modes  for  high-frequency  branches. 
A  code  was  developed  to  find  these  modes 
numerically  for  beams  with  arbitrary 
anisotropy  and  inhomogeneity.  The  big 
surprise  of  the  first  runs  is  that  for  some 
typical  types  of  anisotropy,  the  first  high 
frequency  mode  is  not  the  transverse  shear 
mode,  as  is  assumed  in  Timoshenko’s 
theory,  but  some  in-plane,  cross-sectional 
deformation.  That  explains  why 
Timoshenko’s  theory  does  not  work  well  for 
composite  beams.  A  parametric  study  for 
the  first  high-frequency  modes  of  a  thin 
strip,  anti-symmetric  angle  ply  laminate  has 
been  conducted.  Four  nonclassical  modes 
with  the  lowest  frequencies  are  shown  in 
Figure  1. 
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A  method  of  structural  analysis  based  on  the 
spline-sublayer  approach  is  presented  This  method 
was  developed  to  understand  the  mechanical  behavior 
of  inhomogeneous  structures,  particularly  composites, 
as  existing  commercial  packages  do  not  satisfy  the  fun¬ 
damental  requirements.  In  particular,  the  method  is  ca¬ 
pable  of: 

•  analysis  of  laminated  composite  structures  providing 
continuous  through  the  thickness  distribution  of  normal 
stresses 

•  analysis  of  textile  composites  with  elastic  properties 
changing  through  the  volume  of  the  composite 
•analysis  of  progressive  damage  of  composites 

For  the  correct  prediction  of  the  behavior  of 
structures  both  the  conditions  of  continuity  far  dis¬ 
placements  and  transverse  stresses  must  be  satisfied  at 
each  interface  of  the  structure.  The  continuity  of  trans¬ 
verse  stresses  can  be  satisfied  only  if  transverse  strains 
are  discontinuous  at  each  interface,  as  follows  from 
stress-strain  relationships.  This  leads  to  the  discontinu¬ 
ity  of  the  first  derivatives  of  the  displacements.  Hence, 
only  those  kinematic  models  which  incorporate  these 
discontinuities  are  correct  [  1].  Three  examples  are  con¬ 
sidered  in  this  paper  to  demonstrate  the  power  of  de¬ 
scribed  technique: 

A  three-layer  [0/90]  s  graphite/epoxy  plate 
simply  supported  along  its  four  side  edges  is  considered 
for  progressive  damage  analysis.  The  square  plate  has 
length  a  and  thickness  h.  All  three  layers  are  of  the  same 
thickness.  A  normal  load  is  applied  on  the  top  (z*h)  sur¬ 
face  as 

q*(x,y)  *  -qo  sin(nx/a)  sin(jty/b)  (1) 


where  qo  >  0.  The  first  step  of  the  initial  failure  analy¬ 
sis  is  to  define  the  load,  point  and  mode  of  initial  fail¬ 
ure.  After  a  computer  ’’scanning”  of  the  plate  volume 
is  applied  to  all  six  stress  components  normalized  by 
corresponding  ultimate  stress  values  (and  designated 
further  with  an  asterisk),  the  greatest  is  obtained.  In 
this  particular  example,  the  maximum  stress  criterion 
was  first  satisfied  at  qo=0.779  Msi  for  the  oy*  compo¬ 
nent  at  x  =  y  =  0.5a,  and  z=0.  Then,  in  the  parallele¬ 
piped  zone  around  this  point  in  which  the  inequality 

Oy*>0.98  (2) 

is  satisfied,  the  moduli  E 22,  G12  and  G32  have  been 
reduced  (10  times  for  this  particular  calculation). 
Equation  (2)  defines  the  damaged  zone  extending 
over  the  region  0.45  <  x/a  <  0.55  and  0.45  <  y/a  < 
0.55,  the  damaged  zone  being  taken  as  the  bottom 
layer.  The  next  step  in  the  gradual  failure  analysis  re¬ 
quires  one  to  define  the  new  distribution  of  the 
normalized  stresses  and  the  dimensions  of  tire  zone 
damaged  try  that  stresses.  Then  the  stiffnesses  in  this 
zone  have  to  be  reduced.  After  this  is  done  using  the 
same  approaches  as  in  the  first  step,  the  new  stress 
distributions  are  obtained  around  the  new  damaged 
zone.  The  variation  of  Oy*  on  z  and  x  coordinates  is 
shown  in  Figure  1.  The  stress  concentration  is  higher 
than  after  the  first  step,  but  is  more  localized.  This  in¬ 
dicates  that  the  gradual  failure  process  will  stop  ex¬ 
tending  after  several  steps  of  die  analysis  due  to  de¬ 
creasing  Oy*  in  both  x  and  y  directions.  In  order  to 
continue  the  analysis,  the  applied  load  has  to  be  in¬ 
creased.  Then  at  some  step  of  the  analysis  another 
failure  mods  will  be  revealed. 


Fig.  1.  Variation  af  die  normalized  stress  Oy*  along  z  coadmun  at 
zMISa.  yO.Sa  (a)  and  along  x  coordinate  at  yO 3a,  vO  (b)  for 
the  nnd«naged  plate  (curve  1)  and  the  damaged  plate  (23.4)  with 
3  different  ithnwiiiAw  mcthei 

A  double  box  beam  spar,  exhibiting  quasi- 
unbalanced  behavior  was  considered.  It  exhibits  twist- 
bend  coupling  without  extension-shear  coupling.  Both 
experimental  and  analytical  results  have  been  obtained 
far  twist/bend  ratio  of  die  spar  for  different  reinforce¬ 
ment  angles.  The  geometry  of  the  spar  was  as  follows:  a 
* 900 mm,  b*  46  mm,  h*48mm,d*  1  mm.Theelastic 
properties  of  die  Uniweave  fabric  composite  are:  E]  * 
12*  GF*.  Et  - 12.1  GPU,  vu  -  0.26,  va  -  0.2,  Gu  *  4.56 
GPa,  and  Gu  »  5.01  GPa.  The  analytical  and  experimen¬ 
tal  data  arc  compared  in  Figure  2.  The  analysis  was  per¬ 
fumed  for  die  length-to-wall  thickness  ratio  ISO 
(instead  of 900),  to  reduce  the  number  of  degrees  of  free¬ 
dom  employed,  while  operating  at  the  maximum  accept¬ 
able  element  aspect  ratio.  Despire  this  simplification  of 
the  model,  die  results  match  well  to  die  experimental 
data.  Not  only  are  the  general  trends  die  same,  but  also 


Fig.  2.  Analytical  and  experimental  data  for  the  twist/bend 
ratio  at  the  wiom  reinforcement  angles  of  the  box  beam  spar 


the  twist/bend  ratio  values  for  the  range  of  reinforce¬ 
ment  orientation  values  are  close.  It  should  be 
pointed  out  that  not  one  of  the  previous  theoretical 
models  was  able  to  account  for  a  continuous  varia¬ 
tion  of  elastic  properties  around  a  composite  box 
beam  cross  section  as  well  as  for  different  ’’non-clas¬ 
sical”  effects. 

As  a  final  example,  the  analysis  of  stress 
fields  in  a  cantilever  laminated  composite  plate  un¬ 
der  transverse  and  in-plane  loads  is  considered.  A 
square  plate  with  a/h*5  is  composed  of  three  uni- 
directionally  reinforced  layers  having  identical 
thicknesses  and  fiber  orientation  along  the  x-axis  in 
die  top  and  bottom  layers  and  along  the  y-axis  in  the 
middle  layer.  Some  results  for  the  case  of  transverse 
loading  (1)  are  presented  in  Figure  3.  It  is  seen  that 
longitudinal  normal  stress  ox  suffers  discontinuities 
at  the  interfaces  between  0  and  90  degree  layers, 
while  the  transverse  stresses  and  oz  don't  show 
any  visible  discontinuities. 


Fig.  3.  Variations  of  o*  (a),  ■*»  (b),  and  a*  (c)  along  z  for 
a  three-layer  cantilever  plate  under  transverse  load 
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Rate  Effects  in  Composite  Materials  * 
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Fiber  reinforced  composites  are  presently 
being  used  in  a  wide  range  of  engineering 
structures  susceptible  to  some  form  of  high 
speed  or  dynamic  loading.  To  date,  very  few 
systematic  studies  have  been  undertaken  to 
assess  the  rate  sensitive  fracture  properties  of 
fibrous  composites  and  a  detailed 
understanding  of  these  materials  response  to 
rapid  loading  is  still  lacking.  Welsh  and 
Harding  (1)  used  the  split  Hopkinson  bar  test 
to  characterize  the  rate-dependency  of  the 
tensile  strength  and  modulus  of  a  number  of 
long  fiber  reinforced  systems.  Other  workers 
(2,3)  have  used  the  double  cantilever  beam 
test  to  study  the  influence  of  loading  rate  the 
on  interlaminar  fracture  toughness  of 
thermosetting  and  thermoplastic-based 
composite  materials. 

In  this  study  a  fracture  mechanics-type 
specimen  geometry  has  been  used  to  study  the 
influence  of  loading  rate  on  material 
toughness  and  energy  absorption.  Tests  have 
been  undertaken  on  a  wide  range  of  chopped 
strand  mat  composites  enabling  the  influence 
of  matrix  properties  on  fracture  toughness  to 
be  studied  in  detail.  Finally,  the  relationship 
between  the  measured  toughness  of  the 
material  and  it's  resistance  to  impact 
perforation  was  examined. 

The  panels  were  manufactured  by  a  hand  lay¬ 
up  technique  involving  impregnating  three  E- 
glass  fiber  mats  (300g/m^)  with  either  an 
epoxy-based  vinyl  ester  resin  or  a  polyester 
resin.  Twelve  commercially  available  resins 
were  studied  in  order  to  have  the  widest  range 


of  matrix  mechanical  properties.  The  nominal 
thickness  of  the  molded  panels  was  2mm. 
Single  edge  notch  bend  (SENB)  specimens 
having  dimensions  14.5x58mm  x  thickness 
were  then  cut  from  the  panels  using  a 
diamond  slitting  wheel.  Notches  were  then 
machined  at  the  mid-span  using  a  rotating  saw 
blade.  The  notches  were  then  sharpened  by 
sliding  a  fresh  razor  blade  over  the  tip  region. 
The  specimens  were  supported  on  two 
cylindrical  rollers  positioned  58mm  apart  and 
loaded  at  their  mid-points  directly  above  the 
notch.  Tests  using  crosshead  speeds  between 
0.1  and  lOOmm/min  were  conducted  on  an 
Instron  1122  machine  whereas  impact  testing 
was  undertaken  using  an  instrumented  drop- 
weight  tower.  Material  toughness  was 
characterized  by  determining  Kc  and 
measuring  a  work  of  fracture  Wf  by  simply 
dividing  the  energy  dissipated  by  the  fracture 
area. 


Fig.  1.  Damage  zone  in  a  fractured  CSM 
SENB  specimen  tested  at  1  mm/minute. 

A  typical  damage  zone  corresponding  to  a 
fractured  CM  specimen  tested  at  lmm/min.  is 
shown  in  Fig.  1.  It  is  clear  that  energy  has 


been  dissipated  over  quite  a  large  zone  around 
the  principal  crack.  Figs.  2  and  3  show  the 
variation  of  Kc  and  Wf  with  crosshead  speed 
for -an  epoxy-based  vinyl  ester  resin  CSM 
composite. 
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Fig.  2.  Variation  of  Kc  with  rate  for  an  epoxy 
vinyl-ester  based  CSM. 

It  is  interesting  to  note  that  both  of  these 
parameters  indicate  that  the  material  is  capable 
of  absorbing  greater  energy  at  higher  rates  of 
strain.  Indeed,  the  work  of  fracture  increases 
four  fold  over  the  range  examined.  Previous 
work  (1)  has  shown  that  glass  fibers  exhibit  a 
significant  rate  sensitive  response  at  room 
temperature.  Such  effects  were  investigated  in 
more  detail  in  this  study  where  tests  were 
undertaken  on  unidirectional  glass  fiber 
specimen  loaded  in  flexure.  The  results  of 
these  tests  indicated  that  the  tensile  strength  of 
the  fibers  increased  by  almost  thirty  percent 
over  die  range  of  loading  rates  shown  in  Figs. 
1  and  2.  It  is  believed  that  the  trends  apparent 
in  these  figures  result  directly  from  this  rate 
dependent  behavior. 

A  series  of  drop-weight  impact  tests  were 
undertaken  in  order  to  correlate  the  energy 
required  to  perforate  the  panels  with  Kc  and 
Wf.  Fig.  4  shows  the  variation  of  the 
perforation  energy  normalized  with  respect  to 


specimen  thickness  as  a  function  bf  the  work 
of  fracture  measured  under  impact  loading. 
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Log  (Crosshead  Speed)  mm/min 
Fig.  3.  Variation  of  Wf  with  rate  for  an  epoxy 
vinyl  ester  CSM. 

From  the  figure  it  appears  that  the  level  of 
correlation  is  good.  It  is  interesting  to  note 
that  those  composites  having  brittle  matrices 
offer  the  higher  work  of  fractures  and  in  turn 
superior  perforation  energies. 


Work  of  Fracture  (kj/m  ) 


Fig.4.  Variation  of  perforation  energy  with 
work  of  fracture  for  various  composites. 
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The  use  of  embedded  and  surface  mounted 
optical  fiber  sensors  for  health  monitoring  of  key 
structural  components  in  support  vehicle  systems  has 
begun  to  evolve  in  the  last  decade.  These  novel 
sensor  systems  offer  a  new  methodology  to 
investigate  the  occurrence  of  local  damage  and  alert 
maintenance  personnel  to  die  location  and 
progression  of  damage  in  composites  and 
homogeneous  material  systems.  While  die  benefits 
offered  by  health  monitoring  systems  represent 
substantial  cost  saving  in  terms  of  both  money  and 
lives,  a  viable  methodology  to  employ  the  optical 
fiber  sensors  for  damage  detection  and  evaluation 
has  not  been  developed,  a  detailed  study  of  die 
sensors  robustness  to  survive  in  standard  operating 
conditions  has  not  been  performed,  an  adequate 
technique  to  embed  die  sensors  has  not  been 
presented,  and  appropriate  analytical  models 
depicting  the  effect  of  optical  fibers  on  die 
mechanical  properties  of  the  composite  has  not  been 
developed. 

This  investigations  presents  a  new  technique 
employing  Extrinsic  Fabry  Perot  Fiber 
Interferometric  Fiber  Optic  Strain  Sensor  FP-FOSS 
to  measure  thermal  parameters  adequately 
characterizing  damage  mechanisms  arising  in 
today’s  advanced  materials  (e.g.  polymeric,  ceramic, 
and  metal  matrix  composites).  The  monitoring  of 
thermal  parameters  with  FP-FOSS  to  infer  damage 
accumulation  represents  the  first  attempt  to  employ 
this  methodology  in  health  monitoring  applications 
for  composite  materials.  The  information  obtained 
from  the  embedded  and  surface  mounted  FP-FOSS 
on  die  damage  evolution  represents  a  new  approach 
for  determining  the  health  of  a  heterogeneous 
material  system. 


The  approach  employed  in  this  document  is 
founded  upon  sound  scientific  evidence  concerning 
the  dependence  of  a  material’s  fatigue  strength  on 
die  evolution  of  certain  microstructural  damage 
features,  specifically  damage  networking.  The 
sensitivity  of  thermal  parameters  to  this  significant 
damage  feature  is  studied  in  this  investigation  with 
the  use  of  metal  matrix,  ceramic  matrix,  and 
polymer  matrix  composite  materials  subjected  to 
mechanical  and  thermal  fatigue  cycling.  The  results 
show  that  the  proposed  methodology  is  a  significant 
improvement  over  damage  evaluation  techniques  and 
life  predictive  methodologies  that  are  presently 
based  on  number  of  fatigue  cycles  or  stiffness 
reduction  measurements.  With  die  use  of  the 
proposed  characteristic  damage  metrics, 
measurement  of  distributed  microstructural  damage 
features  is  possible  and  is  conelatable  with  loading 
profile,  fatigue  response,  and  residual  strength 
measurements  of  a  composite. 

An  analytical  model  correlating  die  experimental 
data  obtained  from  sensor  measurements  to  die 
damage  state  present  in  die  material  is  discussed. 
This  simple  model  is  presently  based  cm  classical 
shear  lag  techniques.  Test  data  and  analytical  results 
are  presented  for  polymeric  systems  including 
woven  cross-ply  Celion  G30-500/PMR-15  resin 
system  thermally  cycled  from  room  temperature  to 
450  degrees  F  and  a  AS-4-3501/6  quasi-isotropic 
layup  mechanically  loaded.  Additional  test  data  is 
provided  for  a  titanium  matrix  composite  system  and 
a  silicon-carbide  ceramic  matrix  system  subjected  to 
mechanical  loads.  Experimental  and  theoretical 
results  obtained  from  die  polymeric  matrix  system 
reveal  that  the  thermal  parameter  changes  by  as 
much  as  80%  over  the  fatigue  life  of  the  composite 


(see  Figure  1).  This  further  substantiates  the 
supposition  that  these  parameters  can  be  employed 
to  detect  damage  accumulation  in  a  composite  and 
thus  infer  remaining  strength  and  life.  We  show  that 
there  is  a  distinct  correlation  between  damage 
density  and  the  thermal  parameters  in  these  systems, 
both  of  which  are  path  dependent.  That  is,  even 
though  two  specimens  undergo  similar  load  levels 
the  damage  evolution  are  distinctly  different  in  the 
two  specimens.  The  measurement  of  the  thermal 
parameters  delineates  the  damage  present  in  each  of 
the  specimens. 

In  this  paper  we  also  review  various 
manufacturing  techniques  employed  to  construct 
extrinsic  Fabry-Perot  Interferometric  fiber  optic 
sensors  and  correlate  these  processes  with  measured 
strength  values  [1].  We  show  that  if  properly 
constructed  the  fiber  optic  sensors  strain  to  failure 
(i.e.  3%)  exceeds  that  of  classical  polymeric 
composite  materials.  This  results  suggests  that  die 
FP-FOSS  should  fail  after  die  composite  system  has 
foiled,  at  least  in  quasi-static  loading.  Test  results 
are  also  presented  on  fatigue  loaded  polymeric 
composite  systems  and  metal  matrix  composite 
systems  with  attached  FP-FOSS.  These  results 
demonstrate  the  longevity  of  these  sensors  and  the 
advantages  they  offer  over  classical  resistance  strain 
gauges  with  limited  fatigue  endurance  curves. 

We  also  review  die  effect  of  embedded  optical 
fibers  on  a  composite  materials  properties.  It  has 
been  shown  that  the  compression  strength  and  the 
transverse  tensile  strength  of  a  composite  material  is 
compromised  if  the  laminate  is  not  properly 
constructed  [2,3].  We  hove  recently  developed  a 
working  model  providing  detailed- information  on 
appropriate  manufacturing  techniques  of  composites 
containing  optical  fibers.  These  analytical  models 
demonstrate  fiber  optic  coatings  can  be  employed  to 
minimize  the  degradation  in  a  composites  transverse 
tensile  strength  caused  by  optical  fibers.  Guidelines 
involving  the  effect  of  fiber  size  and  orientation  in 
a  laminate  and  die  expected  degradation  on 
compression  strength  is  also  presented. 

The  information  presented  in  this  study  has 
immediate  and  direct  applications  in  both  real  time 
in-situ  health  monitoring  systems  and  non¬ 
destructive  evaluation  techniques  currently  under 
development  The  work  performed  in  this 
investigation  employs  state  of  the  art  measurement 
techniques  to  measure  the  characteristic  damage 
metrics,  including  optical  fiber  sensors.  The 


Figure  1:  Comparison  of  thermal  damage 
metrics. 


knowledge  massed  from  this  investigation 
represents  die  formulation  of  a  new  and  resourceful 
approach  for  accurately  characterizing  die  damage 
state  in  a  heterogenous  material  system  subjected  to 
fatigue  loading.  By  introducing  a  new  technique  for 
quantifying  a  composite  material’s  damage  state,  we 
provide  a  fundamental  approach  not  previously 
embarked  upon  to  unravel  die  complexities 
associated  with  multifarious  damage  phenomena. 
The  results  obtained  from  this  investigation  provide 
essential  information  to  researchers  studying 
composite  materials  and  to  engineers  designing 
structural  components  manufactured  from  composite 
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at  time  t  may  be  written  in  terms  of  the  reference  configuration 
Introduction  Cr  as  follows. 

dHl: 


The  post-failure  mode  is  a  nonlinear  phenomenon  and  is 
important  in  the  design  of  crashworthy  structures.  The 
material  nonlinearities  are  a  consequence  of  plastic  deformation 
whereas  geometric  nonlinearities  are  due  to  large 
displacements,  nonlinearities  in  strain-displacement  relations 
and  changes  in  boundary  conditions  due  to  possible  local 
buckling.  Since  the  response  is  dependent  on  the  material  and 
geometric  properties  of  the  structure,  optimum  selection  of 
design  variables  such  as  ply  orientations,  stacking  sequences 
and  thicknesses  will  yield  efficient  structural  design  for  pre- 
and  post-failure  loading  conditions.  The  objective  of  the 
current  research  is  to  develop  a  sensitivity  analysis  procedure 
which  can  be  used  to  provide  such  necessary  design  trends. 

The  sensitivity  analysis  for  linear  structural  problems  is  a  well 
studied  area.  Research  has  also  been  performed  on  nonlinear 
sensitivities!  1.2].  However,  only  isotropic  material  has  been 
considered.  In  this  paper,  a  semi-analytical  procedure  is 
described  for  the  design  sensitivity  analysis  of  composite 
structures  including  both  material  and  geometric  nonlinearities. 
A  higher  order  approximation  is  used  in  the  integration  of  the 
rate  constitutive  equations  to  improve  accuracy.  A  direct 
differentiation  approach  (DDA)  is  used  to  calculate  the 
sensitivities.  A  partial  differentiation  of  the  constitutive 
equations  in  the  inelastic  range,  is  used  to  calculate  the  partial 
derivative  of  stress  with  respect  to  the  design  variables.  The 
procedure  developed  is  demonstrated  through  a  composite 
laminated  beam  problem. 

General  Theory 

Response  Analysis:  Consider  a  continuum  initially  occupying 

a  domain  Co  of  volume  0  v  and  bounded  by  the  surface  ®p. 
At  load  level  t,  the  configuration  deforms  to  a  new 

configuration  with  a  domain  Ct  of  volume  and  bounded 
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where  the  superscript  t  on  the  left  denotes  time,  C-  is  the 
Cauchy  stress,  Hi:  is  the  displacement  field.  8  denotes  the 
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An  incremental  solution  procedure  is  used  to  solve  the  above 
deformation  problem. 

Constitutive  Models:  The  procedure  being  developed  is  first 
demonstrated  through  a  nonlinear  material  model.  The 
constitutive  equations  are  expressed  as: 
t-  ,  t- 
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are 

the 


where  £j,j  are  the  infinitesimal  strain  rates, 
stiffnesses  which  are  functions  of  the  stresses  at  time  t  and 

design  variables,  that  is,  Cfydj,-  •  -,dn) 

where  dj,-  •  -,dn  are  the  design  variables. 


by  the  surface  *  I". 

A  material  point  =  (0xi,0x2,0x3)Tin  the  initial 

T 

configuration  moves  to  a  new  position  lx  =  (txi,^,^) 
at  time  t.  In  the  configuration  Ct,  any  volume  element  d*V 

experiences  a  body  force  *P  d*V  where  P  is  the  mass 
per  unit  volume  and  b  is  the  body  force  per  unit  mass,  and  any 

oriented  surface  element  d*T  =  d*r  experiences  a  contact 
force  *t(t*»)  dV,  where  H(^i)  is  the  surface  traction  and  4) 
is  the  unit  outward  normal  on  *p.  The  equilibrium  equation 

*  Assistant  Professor,  Senior  Member  AIAA,  Member  AHS, 
ASME,  SPIE 

**  Graduate  Research  Assistant 


Integration  of  Rate  Constitutive  Equations:  The  rate 
constitutive  equations  are  a  set  of  nonlinear  differential 
equations  with  the  elastic-plastic  yield  stresses  as  initial 
conditions.  Solution  of  this  initial  value  problem  for  each 
load  level  is  CPU  intensive.  A  linear  approximation  of  the 
integration  of  the  rate  constitutive  equations  was  used  by  Tsay 
and  Arora  [1].  Here,  a  higher  order  approximation  of  the 
integration  is  proposed  to  increase  the  accuracy.  Integration  of 
Eqn.(2)  yields  the  following. 


This  can  be  approximated  as  follows. 

t  X  X  X  . 

°.J=  Oy  +  xQijkl  ( ek!  -  ekl) 

+  J  (dXQijmn^Tekl)  <£kl  ~  t£klXtemn  -  \n)  (3) 
where3TQijmn^Xekl  =  (XQijmn^X°pq) 


which  can  be  expressed  in  terms  of  the  stresses  at  time  X  and 
the  design  variables. 

Sensitivity  Analysis 


Direct  Differentiation  Approach:  Direct  differentiation  of  the 
equilibrium  equation  with  respect  to  the  design  variable  dm 
yields 
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In  Eqn.(4).  <JJ,  rJ,  <JJp,  rJp,  rP,  Xbj  and  are  explicit 
functions  of  the  design  variables.  Their  sensitivities  with 
respect  to  the  design  variable  dm  are  easy  to  determine.  Since 
the  responses  are  known  at  time  t,  the  quantities  dtGij^fepq 
are  also  known.  The  quantities  dt£pq/ddn]  can  be  obtained 
from  dtui/ddm.  Therefore  the  unknown  quantities  are 
akjy/ad,,,  and  dkij^dn,.  The  following  section  describes  a 
procedure  to  solve  for  dtOiyddm.  Once  dtGiyddm  are 
known,  dttiiAHn  can  be  solved  using  Eqn.(4). 


Partial  Differentiation  of  Constitutive  Equations:  Fixing  the 
state  variables,  partial  differentiation  of  the  constitutive 
equations  (2)  with  respect  to  the  design  variable  (dm)  yields 


dO‘°,/3dJ  -  O’Qyu Aw)  A.I  (5) 

This  is  a  set  of  linear  differential  equations  in  9  Gt7ddm. 

Their  initial  values  can  be  known  by  differentiating  the 
stresses  at  the  yield  point  with  respect  to  the  design  variable. 
This  linear  problem  can  easily  be  solved  numerically. 

Numerical  Example 

The  above  procedure  is  applied  to  a  (0*  / 10*  /  45’  )2  anti 

symmetric  Glass/Epoxy  composite  beam  with  equal  ply 
thicknesses  (0.2in)  undergoing  pure  bending.  The  ply  angles 
are  chosen  as  the  design  variables.  The  Tsai-Hill  theory  is 
adopted  as  the  yield  criterion.  In  the  rest  *>nse  analysis,  the 
results  of  the  second  order  approximation  of  the  integration  of 
the  rate  constitutive  equations  are  compared  with  the  direct 
solution  of  the  set  of  nonlinear  differential  equations  (Table  1). 
They  show  excellent  agreement  in  both  elastic  and  plastic 
ranges.  However,  the  CPU  time  requirement  for  the  second 
order  approximation  is  much  less  than  that  for  the  differential 
solution.  In  the  sensitivity  analysis,  the  results  of  the  direct 
differentiation  approach  (DDA)  are  compared  with  the  those 
obtained  using  the  finite  difference  approach  (FDA).  The 

sensitivities  of  curvature  with  respect  to  the  ply  angle  63  (20 
deg.)  are  presented  in  Table  1  and  show  excellent  agreements  in 
both  elastic  and  plastic  ranges. 

Table  1  Comparison  of  response  and  sensitivity  analysis 
Load _ , _ Curvature _ Sensitivity  of  curvature 


xlO^Ib-in  Direct  Sol  2nd-orSolu 

FDA 

DDA 

0.1 

.001309 

.001309 

.000922 

.000922 

0.2 

.002618 

.002618 

.001845 

.001845 

0.3 

.003947 

.003968 

.003458 

.003523 

0.4 

.005584 

.005652 

.007279 

.007436 

0.5 

.007377 

.007453 

.010068 

.010238 

0.6 

.009205 

.009290 

.012938 

.013112 

0.7 

.011062 

.011153 

.015842 

.016013 

0.8 

.012940 

.013038 

.018816 

.018991 

0.9 

.014846 

.014951 

.021935 

.022106 

1.0 

.016775 

.016886 

.025047 

.025212 
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VIBRATION  ANALYSIS  OF  FINITE  THIN 
LAMINATED  DOUBLY-CURVED  PANELS 
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ABSTRACT 


Governing  partial  differential  equations 
and  the  associated  boundary  conditions, 
for  the  boundary  -  value  problem  of  large 
amplitude  vibration  of  arbitrarily 
laminated  (of  which  symmetric/ 
antisymmetric  cross-  ply  and  angle-ply 
laminations  are  special  cases)  thin 
doubly-curved  panels,  that  account  for 
von-Karman  type  geometric  nonlinearity, 
are  derived,  starting  from  the  method  of 
virtual  work.  Two  numerical  solution 
approaches,  based  on  finite  element  and 
global  Galerkin  methods,  are  also 
presented.  The  latter  is  restricted  to 
curved  panels  of  rectangular  planform. 
Actual  computation  of  the  numerical 
results  is  currently  underway  at  the 
University  of  Utah,  and  will  be  presented 
in  future.  Next,  the  governing  partial 
differential  equations  and  the  associted 
boundary  conditions  are  linearized  for 
the  special  case  of  small  amplitude 
vibration,  which  are  solved  using  two 
recently  developed  Fourier  type 
approaches.  Details  of  these 
solutions  to  the  boundary-value  problem 
of  free  small  amplitude  vibration  of 
arbitrarily  laminated  thin  doubly-curved 
panels  of  rectangular  planform  are 
available  elsewhere.  Both  boundary- 


discontinuous!  1]  and  boundary- 
continuous-  displacement  [2]  based 
double  Fourier  series  approaches,  are 
employed  to  solve  the  linear 
boundary-value  problems,  involving 
highly  coupled  linear  partial  differential 
equations  with  constant  coefficients, 
resulting  from  CLT  (Classical  Lamination 
Theory)-based  formulations  that  also 
include  surface-parallel  inertias. 

Extensive  numerical  results  that  are 
presented  in  this  study  include  (i) 
convergence  characteristics  of  computed 
natural  frequencies,  and  (ii)  effects  of 
length-to- thickness  ratio,  radius-to-length 
ratio,  fiber  orientation  angle,  lamination 
sequence,  shell  geometry  and  boundary 
constraints  on  the  response  quantities  of 
interest  The  accuracy  of  the  solutions 
are  ascertained  by  comparison  with  the 
available  FSDT-based  analytical  and 
CLT-based  Galerkin  solutions. 
Comparisons  with  the  available  FSDT 
(first-order  shear  deformation  theory) 
-based  analytical  solutions  help  in 
establishing  the  upper  limit  (with  respect 
to  the  thickness-to-length  ratio)  of 
validity  of  the  present  CLT  (classical 
lamination  theory)- based  solutions.  Also 
investigated  is  the  highly  complex 


interaction  among  bending-stretching  type 
coupling  effect,  membrane  action  due  to 
shell  curvature,  and  the  effects  of 
surface-parallel  (membrane)  inertias.  It  is 
also  concluded  that  boundary- 
continuous-displacement  based  Fourier 
solutions  are  superior  to  their 
boundary-discontinuous  counterparts  in 
their  ability  to  avoid  Gibb's  phenomenon, 
in  the  case  of  rigidly  clamped  boundary 
conditions. 
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Fig.  2  Variation  of  normalized  fundamental  natu¬ 
ral  frequencies  of  square  relatively  deep  l  ft /a  * 
5.73)  and  relatively  thin  (ath  »20j  antisymmetric 
angle-ply  two-  and  four-layer  spherical  panels  with 
0.  and  comparison  with  the  FSDT. 


Fla  1.  Variation  of  normalized  fundamental  and  second  (numerically  ordered)  natural  frequencies  of 
square  antisymmetric  angle -ply  rwo^ayer  (457-45*)  and  four-layer  *Phe***1 
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INTRODUCTION 

The  use  of  electrorheological  (ER)  materials 
in  adaprive  structures  has  received  much  attention 
recently.  Electrorheological  materials  experience 
reversible  changes  in  rheological  properties  such  as 
viscosity,  plasticity  and  elasticity  when  subjected  to 
electrical  fields.  These  reversible  changes  are  due  to 
chain  formation  of  micron-sized  dielectric  particles  in 
non  polar  media  under  an  applied  electrical  field, 
resulting  in  up  to  several  orders  of  magnitude  increases 
in  rheological  properties,  hi  essence,  ER  material 
behavior  transforms  from  dun  of  a  liquid  to  dial 
characteristic  of  a  solid-like  gel  when  an  electric  field 
is  applied.  Adaptive  structures  are  based  on  controlling 
the  pre-yield  theology  of  ER  materials. 

In  this  study  the  dynamic  behavior  of  an  ER 
based  adaptive  beam  was  modeled,  based  on  shear 
configurations.  The  beam  was  composed  of  three 
layers:  an  ER  material  controllable  damping  layer  and 
surrounding  upper  and  lower  elastic  plates.  The 
existing  models  for  viscoelastically  damped  composite 
three  layer  beams  were  examined  for  their 
applicability  to  ER  adaprive  beams.  From  the 
rheological  findings  it  was  observed  that  ER  materials 
exhibit  linear  rheological  behavior  similar  to  many 
common  viscoelastic  materials.  Based  on  this 
similarity  in  material  behavior  a  modification  of 
existing  models  used  for  viscoelastically  damped 
sandwich  structures  was  proposed.  . 

The  structural  model  of  the  assembly  in  a 
transverse  continuous  vibration  mode  subjected  to 
simply-supported  boundary  conditions  and  actuation  at 
a  single  point  on  die  adaptive  beam  surface  was 
analyzed.  Theoretical  structural  natural  frequencies, 
and  mode  shapes  for  a  0  -300  Hz.  bandwidth  range 
were  calculated  and  the  transverse  displacement  of  a 
specified  point  onanER  material  adaptive  beam  of  the 
experimental  setup  is  predicted.  The  earlier 
experimental  results  of  Don  and  Coulter  (Don,  1993) 


were  used  for  comparison  purposes.  The  analytical 
results  are  compared  with  experimental  results  under 
the  same  physical  conditions.  Qualitative  agreement 
between  theory  and  experimentation  resulted.  In 
addition  an  effort  was  made  to  reduce  the  vibration  of 
the  structure  by  selecting  the  optimum  electrical  field 
which  yields  minimized  vibration  for  each  frequency. 

MODEL  DESCRIPTION 

In  the  Ross-Kerwin-Ungar  (RKU)  analysis, 
the  transverse  vibration  equation  derived  for  Euler- 
Bemoulli  beams  is  considered.  The  total  effective 
flexural  rigidity  is  found  by  assuming  that  all  of  the 
layers  are  composed  of  ideal  elastic  materials.  Once 
the  expression  is  obtained  for  the  transverse  vibration 
of  the  beam,  the  real  modulus  in  die  governing 
vibration  equation  is  replaced  by  a  complex  modulus 
for  the  viscoelastic  layer  only.  The  final  form  of  the 
flexural  rigidity  is  expressed  in  the  form  of  real  and 
imaginary  components.  A  detailed  derivation  of  the 
complex  flexural  modulus  can  be  found  in  the 
literature  (Ross  et  al.,  1959). 

The  basic  model  developed  by  RKU  can  be 
further  extended  by  including  boundary  conditions,  an 
external  forcing  function,  and  the  controllable  nature 
of  ER  materials.  Structural  adaptability  is 
incorporated  by  assuming  that  the  complex  shear 
modulus  of  the  ER  material  is  a  function  of  electrical 
field.  The  Young’s  modulus  of  the  ER  material  is 
considered  negligibly  small  compared  to  Young's 
moduli  of  elastic  plates. 

The  transverse  displacement  of  die  adaptive 
beam  was  obtained  by  die  use  of  the  expansion 
theorem.  Considering  the  orthogonality  condition  of 
modes,  the  differential  equation  of  the  continous 
system  lead  to  an  infinite  set  of  uncoupled  ordinary 
differential  equations  which  is  in  summation  resulted 
with  die  total  transverse  displacement 

Thus,  die  extended  model  is  capable  of 
determining  natural  frequencies,  mode  shapes  and 
transverse  vibrational  response  at  a  point  on  the 
adaptive  beam  for  an  applied  point  force  at  a  specified 
location. 

RESULTS  AND  DISCUSSIONS 

Theoretical  transverse  vibration  response 
predictions  were  obtained  and  recorded  for  varying 
electric  field  levels  from  0  -  3.5  kV/mm  with  0.1 
kV/tnm  increments,  with  a  frequency  range  of  0  -  300 
Hz.  The  transverse  vibration  amplitude  as  a  function 
of  frequency  for  electrical  fields  of  0  kV/mm,  1.5 


kV/rrnn,  2.5  kV/mm,  and  3.5  kV/mm  is  plotted  in 
Figure  1.  In  this  figure  the  peak  values  of  each  curve 
represent  the  resonance  frequencies  of  the  structure 
corresponding  to  that  specified  electrical  field.  As  the 
electric  field  increases,  the  resonance  frequencies  shift 
to  higher  values  while  decreasing  the  general 
transverse  response  of  the  structure.  Corresponding 
experimental  results  shows  the  same  trend  of  vibration 
response.  Also  the  experimental  results  validate  the 
theoretically  predicted  controllability  of  the  structure. 

The  final  step  taken  during  the  present  study 
was  to  predict  electrical  field  levels  for  optimal 
frequency  dependent  control  of  the  prototype  ER 
material  based  adaptive  beam.  For  this  purpose  the 
electrical  field  applied  to  the  HR  material  damping 
layer  which  caused  the  smallest  transverse  vibration 
response  of  the  beam  was  chosen  as  optimal.  At  each 
forcing  frequency,  only  oat  electric  field  level  was 
optimal.  Theoretically  predicted  and  experimentally 
observed  optimal  electrical  field  levels  are  plotted  as  a 
function  of  forcing  frequency.  Figure  2  illustrates  the 
theoretical  predictions.  In  both  cases  the  optimal 
electrical  field  pattern  resembled  a  repeating  ramp 
type  curve  with  a  visible  period  of  frequency.  If 
electrical  fields  higher  than  3.5  kV/mm  were 
permitted,  the  upper  limit  would  be  expected  to  have 
sharp  peaks  rather  than  the  flat  plateaus  shown.  The 
only  deviation  between  the  theoretically  predicted  and 
experimentally  observed  optimal  electric  field  levels 
was  die  period  of  this  repeating  pattern.  This 
deviation  could  be  explained  by  inaccurate  theoretical 
ER  material  rheological  properties. 
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CONCLUSIONS 

In  this  study,  the  theoretical  modeling  and 
experimental  analysis  of  ER  material  adaptive 
structures  based  on  shear  configurations  was 
performed.  A  model  was  proposed  and  tested  under 
varying  farcing  frequency  and  applied  electrical  field 
levels.  The  results  were  compared  with  experimental 
tests  with  an  actual  ER  material  adaptive  beam.  In 
both  cases  modal  natural  frequencies  increased  and 
modal  vibration  amplitude  levels  changed  as  the 
electric  field  applied  to  the  ER  material  was  increased. 
This  controllable  characteristic  of  ER  material 
adaptive  structures  is  useful  when  variable 
performance  of  the  structure  is  desired. 
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Figure  1  Theoretical  frequency  response  of  adaptive  beam 
at  Electric  Fields  of  0, 1J,  25,  and  35  kV/mm. 
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Figure  2  Theoretical  Optimal  Electric  Field  Corresponding 
Each  Excitation  Frequency. 
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Aircraft  vulnerability  codes  rely  upon  several  critical 
parameters  for  accuracy  of  the  codes  predictions.  One  such 
parameter  is  the  ballistic  limit  (also  referred  to  as  the  Vw). 
Probabilistic  in  nature,  the  VM  describes  a  target's  penetration 
resistance  and  is  defined  as  the  impact  velocity  where 
penetration  will  occur  exactly  50%  of  the  time.  The  V„  is 
commonly  determined  experimentally  (or  estimated  based  on  an 
empirical  equation)  due  to  the  parameter's  variation  with 
changes  in  projectile  or  target  geometry1-1*4.  Penetrated 
composite  laminates  have  been  observed  to  eject  massive 
amounts  of  spall  particles,  often  numbering  in  the  thousands1. 
Composite  spall  are  assumed  to  exit  the  laminate  over  a  range 
of  velocities  based  on  the  inspector's  velocity  and  the  distance 
of  the  spall's  origin  from  the  shotline.  When  the  energy 
associated  with  spallation  becomes  significant,  simplistic 
equations  (such  as  those  to  predict  the  ballistic  limit)  become 
confounded. 

Experimental  data  have  shown  that  a  nearly  linear 
relationship  exists  between  a  projectile’s  initial  (pre-impact)  and 
residual  (post-impact)  kinetic  energies*7.  Deviations  from 
linearity  are  the  result  of  two  primary  factors.  During  the 
penetration  event,  energy  is  absorbed  within  the  laminate  in  the 
form  of  delammadoa,  fiber  fracture,  matrix  cracking,  and  fiber- 
retin  interface  failure.  Energy  absorbed  during  penetration  can 
also  be  expected  to  include  target  displacement,  heat,  noise,  and 
even  projectile  deformation  (in  the  case  of  semi-rigid 
impactors).  Upon  completion  of  the  penetration  event, 
continued  deviation  from  linearity  is  believed  the  result  of 
spallation.  Although  energy  absorbed  during  penetration  can 
easily  be  assessed  and  equated  to  foe  E*  (defined  as  foe  kinetic 
energy  absorbed  at  the  V„),  the  energy  associated  with 
■paiUfi^i  icrnaim  questionable. 

The  objective  of  this  study  is  to  determine  mass  and  velocity 
relationships  between  spherical  metallic  penetratare  and  spall 
discharged  from  laminated  composites.  Of  particular  interest  is 
defining  foe  critical  impact  energy  requirement  to  achieve  a 
threshold  level  of  spallation.  Also  of  value  is  foe  initial  rate  of 
spall  mass  generation  and  how  latent  reductions  in  this  rate 
affect  the  average  spall  velocity. 

The  target  material  used  in  this  series  of  experiments  was 
grephite/epoxy  (AS4/3501-6).  All  experiments  were  performed 
on  32-ply  [(0/90lM5/-43)(]t  laminates  having  a  total  thickness 
of  0.180- inch.  The  projectile’s  mass  was  adjusted  using  1/2- 
inch  diameter  aluminum,  steel,  and  tungsten  spheres 
(6.437x10*  18.919x10*.  and  35.231x10*  Ibtn,  respectively). 
All  shotlines  were  normal  to  foe  laminate’s  surface  with 
impacts  occurring  in  the  geometric  center  of  each  specimen. 
Projectile  launch  velocities  were  tightly  controlled  using  a  light- 
gas  gun.  Several  tests  were  performed  with  each  sphere  mass 
at  impact  velocities  near  the  ballistic  limit  to  establish  the  on¬ 
set  of  spallation  as  well  ax  a  precise  Vn  velocity.  Higher 
velocity  impact  testa  were  conducted  to  determine  the  change 


in  spall  mass  and  velocity  as  a  function  of  impact  velocity. 

As  shown  in  Figure  1,  foe  relationship  between  the 
projectile's  initial  (EJ  and  residual  (E,)  energies  remains 
roughly  linear  over  foe  majority  of  foe  energy  spectrum  of 
interest  Although  foe  slopes  vary  slightly  for  data  generated  by 
projectiles  of  differing  masses,  all  data  roughly  converges  on  a 
single  energy  required  for  penetration  (between  37 J  and  43.5 
ft-Ibs).  This  kinetic  energy  requirement  for  penetration 
(referred  to  as  the  Ex)  is  i*~%ely  independent  of  the  projectile’s 
mass  and  can  be  considered  constant  for  the  laminate  on  which 
it  is  established.  Other  parameters  such  as  the  projectile's 
geometry  (note  shape  and  cross  sectional  area)  and  target 
specific  parameters  (material  and  lay-up)  are  expected  to 
influence  foe  E».  but  are  outside  foe  scope  of  this  paper. 

For  steel  and  tungsten  impactors,  the  E,  vs.  E,  slope  changes 
from  1.00  to  0.82.  This  change  is  hypothesized  as  resulting 
from  the  influence  of  spallation.  A  differing  slope  transition  for 
the  softer  aluminum  impactor  (to  0.63)  is  presumed  due  in  part 
to  the  projectile’s  deformation. 

To  maintain  conservation  of  energy,  the  projectile’s  initial 
energy  must  be  perfectly  transformed  into  several  components 
upon  impact  Data  deviating  from  the  theoretical  1:1  slope  are 
the  result  of  many  factors,  all  of  which  produce  a  net  effect  of 
reducing  foe  projectile’s  residual  velocity.  A  designated  amount 
of  energy  (equivalent  to  the  E*)  is  absorbed  by  foe  laminate 
during  foe  perforation  process.  After  surpassing  the  penetration 
threshold,  only  two  major  components  of  kinetic  energy 
remain._foat  associated  with  foe  projectile’s  residual  velocity 
and  with  spall  particles  ejected.  Any  change  of  impact  energy 
above  foe  E»  must  therefore  equal  the  change  of  foe  projectile’s 
residual  energy  plus  foe  spall  energy.  By  setting  the  x-intercept 
equal  to  foe  E*.  foe  1:1  relationship  between  energy  "in"  and 


Figure  1.  Deviations  in  foe  E,  vs.  E*  slope  as  a  function  of 
kinetic  energy  posessed  by  aluminum,  steel,  and  tungsten 
impactors. 


energy  "out"  becomes  a  useful  tool  for  estimating  spall  energy. 
The  energy  associated  with  spallation  (E,)  can  therefore  be 
calculated  according  to  the  equation  E^EpE^-E*,. 

Spallation  energy  is  plotted  as  a  function  of  impact  energy 
in  Figure  2.  As  impact  energy  increases,  the  spall  energy 
becomes  a  significant  factor  capable  of  diversely  infWring 
VM  predictions. 

Figure  3  describes  the  generation  of  spall  mass  as  a  function 
of  impact  velocity.  For  the  rigid  steel  and  tungsten  impactors, 
the  critical  velocity  far  spall  initiation  coincides  exactly  with 
the  V*.  Only  for  the  softer  aluminum  peneCrator  does 
spallation  begin  at  a  lower  velocity  (0.87  times  the  V*).  [Note: 
The  ciive  shift  associated  with  the  change  in  velocity  between 
0.87  and  1.00  rimes  the  V„  (corresponding  to  93  ft-lbs)  is 
attributed  to  the  projectile's  deformation  energy.]  Although 
spall  initiation  coincides  roughly  with  the  V,,,  spallation  is  not 
maximized  until  velocities  are  in  excess  of  2.8  times  the 
As  impact  velocities  increase,  the  slope  transitions  to  zero  as 
the  mass  subject  to  spallation  attains  a  maximum  limit.  For 
impact  velocities  below  the  transition,  shear  plugging  is 
incomplete,  whereas  for  velocities  above  the  transition,  an  open 
hole  (roughly  the  diameter  of  the  projectile)  is  obtained. 


Impact  Energy  (tt-lbs) 


Figure  2.  Spallation  energy  as  a  function  of  impact  energy 
using  aluminum,  steel,  and  tungsten  impactorr. 
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Figure  3.  Spall  mass  generation  as  a  function  of  the  normalized 
impact  velocity  using  aluminum,  steel,  and  tungsten  impactors. 


Although  not  directly  measurable,  the  average  spall  velocity 
(V|)  can  be  estimated  according  to  the  equation 

VI-K2XE,Xm,-,)r 

where  m,  is  the  spall  mass  (equal  to  the  change  in  target 
weight).  Figure  4  describes  the  average  spallation  velocity  and 
projectile  residual  velocity  as  a  function  of  the  impact  velocity. 
The  spall  velocity  is  dependent  on  the  projectile's  mass  and 
quickly  becomes  large  for  impact  velocities  above  the  V^. 
Note  that  over  the  entire  impact  velocity  spectrum,  composite 
spall  generated  fay  the  aluminum  projectile  has  a  velocity 
identical  to  that  of  the  projectile's  residual  velocity.  For  the 
heavier  impactors  (having  densities  significantly  greater  than 
that  of  the  composite  spall),  the  spall  velocity  deviates 
significantly  from  the  projectile's  residual  velocity.  Differences 
between  the  spall  and  projectile  velocity  curves  sre  seen  to 
increase  with  projectile  mass.  Low-mass  high-velocity 
projectiles  produce  a  greater  spall  velocity  than  high-mass  low- 
velocity  projectiles. 


Impact  Velocity  (Us)  (normalized  to  the  VSO) 


Figured.  Average  spall  velocity  and  projectile  residual  velocity 
as  functions  of  normalized  impact  velocity  using  aluminum, 
steel,  and  tungsten  impactors. 
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The  behavior  of  pretwisted  rods  is  distinguished  by  the 
ever  presence  of  two  coupling  effects:  (1)  extension  with 
twist  and  (2)  bending  about  one  axis  with  that  of  the  other 
(even  when  principal  axes  are  adopted).  This  paper  is 
concerned  with  a  finite  element  method  of  analysis  of  the 
natural  vibrations  in  pretwisted  rods  based  on  on  three- 
dimensional  elasticity.  The  finite  element  modeling 
occurs  over  the  cross-section,  which  accommodates  a 
cross-section  with  arbitrary  shape,  inhomogeneity, 
mechanical  and  inertial  properties  and  location  of  the 
pretwist  axis.  The  relevant  equations  from  linear  elasti¬ 
city  are  those  in  terms  of  a  rotating  coordinate  system 
that  tracks  the  cross-sectional  geometry  along  the  axis  of 
twist  The  cross-sectional  finite  element  modeling  leaves 
the  axial  dependence  £  and  time  t  undetermined  at  the 
outset. 


The  equations  of  motion  derived  by  Hamilton’s 
principle  are  partial  differential  equations  in  the  axial 
coordinate  and  time. 


Only  real  k's  are  considered  in  this  paper,  and  they 
represent  propagating  waves  and/or  harmonic  standing 
vibrations.  Complex  eigenproblem  (2)  can  be  rendered 
completely  real  by  combining  it  with  its  multiplication  by 
-i  to  give 


K3  Hr  kZK,  kK2 
-kK2  K3  +  k2K, 


UD 


-iU„ 


M 

u„  ‘ 

=  (I)2 

M 

-iU„ 

(3) 


From  symmetry  of  [K,]  and  [K3]  and  antisymmetry  of 
[K2J.  algebraic  eigensystem  (3)  involves  real,  symmetric, 
positive  definite  matrices,  so  that  all  cd2  are  real  and  posi¬ 
tive.  Further  study  of  Eq.  (3)  shows  it  to  consist  of  two 
identical  subsystems,  yielding  identical  solution  pairs. 
They  represent  the  same  wave  form  except  for  a  rt/2 
phase  difference  in  the  axial  direction. 


[Ki]{  U"  }  +  [KjlIU'l  -  [KjHU| 

=  ( M 1  {  U  |  (1) 

where  primes  and  dots  denote  differentiation  with  respect 
to  £  and  t.  respectively.  The  stiffness  matrices  [Kj],  [K3] 
and  mass  matrix  [M]  are  symmetric  and  [K2]  is  antisym¬ 
metric. 


A  computer  code  prepared  was  prepared 
employing  8-node  quads  and  6-node  triangles  (both  both 
quadratic  interpolations).  Isoparametric  finite  element 
methodology  was  used  in  element  matrix  formations. 
Because  of  the  large  size  of  the  algebraic  eigensystem. 
subspace  iteration  was  used  to  extract  the  lowest  subset 
of  the  eigendata. 


Introducing  a  harmonic  wave  form,  U  =  U„ 
exp{i(k£  +  wt)|  into  Eq.  (1).  where  k  =  n/X  is  the  axial 
wave  number  (X  the  wave  length)  and  &>  the  natural  circu¬ 
lar  frequency,  results  in  an  algebraic  eigensystem  in 
terms  of  k  and  oi  in  the  form 

([Kj]  -  ik[K2]  +  k2 [K,] )  (U.J 


Two  examples  were  given  in  the  full  length 
paper  to  illustrate  the  method  of  analysis  and  vibrational 
characteristics,  viz.,  (1)  a  homogeneous,  isotropic  bar  and 
(2)  a  two-layer  ±30°  angle  ply  composite,  both  cross- 
sections  with  the  aspect  ratio  of  ten  (width  b  to  height  h). 
Herein,  only  results  for  the  two-layer  composite  are 
presented  (see  Fig.  1). 


(2) 


The  properties  of  the  composite  are 

^*0.5:  ~~  =  0.38 ; 

Br  Br  Et 

(4) 

Vlt  *  0.21 ;  v-jt  ■  0.3 

All  frequencies  o>  were  normalized  as  follows: 
fl2  *  (i^Er/ph2 

where  p  is  the  unit  mass  density.  Five  cases  of  pretwist 
were  considered  for  both  bars:  (1)  a  *  0.3°  about  the 


center.  (2)  o«0.5°  about  the  quarter  point:  (3)  a*  l.(f 
about  the  center.  (4)  a  =  1.0°  about  the  quarter  point:  (3) 
no  pretwist.  The  pretwist  angle  a  refers  to  a  twist  rate 
relative  to  the  thickness  h  of  the  bars. 

Shown  in  Fig.  2  are  the  spectral  curves  for  the 
lowest  six  modes  of  vibrations  for  cases  (1)  and  (3).  In 
Fig.  3  is  shown  the  cross-sectional  modal  patterns  of  the 
2nd  lowest  (torsion-extensional)  mode  of  vibrauon. 
These  plots  reveal  the  complicated  but  very  intriguing 
nature  of  the  behavior  in  pretwisted  rods. 


The  uial  and  inpUne  displacement]  may  be  of  different  scales 
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ABSTRACT 


The  development  of  constitutive  models  for 
composite  materials  is  motivated  in  large  pan 
by  their  use  in  computer-aided  structural 
analysis  and  design.  Finite  element  methods 
are  by  far  the  most  commonly  applied  tool  in 
such  applications.  Due  to  the  large  CPU  time 
and  memory  requirements  of  dynamic,  three- 
dimensional  finite  element  analysis,  most 
micromechanical  composite  material  models 
see  numerical  implementation  in  the  form  of 
anisotropic  continuum  models  with  internal 
state  variables  (Lagoudas  and  Huang,  1992). 
The  latter  farm  is  both  compatible  with 
general  purpose  finite  element  codes  and 
acceptably  computer  resource  intensive  for 
large  scale  structural  simulations.  Typically 
first  order  evolution  equations  for  the  internal 
state  variables  are  appended  to  a 
displacement-based  finite  element  model  in 
order  to  simulate  the  composite  structural 
dynamics. 

Although  displacement-based  finite  element 
methods  are  well  suited  to  isothermal  solid 
dynamics  modeling,  many  composites 
engineering  applications  place  a  particular 
emphasis  on  thermomechanical  coupling 
effects.  Included  in  this  category  are  studies 
of  fabrication  processes  (Pindera  and  Freed, 
1992),  the  simulation  of  high  velocity  impact 
dynamics  (Drumheller,  1987),  and  design 
with  shape  memory  alloys  (Gandhi  and 
Thompson,  1992).  For  mixed  energy  domain 
problems  of  the  thermomechanical, 
thermochemical-mechanical,  or 
thermoelectric-mechanical  type,  diffusion 
models  or  other  evolution  equations  are 
normally  appended  to  the  finite  element 
structural  model,  in  order  to  simulate  die 
overall  system  dynamics.  Although  this 
method  has  met  with  some  success,  the 


disadvantages  of  a  somewhat  ad  hoc  or 
problem  dependent  modeling  approach 
become  apparent  as  the  complexity  of  the 
problems  of  interest  increases.  Clearly  a 
unified  and  general  purpose  approach  to  the 
development  of  numerical  models  for  coupled 
solid  dynamics  and  diffusion  problems  is  of 
interest  for  a  significant  class  of  composites 
engineering  problems. 


Fig.  1  Bond  graph  element  representing 
internal  energy  storage  in  an  anisotropic 
thermoelastic  material  with  continuum 
damage. 


The  bond  graph  approach  (Rosenberg  and 
Kamopp,  1983)  to  system  dynamics 
modeling  offers  a  unified,  energy  based 
method  for  the  formulation  of  state  space 
models  of  nonlinear  systems,  incorporating 
complex  coupling  of  mechanical,  thermal, 
electrical,  chemical,  and  other  energy 
domains.  Although  bond  graphs  have  seen 
extensive  application  to  low  order  lumped 
parameter  systems,  only  recently  has  the 
method  been  extended  to  continuum  solid 


dynamics  of  the  elastic-plastic  (Fahrenthold 
and  Wu,  1988)  or  thermoelasdc  (Ingrim  and 
Masada,  1988)  type.  No  previous  bond 
graph  work  has  provided  a  general  modeling 
approach  to  the  anisotropic,  thermoelasdc, 
solid  dynamics  and  heat  diffusion  problem 
with  internal  state  variables. 


Fig.  2  Bond  graph  element  representing 
kinetic  energy  storage  in  a  tetrahedral  finite 
element. 
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Fig.  3  Bond  graph  element  representing 
thermomechanical  damage  evolution  in  a 
tetrahedral  finite  element. 


Recognizing  the  importance  of  the  latter  class 
of  problems  to  composite  structural  analysis 
and  design,  the  present  paper  develops  a 
three-dimensional,  state  space  description  of 
coupled  heat  diffusion  and  structural 
dynamics,  for  the  case  of  anisotropic 
thermoelasdc  materials  with  continuum 
damage.  Combining  bond  graph  techniques 
for  mixed  energy  domain  modeling  with  well 
established  finite  element  discretization 
methods  (Fahrenthold  and  Wargo,  1992) 
provides  a  systematic  and  thermodynamically 


rigorous  treatment  of  the  thermomechanical 
problem.  The  method  is  well  adapted  for 
extension  to  include  electrical,  chemical,  and 
other  effects  of  particular  interest  to 
composites  engineers. 
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Abstract 

Thick  components  lt>  structures  such  as  rocket 
aotor  casos,  prsssurs  vassals  and  hulls  for  tanks, 
ships  and  subearlnas  arc  often  subjected  to  vibratory 
loading,  but  little  Is  known  of  the  vibration  deeping 
characteristics  of  such  coaponants.  This  Is 
partlculsrly  true  of  thick  coaponants  made  of 
coeposlte  Materials. 

In  general,  conventional  Metallic  structures 
have  poor  Internal  doping,  and  Most  of  the  daaplng  in 
such  structures  comss  frou  joint  friction.  The 
enhancanent  of  vibration  doping  characteristics  of 
thin  Metallic  panels  Is  typically  accoapllshed  by  the 
application  of  constrained  viscoelastic  layer  daaplng 
treataenta  to  the  panel  surfaces  (Nashlf,  at  el., 
1985).  Such  daaplng  treataenta  have  proved  to  be  very 
effective  In  Inproving  daaplng,  but  at  the  expense  of. 
added  weight  due  to  the  daaplng  layer  and  the  netal 
constraining  layer.  Host  of  the  literature  on 
constrained  layer  daaplng  deale  with  thin  beeas  and 
plates,  however,  and  the  effects  of  thick  base 
structures  have  not  been  investigated.  Uhtle  advanced 
polyaer  conposlte  materials  offer  the  potential  for 
the  developaent  of  strong,  lightweight  and  highly 
daaped  structures,  previous  research  has  again  dealt 
aalnly  with  thin  structures  with  and  without  surface 
daaplng  treataenta  (Han tens,  at  al. ,  1991)  and  not 
auch  Is  known  of  the  daaplng  characteristics  of  thane 
materials  whan  used  In  thick  structures. 

This  paper  presents  the  preliminary  results  of 
an  investigation  of  the  vibration  characteristics  of 
thick  cantilever  beams .  The  principal  objective  Is  to 
study  the  effects  of  Increased  beaa  thickness  on 
daaplng  and  resonant  frequency  characteristics  of 
thick  structural  eleaents  of  two  types:  (1)  alualnua 
beaa  structures  with  constrained  viscoelastic  layer 
dealing  treataenta,  and  (2)  graphite  fiber -reinforced 
epoxy  laminated  beaa  structures  without  daaplng 
treatments . 

The  analysis  of  daaplng  was  carried  out  by  using 
a  strain  energy/f inite  clcnent  approach,  which  is 
often  referred  to  as  the  aodal  strain  energy  approach. 
The  theoretical  basis  of  this  approach  is  the  Ungar- 
Kervln  equation  (Ungar  and  Kervln,  1962)  for  the  total 
structural  loss  factor  In  tens  of  eleaent  loss 
factors  and  the  fraction  of  strain  energy  scored  in 
each  eleaent: 
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where  m  Is  the  loss  factor  of  the  1th  eleaent,  V*  Is 
the  strain  energy  stored  in  the  1th  element  at  aaxlaua 
vibratory  dlsplaeeaent  and  n  is  the  number  of 
eleaents.  Finite  eleaent  analysis  is  used  to  determine 
the  strain  energy  distribution  In  the  structure,  and. 
In  this  case,  n  Is  the  number  of  finite  eleaents.  In 
general,  this  equation  la  evaluated  for  each  aode  of 
vibration,  resulting  in  a  aodal  loss  factor  for  each 
aode.  A  previously  developed  three-dimensional  strain 
energy/flnlce  elaaenc  analysis  (Hwang  and  Gibson. 
1991)  based  on  the  fully  expanded,  three-dimensional 
version  of  this  equation  was  used  to  evaluate  the  loss 
factors  for  the  beams  In  the  first  bending  aode  of 
vibration.  The  corresponding  first  mode  natural 
frequencies  were  determined  from  the  eigenvalue  solver 
In  the  finite  eleaent  prograa. 

Experlaents  have  only  been  done  so  far  on  the 
constrained  layer  beeas,  which  consist  of  2026  alloy 
alualnua  beams  of  various  thicknesses  with  3M  SJ-2052X 
constrained  viscoelastic  layer  daaplng  cape  of 
constant  thickness  on  one  side  of  the  beaa.  The 
viscoelastic  adhesive  layer  Is  made  of  3M  ISO- 112 
acrylic  polymer,  end  the  constraining  layer  la  deed 
soft  1100  alualmai.  Unidirectional  T- 300/934 
graphite/epoxy  composite  beams  with  fiber  orientations 
of  0°  and  45®  wero  only  evaluated  analytically,  and 
the  corresponding  experlaents  have  not  been  conducted 
yet. 

A  previously  developed  Impulse -frequency 
response  technique  (Suarez  and  Gibson,  1987)  was  used 
to  evaluate  daaplng  and  frequencies  of  the  beaa 
specimens.  This  technique  is  based  on  the  enalysis  of 
the  frequency  response  spectrua  for  s  specimen  which 
has  been  excited  by  a  saall  electromagnetic  lapulse 
haaaer.  The  cantilever  beaa  version  of  this  apparatus 
is  shown  in  Fig.  1.  Displacement  response  is  monitored 
bv  a  non-contacting  eddy  current  transducer,  and  a 
Fast  Fourier  transform  analvter  Is  used  to  compute  ar.d 
display  the  frequency  response  spectrua.  Once  the 
first  aode  peak  in  the  frequency  response  spectrua  is 
isolated,  daaplng  Is  determined  by  the  half-power 
bandwidth  method,  and  the  natural  frequency  is  equal 
to  the  peak  frequency. 


rtf.  2  show*  ch*  predicted  variation  in  tha 
flrat  flexural  node  loaa  factor  of  tha  beams  aa  a 
function  of  haaa  layer  thleknoaa.  High  daaplng  la 
found  for  thin  conatralaad  layar  baana  duo  to  largo 
ahaer  dafornatlooa  gana rated  In  tha  viscoelastic 
layar.  Aa  tha  baaa  layar  lncraaaaa,  howavar,  tha 
dealing  af  tha  eooatralnad  layar  ban  dacraaaaa  aa  tha 
affact  of  tha  aurfaea  daaplng  treatment  dlalnlahaa. 
While  tha  dating  af  tha  eooatralnad  layar  baaoa  la 
graatar  than  due  of  tha  ccapoilt*  baaoa  for  aoall 
boon  thleknoaa,  tha  c pop o alto  boon  daaplng  lncraaaaa 
and  aurpaaaea  due  of  the  eonatralnod  layer  baaoa  aa 
the  beao  thleknoaa  lncraaaaa.  The  lncroaaa  In  tha 
canpaalto  dapping  with  lncraaa lng  baao  thleknoaa  ta 
believed  to  bo  duo  to  tha  frequency  depandanca  of  tha 
daaplng  in  tha  a  poxy  aatrlx  notarial,  and  to  lncreaaed 
through*  ehlcknaaa  a  hoar  dafornaelona  In  tha  ehlck 
conpoalta  baaoa. 

Additional  raaulea  have  boon  obtained  regarding 
the  contribution  to  total  daaplng  by  tha  varleua 
naterlala  and  tha  different  atreaa  cooponenta  In  tha 
baaoa.  For  axaople,  ahaar  coupling  in  tha  45° 
coapoaite  baaa  and  the  ahaar  atreaa  In  tha 
vlaeoalaatlc  layar  of  tha  conatralnad  layer  baaa  are 
found  to  be  Important  factors  affecting  tha  total 
daaplng.  Experimental  results  for  the  constrained 
layar  baaoa  shoved  good  agreement  with  predicted 
values,  and  the  affect  of  Increasing  baaa  thickness  on 
the  daaplng  was  verified  experimentally,  aa  shown  In 
rig.  3.  In  conclusion,  those  resulta  dearly  show 
that,  while  natal  beams  with  surface  daaplng 
treatments  offer  the  highest  damping  for  chin  beams, 
composite  materials  offer  superior  daaplng  as  tha  baaa 
thleknoaa  increases.  The  Implications  for  real 
structures  having  thick  cooponenta  camot  be  Ignored, 
and  further  research  on  chick  places,  thick  shells  and 
ocher  chick  structures  having  more  complex  geooeertes 
la  warranted. 
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Figure  1.  Impulse-frequency  response  apparatus  for 
measurement  of  frequencies  and  daaplng  of  cantilever 

bean  specimens. 


Figure  2.  Predicted  variation  of  first  node  flexural 
loss  factor  with  base  layer  thickness  for  cantilever 
beam  specimens  constructed  of  alualnua/conatralned 
layer  and  graphlte/epoxy . 


Figure  3.  Variation  of  predicted  and  measured  first 
soda  flexural  loss  factor  with  base  layer  thickness 
for  cantilever  beaa  specimens  constructed  of 
alumlnua/constrained  layar. 
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This  project  is  a  joint  research  program  between  Georgia 
Tech  and  Clark  Atlanta  University  to  investigate  the  elastic 
strain  and  stress  respoose  of  a  filament  wound  gun  barrel  during 
manufacture  cure  and  firing-  Included  in  the  project  is  the  actual 
fabrication  of  various  overwrapped  steel  tubes  with  two  types  of 
fiber  and  resin  at  two  fiber  tensions.  Also  demonstrated  are  the 
advantages  of  joint  projects  between  traditional  engineering 
schools  end  historically  black  institutions  that  are  strong  in 
science.  This  pngect  has  been  fimded  by  die  ARMY  through 
fee  Benet  Labs  in  Wstervliet,  New  York. 

This  paper  will  concentrate  on  the  analytical  solution  of  the 
stress  during  cure  and  firing.  Numerical  results  are  presented. 
The  results  show  that  a  significant  residual  stress  state  will  exist 
because  of  the  processing  technique  Furthermore,  depending 
on  lay-up  and  cure  cycle,  large  localized  stresses  can  be 
generated  which  could  cause  defects  such  as  longitudinal 
wrinkles.  Also  diacussed  are  several  possible  solutions 

Large  bore  gun  can  be  significantly  improved  by 
overwrapping  the  barrel  with  advanced  composite  materials. 
However  winding  pattern,  winding  tension,  cure  cycle,  and  resin 
selection  are  ell  considered  key  factors  affecting  quality  end 
consistency  of  highly  structural  efficient  filament  wound  vessels 
(Lark.  1977).  Selecting  these  is  only  part  of  the  process. 
Concerns  specific  to  a  particular  application  must  be  addressed, 
both  in  terms  of  manufacturing  and  utilization.  This  paper  will 
look  at  the  winding  tension  during  fabrication,  and  then  the 
curing  temperature  and  curing  pressure  during  ramp  up  and 
ramp  down  on  the  residual  stress  after  cure.  Finally,  these 
stresses  are  examined  as  a  source  of  wrinkling  in  the  bands. 
The  basic  problem  is  illustrated  in  Fig.  1  as  a  multilaycd  tube. 


Fig.  1.  Basic  problem  of  multilayered  cylindrical  tube. 

One  problem  with  filament  winding  thick  composites  is  the 
loiyitudmal  wrinkles  that  form  during  the  curing  process.  It  is 
thought  that  this  is  due  to  the  bleeding  of  excess  resin  as  well  as 
fiber  compaction.  The  fibers  accommodate  the  loss  of  volume 
by  folding  over,  farming  wrinkles  (Calabrese,  et  a],  1987).  Even 
an  almost  completely  (by  wind  does  produce  these  wrinkles 


indicating  that  the  problem  is  more  than  just  resin  related. 
Additional  tension  applied  during  fee  winding  could  increase 
compaction  and  help  to  reduce  wrinkling.  Unfortunately,  high 
tension  often  causes  filaments  wound  circumferentially  to  cut 
into  lower  layers,  destroying  them.  Future  work  toward  solving 
this  problem  will  involve  optimizing  fee  pwJcing  of  the  fibers 
during  tiie  winding  process  and  investigating  the  effect]  of  the 
cure  cycle. 

The  stress  distribution  in  e  compound  cylinder  may  be 
found  by  decomposing  the  problem  into  a  number  of  simpler 
ones,  solving  each  of  them,  and  superimposing  the  solutions. 
For  each  layer,  a  spatially  uniform  temperature  gradient  is 
assumed.  By  breaking  the  cylinder  down  into  a  large  number  of 
these  elements,  any  general  temperature  profile  can  be 
approximated.  Equations  have  been  developed  to  allow  the 
calculation  of  stresses  and  displacements  on  orthotropic 
elements.  Additional  boundary  conditions  must  be  enforced  at 
the  interface  of  two  layers  within  the  cylinder,  such  as  matching 
pressure  and  strains  at  the  interface. 

For  fee  winding  tension,  the  method  given  by  Faupel  and 
Fisher  (1981)  is  used  to  calculate  the  tension  stress  on  filament 
wound  cylinders. 


The  analysis  for  an  orthotropic  cylinder  subjected  to 
pressure  loads  closely  parallels  feat  of  an  isotropic  cylinder  and 
is  presented  in  Lddmitritii  (1981).  Witherell  (1990a),  after 
some  simplification  of  Lekhnitskii's  solution,  produced  the 
following  set  of  equations  for  a  monolayered  orthotropic  cylinder 
under  plane  strain  condition  subjected  to  internal  pressure  q  and 
external  pressure  P0. 


For  an  orthotropic  cylinder  subjected  to  a  temperature 
difference,  several  authors  including  Hyer  (1987)  and  Witherell 
(1993)  have  outlined  elasticity  solutions  for  calculating  the 
thermal  stresses  in  multilayered  composite  tubes.  The 
assumptions  are  that  the  stresses  and  strains  are  not  functions  of 
fee  axial  coordinate  and  feat  temperature  is  spatially  uniform 


within  each  layer.  This  implies  that  the  responses  are 
independent  of  theta.  The  equations  for  these  results  are 
somewhat  lengthy  are  not  presented  here. 

For  each  sub-problem,  expressions  exist  for  the  redial  and 
circumferential  stresses.  The  last  set  of  unknowns  that  must  be 
determined  are  the  interlace  pressures.  This  is  accomplished  by 
requiring  that  die  hoop  strain  for  adjacent  layers  be  equal  at 
their  interface.  In  Witherell's  notation,  the  hoop  strain 
equivalence  between  layers  i  and  r+1  is  given  by 

where  the  superscript  T  denotes  thermal  effects  and  'ptf  denotes 
pressure  effects  If  this  condition  is  enforced  for  a  cylinder  with 
n  layers,  there  will  be  n-1  interfaces  and  thus  n-1  equations  to  be 
solved  for  the  unknown  pressures  at  each  interface. 

The  models  described  above  can  be  used  to  determine 
stresses  induced  on  composite  cylinders  due  to  manufacturing  or 
service,  fa  the  following  section,  numerical  results  are  obtained 
and  their  impact  discussed  for  a  cylinder  with  a  steel  liner  end  a 
composite  overwrap  of  IM7  as  an  example.  The  following 
material  properties  (IM7/epoxy)  and  geometry  were  used. 


Composite  Properties 

Mandrel  Properties 

E,  « 146.8  GPa 

Er  -9.9  GPa 

Ex  -9.1  GPa 
viS  “vra“0.3 
vfe-0.49 

a,  -  04  -  33.7  x  lO'V’c 

Oe  « -0.0774  xlO'V’c 

E  -  200  GPa 
v-0.3 

Geometry,  radius 

inside,  a- 34.8mm 
interface,  b- 38.0  mm 
outside,  c  -45.7  mm 

Although  these  properties  represent  a  hoop  lay-up,  die  solution 
technique  is  not  restricted  to  this  type  of  lay-up.  Any  orthotropic 
lay-up  can  be  handled,  even  ±®  wraps  like  those  produced  by  die 
filament  winding  process,  so  long  as  the  angle  pair  is  view  as  a 
single  orthotrepic  layer. 

Winding  tension  far  the  IM7  fibers  (12K)  was  set  at  6  lb  in 
this  study,  which  equates  to  a  113.1  MPa  tensile  stress  During 
the  cure,  the  resin  changes  from  a  liquid  to  a  solid  plastic 
capable  of  transferring  load.  The  stress  due  to  die  winding 
tension,  temperature  increase  in  the  autoclave,  pressure  in  the 
autoclave  and  decrease  of  temperature  in  the  autoclave  are 
modeled  separately  and  are  superimposed  to  give  the  total 
residual  stress.  As  an  example.  Figure  2  shows  die  stress  state 
of  a  cylinder  wound  with  IM7  fibers  at  a  tension  of  6  pounds  and 
cured  using  temperature  only.  Most  authors  have  assumed  that 
the  residual  stresses  can  be  estimated  by  considering  only  the 
temperature  downramp  during  cure.  This  figure  shows  that 
theae  stresses  are  in  fact  the  dominant  ones,  but  that  the  other 
atresaes  are  not  negligible  It  is  interesting  to  note  that  there 
exists  a  significant  tensile  radial  stress  at  die  liner-jacket 
interface,  hi  die  lab,  it  has  been  observed  that  the  jacket  often 
separates  from  die  liner  because  of  die  thermal  mismatch 
between  a  steel  and  composite.  B  may  be  difficult  to  produce  an 
adhesive  bond  strong  enough  to  handle  this  large  load.  The  load 
could  be  reduced  by  increasing  the  winding  tension.  Other 
options  include  autofrettaging.  or  plastically  deforming,  the 
inner  liner  to  force  a  fit  Conversely,  it  may  be  possible  to 
shrink  fit  die  liner  by  taking  it  to  a  very  low  temperature.  B 
would  also  be  worthwhile  to  investigate  die  effect  of  surface 
finish  on  die  interface  properties. 


In  conclusion,  the  equations  that  have  been  developed  are 
useful  to  both  the  design  engineer  as  well  as  the  manufacturer. 
Overall  residual  stress  states  can  be  predicted.  Just  as  important 
is  the  ability  to  predict  the  state  of  stress  during  curing  to 
identify  the  causes  of  manufacturing  problems  such  as 
wrinkling.  Due  to  the  interaction  between  the  winding  tension, 
the  cure  temperature,  material  properties,  and  cure  pressure,  it 
has  been  found  that  significant  compression  can  occur  on  the 
outer  edge  of  the  tubes,  fa  particular,  low  winding  tension,  and 
high  cure  pressure  can  result  in  compress  wrinkles  on  the 
outside  layers.  Future  work  should  involve  accounting  for 
plastic  deformations  because  of  the  potential  to  eliminate  the 
separation  problem  and  to  correlate  the  results  with  the  ongoing 
experimental  phase  of  this  project. 
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To  date,  most  of  the  efforts  that  are  directed  towards  controlling  the  delamination  are  passive 
and  are  focussed  on  preventing  the  formation  of  delaminations.  Several  passive  options  that  have 
been  proposed  include  an  optimization  of  the  lamina  sequence,  an  introduction  of  special  adhesive 
layers,  and  to  serrate  the  fiber  -  matrix  interface.  Even  when  all  these  passive  options  are 
implemented,  delaminations  do  develop  in  practice.  Active  control  of  delaminations  is  an  alternative 
in  such  cases.  In  the  rotorcraft  technology,  some  of  the  applications  include  delamination  control  in 
the  flex  beam  of  bearing  less  rotor  system,  a  composite  hub  and  critical  areas  of  the  airframe.  The 
objective  of  our  studies  are  to  explore  the  feasibility  of  real  time  detection  of  delaminations  in  a 
laminated  composite  structure  and  use  this  information  to  develop  active  controllers  to  prevent  these 
delaminations  from  further  growth.  Steps  that  are  necessary  in  achieving  this  objective  are  as  follows: 
Techniques  are  needed  to  sense  and  detect  the  existence,  size  and  location  of  delaminations.  We  also 
need  dynamic  models  for  the  delaminated  system  that  can  be  used  as  the  plant  model  to  design 
controllers.  Other  items  of  investigation  are  structural  mechanics  foundation  for  controller  design 
and  actuator  options. 

In  this  paper,  we  will  discuss  our  work  in  the  area  of  delamination  detection  techniques  by 
using  the  change  of  dynamic  characteristics  due  to  delaminations.  These  techniques  offer  the 
potential  of  being  used  in  real  time  detection  techniques.  We  have  developed  an  analytical  dynamic 
model  for  delaminated  beams.  Theoretical  solutions  of  the  beam  equations  that  natural  frequencies 
of  delaminated  beam  are  different  from  the  frequencies  of  a  non-delaminated  beam.  The  solutions 
also  indicate  the  existence  of  an  additional  frequency  and  a  corresponding  mode  due  to  the  presence 
of  a  delamination.  This  can  be  called  as  an  opening  mode  of  the  delamination  or  a  delamination 
modes.  We  have  validated  some  of  our  theoretical  studies  by  means  of  experimental  modal  analysis. 
The  presentation  will  include  the  preliminary  work  in  the  area  of  active  control. 
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ABSTRACT 

Bolted  graphite/epoxy  plates  were 
subjected  to  in-plane  biaxial 
loading.  Cruciform  shaped  samples 
were  used.  Testing  was  performed  on  a 
biaxial  siachine  developed  and  built 
at  Concordia  University.  Four  arms  of 
the  cruciform  sample  were 
independently  loaded.  Different 
lateral  loads  were  applied  along  one 
direction.  Along  the  other  direction, 
load  was  applied  on  one  arm  while  the 
opposite  arm  was  held  in  the  grip 
without  applied  load.  This  opposite 
arm  measured  whatever  load  that  was 
bypassed  around  the  constrained  bolt 
hole.  Acoustic  emission  was  used  to 
detect  the  onset  of  failure.  Results 
show  that  there  is  a  significant 
decrease  in  joint  strength  with  the 
application  of  the  lateral  load.  This 
effect  is  more  severe  for  plates  with 
smaller  size  hole  (6.35  mm  diameter) 
as  compared  to  plates  with  larger 
size  hole  (19.05  mm). 

INTRODUCTION 

In  designing  using  composites,  it  is 
advisable  to  avoid  the  joint  as  much 
as  possible.  This  is  because  the 
joint  is  usually  the  weakest  link  in 
a  composite  structure.  However, 
there  are  situations  where  the  joint 
is  not  avoidable.  One  example  is  the 
case  of  patching  repair  of  airframes 
where  patches  need  to  be  joined  to 
the  parent  structure.  Another  example 
is  the  attachment  of  a  composite 
structure  to  other  structural 
elements,  usually  by  rivets. 

The  majority  of  composite  structures 
in  airframes  are  in  the  form  of  thin 
plates  or  shells.  These  structures 
are  subjected  to  biaxial  loading 


conditions  in  real  applications.  It 
is  also  well  known  that  laminated 
composites  are  anisotropic,  this 
means  that  their  response  depends  on 
the  direction  of  the  loads. 

Most  work  on  joints  in  composite 
plates  or  shells  have  considered  only 
uniaxial  loading.  The  joint  is 
allowed  to  float  freely  while  loads 
are  applied  on  the  two  ends  of  the 
assembly  until  the  joint  fails.  This 
type  of  test  may  produce  useful 
information  but  it  may  not  produce 
data  that  are  representative  of  what 
the  structure  is  subjected  to  in  real 
applications . 

In  this  study,  joint  strength  of 
graphite/epoxy  composites  under 
biaxial  loading  condition  is 
considered.  The  configuration  of  the 
joint  has  the  bolted  hole  fixed  with 
the  loads  applied  around  it.  Figure  1 
shows  a  schematic  of  the  load  and 
constraint  arrangement . 


Figure  1:  Biaxial  bearing/bypass  load 
and  constraint  arrangement 


With  the  fastener  held  fixed, 
different  loads  can  be  applied  in 
both  the  x  direction  (right  hand 
side)  and  y  direction  (upper  side) 
while  the  x  direction  (left  hand 
side)  and  y  direction  (lower  side) 
are  allowed  to  remain  free  or  are 
subjected  to  yet  different  loads.  In 
other  words  four  different  loads  can 
be  applied  at  four  loading  arms 
simultaneously.  This  is  possible 
because  the  loads  can  be  supported  by 
the  fixed  fastener. 

Depending  on  the  degree  of  rigidity 
of  the  constraint  of  the  fastener, 
some  of  the  applied  loads  can  be 
bypassed  around  the  fastener.  These 
are  called  bypass  loads.  If  the 
fastener  is  rigid,  there  is  little  or 
zero  bypass  load,  if  the  fastener 
loses  its  rigidity  such  as  the  case 
after  failure  at  the  fastener  hole, 
the  bypass  load  can  be  larger. 

SPECIMEN  DESIGN 


There  appeared  to  be  no  form  of 
specimen  available  for  this 
particular  type  of  test  in  the 
literature.  It  was  therefore 
necessary  to  design  a  specimen  for 
this  particular  test.  The  development 
for  the  design  of  the  cruciform 
specimen  used  in  this  study  was  given 
in  [  1  ]  and  [  2  ] .  This  specimens  has 
the  configuration  shown  in  Figure  2. 


Figure  2:  Configuration  of  biaxial 
specimen. 


RESULTS 

Table  1  shows  the  joint  strength 
results  for  plates  containing  of 
three  different  hole  sizes,  five 
different  lateral  loads  from  0  kN  to 


30  kN  and  a  case  of  no  clamp  load. 
Consider  first  the  lateral  load  as 
the  varying  parameter.  For  plates 
with  small  holes  (6.35  mm  diameter) 
with  a  clamp  load  of  7464  N,  as  the 
lateral  load  is  increased  from  0  kN 
to  15  kN,  the  joint  strength 
decreases  from  about  48  kN  down  to 
about  23  kN,  a  decrease  of  more  than 
50%.  This  is  very  significant  because 
it  has  the  implication  that  design 
data  obtained  from  uniaxial  loading 
is  non-conservative.  For  the  case  of 
plates  containing  holes  of  19.05  mm 
diameter,  as  the  lateral  load 
increases  from  0  kN  to  30  kN,  the 
joint  strength  also  decreases  but  not 
as  significantly  as  for  the  case  of 
plates  containining  a  smaller  hole. 
Consider  now  the  size  of  the  hole  as 
the  varying  parameter.  In  this  case 
it  iB  necessary  to  normalize  the  hole 

size  by  dividing  the  load  by  the 

diameter  of  the  hole.  For  0  kN 
lateral  load,  the  joint  strength  is 
7.56  kN/m  for  plates  with  a  hole 

diameter  of  6.35  mm  and  a  strength  of 
3.25  kN/m  for  plates  with  a  hole 

diameter  of  19.05  mm.  There  is  a 
reduction  of  more  than  50%  of  the 
normalized  strength  as  the  hole  size 
is  increased.  At  a  lateral  load  of  18 
kN,  there  is  also  a  reduction  in  the 
normalized  strength  as  the  hole  size 
is  increased  from  12.7  mm  to  19.05 
mm.  However,  this  reduction  is  not  as 
eevere  as  the  case  of  0  kN  lateral 
load. 

The  case  where  there  is  no  clamp  load 
gives  a  joint  strength  of  only  20  kN 
for  plates  containing  6.35  mm 
diameter  hole.  This  is  the  lowest' 
joint  strength  in  all  of  the  plates 
containing  this  hole  size. 
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Extended  Abstract 


Much  of  the  early  work  on  the 
dynamic  mechanical  response  of  fiber - 
reinforced  composites  used  impact  bend 
tests  of  Izod  or  Charpy  type  [1,2]. 

The  complex  loading  system  in  such 
tests  does  not  allow  a  fundamental 
Interpretation  of  the  data  obtained 
and  prevents  any  simple  analysis  of 
the  effect  of  loading  rate  on  the 
failure  processes.  For  this  reason 
over  recent  years  attempts  have  been 
made  at  the  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory, 
Hanover,  New  Hampshire,  to  develop 
impact  tests  for  fiber-reinforced 
composite  materials  in  which  a  simpler 
loading  system  is  applied  and  account 
can  be  taken  of  stress  wave  effects. 
Version  of  the  the  split  Hopkinson's 
pressure  bar  apparatus  for  simple 
uniaxial  compression  has  been  designed 
and  large  strain  rate  effects  are 
being  investigated  on  48 -ply  GERP, 
where  mechanical  response  is  strongly 
rate  dependent. 

The  CRREL  Hopkinson  Bar  Impact 
Apparatus  consists  of  a  short  striking 
bar  driven  by  compressed  air  and  two 
long,  solid,  round  bars  instrumented 
with  strain  gauges.  The  striker's 
impact  on  the  first  bar  sets  up  a 
compressive  stress  pulse,  whose 
amplitude  depends  on  the  impact 
velocity  and  whose  duration  depends  on 


the  striker's  length.  The  test 
specimen  is  held  between  the  two  bars 
in  a  chamber  cooled  by  liquid 
nitrogen.  At  the  interface,  some 
energy  of  the  stress  pulse  from  the 
first  bar  is  reflected  because  of  the 
mechanical  impedance  mismatch  of  the 
bar  and  the  specimen.  The  rest  is 
transmitted  into  the  specimen, 
propagates  through  it,  and  reaches  the 
second  bar's  interface  with  the 
specimen.  At  the  second  interface, 
part  of  the  stress  pulse  is  again 
reflected  back  into  the  specimen  and 
the  remainder  is  transmitted  to  the 
transmitter  bar.  If  the  specimen  is 
much  shorter  than  the  wavelength  of 
the  loading  pulse,  the  wave- 
transmitting  time  will  be  small 
compared  with  the  duration  of  the 
loading  stress  pulse,  and  after  a  few 
reflections  within  the  specimen,  the 
stress  and  strain  along  the  specimen 
become  approximately  uniform. 

Consideration  of  equilibrium  at 
the  interface  between  the  specimen  and 
the  transmitter  bar  shows  that  the 
forces  in  the  specimen  and  in  the 
transmitter  bar  are  equal.  The  force 
is  determined  by  measuring  the  strain 
on  the  elastic  transmitter  bar  using 
strain  gauges  and  the  average  strain 
in  the  specimen  is  calculated  from  the 
displacements  at  the  end  of  the 
specimen  using  equations  (1) ,  (2) ,  and 
(3). 


specimen 

strain 

Es(t)-L2C 

(1) 

s  0 

spec imen 

strain  rate 

de  (O  - 

=  e  =  e  (t),  and 

ut  s  L  v 

s 

(2) 

specimen 

s  t  ress 

o  ( t )  =  E  e  ( t ) 
s  b  t 

(3) 

where 

c  =  wave  propagation  velocity  in  the  Hopkinson  bars 
Ls  ■  length  of  specimen 
er(t)  =  instantaneous  reflected  strain 
et(t)  *  instantaneous  transmitted  strain 
Eb  *  elastic  modulus  of  the  bar  material. 


Figure  1  illustrates  the  test 
system. 

Currently  an  investigation  is  in 
progress,  supported  by  U.S.  Army  CRREL 
to  study  the  influence  of  low 
temperature  on  the  compression  impact 
response  of  a  48 -ply  Fiberite  974 
epoxy  reinforced  with  T-300  graphite 
fiber,  in  an  attempt  to  determine  how 
far  the  Increased  strength  and  energy 
absorbing  potential  of  different 
orientations  and  layers  of 
reinforcements  at  all  rates  can  be 
optimized.  Figure  2  shows  some 
initial  data  from  high -strain -rate 
impact  on  two  graphite/epoxy  laminates 
impacted  transversely  at  lower  than 
fracture  load  in  the  Hopkinson  bar. 

In  the  A-l  laminates  all  plies  are 
oriented  in  the  same  direction 
< CO] aa  1  >  -  B-l  laminate  also 

consisted  of  48  plies  but  arranged  in 
six  layers  as  [0#/90,/0, Jg .  Note  that 
the  stress -strain  curves,  shown  in 
Figure  2(e)  and  2(f)  for  the  two  types 
of  laminate  are  significantly 
different.  The  average  elastic 
modulus,  E  .  of  the  [0]4,  laminate  is 
0.53  x  10  *  psi  (3.66  GPa)  and  that  of 
the  [04/90,/0,]  is  0.95  x  10«  psi 
(6.55  GPa).  As  a  result,  although 
both  laminates  were  subjected  to 
nearly  the  same  magnitude  of  impact 
stress,  17,500  psi  (120.7  MPa),  the  A- 
1  laminate,  I0]4i,  developed  a  peak 


strain  of  e  -  30  x  10"’  at  the  strain 
rate  c  -  700  strains  s'1;  whereas  for 
B-l  laminate,  [0,/90S/0*J  ,  e  -  19  x 
10~*  at  e  -  550  strain  s  *.  For  any 
particular  laminate  type  the  strain 
rate,  c,  can  be  systematically 
Increased  by  increasing  the  impact 
velocity  or  geometry  of  the  specimen. 
The  above  results  clearly  indicate 
that  the  impact  response  of  the 
layered  composites  is  dependent  on  the 
ply  geometry  and  ply  properties. 

Since  the  properties  do  change  under 
environmental  conditions,  for  example 
high  and  low  temperatures  and 
humidities,  we  propose  to  study  these 
impact  responses  under  the  varying 
environmental  conditions  for  which 
facilities  have  been  developed  at 
CRREL. 
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Abstract 


In  a  series  of  impact  experiments  on  composite 
plates  which  were  aimed  at  testing  the  conjecture 
that  the  perforation  threshold  represented  a  critical 
phenomenon  and  that  Percolation  Theory  could 
eventually  provide  an  adequate  description  of  the 
spallation  or  fragmentation  process  accompanying 
the  perforation  of  anisotropic  layered  materials,  the 
fragments  emerging  from  a  0/90  layup  fiber  epoxy 
plate  target  were  collected  and  photographed.  Image 
analysis  techniques  were  used  to  deduce  the  size 
distribution  of  the  backside  fragments.  It  was  found 
that  the  number  N  of  fragments  having  area  S  may 
be  represented  as  a  power  law  N  -  1  /  S  J  .  The 
average  exponent  for  the  experiments  which  were 
performed  at  nine  different  values  of  impact  velocity 
over  the  range  from  300  to  6000  fps  had  the  value  J 
■  2.17  i  18  %  ,  which  lies  between  hie  previously 
predicted  theoretical  values  from  Percolation  Theory 
of  J(2D)  *  2.05  and  J(3D)  *  Z2.  For  the  velocity  V, 
*  397  fps  which  is  dose  to  the  perforation  threshold 
velocity,  J  had  the  value  2.03  ,  which  is  surprisingly 
dose  to  the  two  dimensional  Critical  Percolation 
Limit. 

Introduction 

The  mechanical  properties  of  a  composite  are 
severely  degraded  by  impad  induced  damage.  A 
significant  fraction  of  the  projectile  kinetic  energy  is 
channelled  into  rival  modes  of  damage,/  i.e. 
delamination,  fiber  fracture,  matrix  cracking  and 
fragmentation.  Some  disintegration  appears  around 
and  inside  the  penetration  hole,  and  at  the  backside 
of  the  target  spall  is  produced  even  at  velocities  less 
than  the  ballistic  limit.  More  energy  is  consumed 
above  this  limit,  as  intense  fragmentation  emerges 
from  the  target.  This  process  may  be  studied 
quantitatively  by  measuring  the  fragment  size 
distribution. 

Under  a  short  impulse  of  pressure,  a  situation  typical 
of  either  an  impact  or  an  explosion  is  created.  The 
stress  acts  in  a  local  way  and  is  'unable”  to  scan  the 
medium  for  its  weakest  points.  Therefore,  many 
cracks  grow  simultaneously  and  independently  in 
different  regions  of  the  target.  These  conditions  lead 


to  the  well  known  result  (4)  that  the  average  fragment 
size  becomes  smaller  the  higher  the  strain  rate. 
Moreover,  they  provide  the  physical  justification  for 
using  percolation  theory  in  predicting  fragment  size 
distributions.  (2,3) 

Percolation  theory  (1)  has  been  recognized  as  a 
statistical  physics  model  for  critical  phenomena.  By 
applying  it  to  fracture  and  fragmentation  problems  we 
were  able  to  show  that  fragmentation  (as  well  as  the 
coalescence  of  many  microcracks  forming  a 
macroscopic  fracture)  is  a  critical  phenomenon. 
Fragmentation  occurs  whenever  the  density  of  cracks 
reaches  its  critical  value.  That  is  to  say,  that  when 
the  probability,  p,  of  creating  a  microcrack  achieves 
the  critical  value,  pc.  For  crack  densities  higher  than 
pe,  the  critical  density  of  microcracks,  the  number  of 
small  particles  in  the  population  grow  at  the  expense 
of  the  reduction  of  th f.  number  of  laige  ones. 
Consequently,  the  total  surface  area  of  the  fragments 
is  growing.  It  is  therf ore,  expected,  that  the  minimal 
surface  energy  required  for  fragmentation  is  received 
at  pe  (  below  pc  a  complete  shattering  of  the  finite 
macroscopic  piece  of  material  is  not  possible). 

Coming  back  to  impact  physics,  the  determination  of 
the  limit  velocity  below  which  perforation  of  a  given 
target  by  a  given  projectile  is  not  possible  continues 
to  be  of  much  interesL  Even  the  definition  of  the 
ballistic  limit  (5)  (the  critical  value  of  striking  velocity 
at  which  the  probability  of  perforation  is  50%) 
suggests  that  stochastic  mechansms  are  at  play. 
Below  the  ballistic  limit,  it  is  expected  that  the  spalled 
fragmented  zone  at  the  backside  of  the  sample  and 
the  shattered  and  penetrated  zone  in  the  front  side 
will  be  well  separated  by  a  zone  of  delaminated 
material  strong  enough  to  prevent  the  projectile  from 
perforating  the  target.  In  the  limit  as  this  intact  zone 
becomes  very  thin,  these  two  zones  will  touch  each 
other,  and  the  sample  will  then  become  disintegrated 
throughout  its  thickness.  It  is  of  interest  to  inquire 
whether  the  critical  condition  for  fragmentation 
supplied  by  percolation  theory  is  able  to  yield 
additional  information  about  the  ballistic  limit  and  the 
physical  conditions  of  its  occurrence. 

Description  of  the  Experiment 


Image  Analysis 

The  population  of  fragments  from  each  experiment 
was  collected  on  the  inner  surface  of  a  cylindrical 
sheet  of  paper  coated  with  a  sticky  adhesive.  After 
the  experiment,  the  cylinder  was  rolled  flat  and 
photographed.  Such  a  procedure  has  to  be  carried 
out  with  caution  in  order  to  avoid  any  optical  noise 
resulting  from  shadows,  unequal  illumination  or 
spurious  particles. 

The  photographs  wens  scanned  and  a  digitized  image 
was  obtained.  The  next  step  is  the  division  of  the 
digital  image  into  dusters  of  material  (fragments)  and 
the  background.  This  procedure  is  not  easy  when 
the  gray  levels  of  the  fragments  and  local 
background  overlap.  Also,  it  is  not  easy  to 
distinguish  two  fragments  partially  overlapping  each 
other.  Special  data  treatments  were  adopted  to 
reduce  distortion  and  deterioration  of  the  information 
recorded  in  the  digital  picture,  which  was  scanned  by 
a  S12  X  S12  black  and  white  scanner.  Higher 
resolution  of  the  original  photograph,  of  course,  will 
result  in  higher  final  resolution  and  consequently 
more  data  on  possible  lost  information  with  regard  to 
the  very  fine  fragments. 

The  measured  fragment  size  histograms  vary 
continuously  in  an  irregular  way.  In  order  to 
overcome  the  effects  of  random  fluctuations,  the 
results  were  accumulated  in  bins  having  unequal 
width.  Each  bin  further  along  on  the  size  axis  was 
made  twice  as  large  as  the  preceding  one.  The  total 
number  N  of  occurences  in  a  bin  was  divided  by  the 
bin  width,  to  obtain  the  average  distribution  density  of 
sizes.  A  linear  least  squares  fit  of  log  N  against  log 
S  was  performed  in  order  to  deduce  the  parameters 
of  the  expected  power  law  relationship.The  results 
are  presented  hi  Figure  1.  The  slopes  and  intercepts 
of  the  linear  fits  are  presented  in  Table  1  for  each 
velocity.  , 

The  role  of  strain  rate  in  controlling  fragmentation 
was  recognized  by  Grady  &  in  the  early  eighties. 
The  highest  strain-rate  found  during  an  impact  may 
be  estimated  by  tire  ratio  of  the  projectile  to  its 
radius,  de/dt  ~  Vp  /  Rp  strain  rates  at  the  backface 
may  be  estimated  by  multiplying  this  value  by  some 
attenuation  function  for  the  particle  velocity,  in  any 
case,  it  is  expected  that  the  higher  the  impact  energy 
is,  the  smaller  will  be  the  average  fragment  size.  It  is 
not  surprising  therefore  that  the  trend  of  intercepts  in 
Table  1  really  represents  the  fraction  of  fine  particles 
in  the  distribution.  The  behavior  of  the  slopes  needs 
further  consideration.  The  mass  of  a  fractal  structure 
having  a  length  L  (l»1)  is  given  by  M(L)  -  A  L  D 
where  D  is  the  fractal  dimension.  At  the  critical  point, 
p«pe  the  value  of  D  is  known.  For  homogeneous 


percolation  in  the  plane,  D  *  1.9  (2D).  In  three 
dimensional  space  6  =  2.5  (3D).  An  additional  result 
of  percolation  theory  is  that  the  number  n,  of 
clusters  having  s  sites  (per  site)  is  given  by  n,  - 1  /  S 

J  ,  where  J  =  (d+D)  /  D.  and  where  the  exponent  J 
has  the  value  J  =  2.05  in  two  dimensions  and  J*  22 
in  three  dimensions  (d=3).  All  the  slopes  in  Table  1 
are  around  these  values  with  an  error  not  exceeding 
about  20%.  Furthermore,  in  the  introduction,  a 
possible  relation  between  the  perforation  threshold 
and  critical  fragmentation  was  conjectured.  The 
lowest  velocity  in  Table  1  is  not  far  from  the 
perforation  threshold  found  experimentally  namely 
V0  *  375  fps.  The  slope  of  the  corresponding 
fragment  size  distribution  is  extremely  dose  to  the 
theoretical  value  of  J  for  the  two  dimensional  case.  If 
the  fragments  of  the  impacted  composite  plate  have 
a  fixed  thickness  (  as  would  be  the  case  if  it  they 
resulted  from  the  fragmentation  of  previously 
delaminated  plies),  then  the  measured  area 
distributions  would  represent,  in  fact,  two  dimensional 
mass  distributions.  In  a  laminated  material  with 
planar  weaknesses  between  the  plies,  such  a 
hypothesis  appears  plausible.  In  general  isotropic 
solids,  a  true  three  dimensional  behaviour  is 
expeded  with  J  =  22.  It  is  interesting  to  check 
whether  fragment  size  distributions  from  impacted 
brittle  soldis  also  behave  in  the  same  way  near  the 
perforation  limit.  More  work  is  needed  in  order  to 
correlate  other  statistical  Physics  features  either  with 
dynamical  predictions  or  with  experimental  data. 
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Low  velocity  impact  induced  damage  as¬ 
sessment  requires  knowledge  of  contact  behavior, 
transient  elastic  wave  propagation,  crack  nucle- 
ation,  and  propagation  and  interaction.  The  initial 
crack  formation  takes  place  in  and  near  the  contact 
region  and  therefore  is  sensitive  to  the  indentation 
behavior.  A  detailed  dynamic  stress  field  based  on 
the  realistic  contact  laws  [1]  provides  information 
about  highly  stressed  zones  in  the  undamaged  tar¬ 
get  composite  structure.  These  zones  can  be  iden¬ 
tified  as  crack  initiation  domains.  The  stress  field 
in  these  crack  initiation  domains  is  used  to  identify 
the  stress  components  responsible  for  cracks  nu- 
cleation.  Once  the  crack  nucleation  is  recognized, 
propagation  can  be  studied  with  the  help  of  the 
fracture  mechanics  principles.  However,  it  should 
be  recognized  that  the  dynamic  stress  field  at  every 
stage  of  propagation  is  very  complex  and  therefore 
intractable  to  detailed  study.  It  is  clear  that  the 
elastic  stress  field,  although  complex,  is  well  de¬ 
fined  up  to  the  point  of  initial  fracture.  Joshi  and 
Sun  [2,3]  presented  a  simple  study  incorporating 
the  above  mentioned  notions  of  the  impact-in¬ 
duced  damage  initiation  in  laminated  composites. 

The  impact  damage  in  laminated  compos¬ 
ites  has  been  an  active  research  area  for  a  long 
time  which  is  evident  from  the  vast  amount  of  lit¬ 
erature  available  in  this  area.  The  purpose  of  this 
paper  is  to  discuss  some  important  features  of  con¬ 
tact  behavior,  transient  incident  energy  propaga¬ 
tion  and  crack  nucleation  with  the  help  of 
numerical  and  experimental  results. 

Determination  of  spatial  and  temporal  dis¬ 
tribution  of  contact  traction  is  important  in  impact 
studies.  It  is  apparent  that  the  traction  distribution 
over  the  contact  area  can  have  a  profound  influ¬ 
ence  on  the  nature  of  the  near  contact  stress  field. 
For  example,  compare  the  Boussinesq  stress  filed 
[4]  in  a  half  space  subjected  to  point  load  to  a  Hu¬ 
ber  stress  field  [5]  with  distributed  load  according 
to  Hertz  analysis.  The  spatial  distribution  of  the 
impact  load  should  not  be  ignored  when  studying 
impact-induced  fracture  initiation  occurring  in  the 
near-contact  region. 


The  temporal  characteristics  of  the  impact 
loading  are  also  important.  The  difference  in  rates 
at  which  the  incident  energy  is  imparted  to  the 
structure  and  propagated  away  from  the  impact 
point  is  responsible  for  the  high  stress  build-up  in 
the  near-contact  region.  The  contact  area  and  the 
traction  distribution  changes  with  change  in  re¬ 
sultant  contact  load  with  time.  In  a  dynamic  con¬ 
tact,  spatial  and  temporal  characteristics  of  the 
contact  are  coupled.  The  use  of  static  indentation 
laws  may  be  considered  and  is  justified  under  the 
following  assumptions: 

1 .  The  total  contact  duration  is  a  order  of  magni¬ 
tude  longer  than  the  time  taken  by  the  longitudinal 
waves  to  travel  through  the  thickness  of  the  target 
composite,  assuming  thickness  is  the  smallest  di¬ 
mension. 

2.  The  rate  of  change  of  contact  area  is  a  order  of 
magnitude  smaller  than  the  characteristic  inplane 
longitudinal  wave  velocity. 

3.  The  impactor  dimensions  are  comparable  to  the 
thickness  dimension  of  the  target  or  it  is  very  stiff 
compared  to  the  target.  This  rules  out  effect  of 
elastic  wave  propagation  in  the  target. 

The  numerical  results  in  this  study  are  ob¬ 
tained  by  incorporating  static  indentation  laws  in 
conjunction  with  the  plate  and  impactor  dynamics. 
The  temporal  distribution  of  contact  force  is  ob¬ 
tained  by  employing  a  plate  finite  element  analysis 
[5]  with  the  static  experimentally  obtained  inden¬ 
tation  laws  [6].  The  effect  of  spatial  variation  of 
the  impact  load  on  through-the-thickness  distribu¬ 
tion  of  stresses  is  captured  by  considering  a  plane 
strain  problem  of  an  infinitely  wide  plate  impacted 
by  a  rigid  cylindrical  impactor.  The  contact  length 
during  the  impact  is  calculated  by  modified  Hert¬ 
zian  indentation  law,  providing  elliptical  force  dis¬ 
tribution  over  the  contact  length.  The  previously 
obtained  (from  the  plate  analysis)  temporal  distri¬ 
bution  is  implemented  as  resultant  impact  load  in 
the  cylindrical  impactor  problem.  The  quantitative 
information  about  the  loading  is  lost  in  doing  so. 
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while  the  qualitative  nature  of  the  impact  loading 
is  retained.  The  contact  length,  obtained  using  the 
Hertzian  law,  is  scaled  to  match  the  experimental¬ 
ly  obtained  contact.  The  spatial  distribution  of  the 
load  is  assumed  to  be  elliptical  and  is  applied  as 
consistent  nodal  forces  within  the  contact  length. 
The  nodal  forces  are  modified  at  every  time  step 
integration,  depending  on  the  contact  length  and 
resultant  impact  load. 

Elastic  waves  travel  at  different  phase 
speeds  in  different  directions  because  of  anisotro¬ 
pic  material  properties  of  the  target.  The  initial 
stress  field  in  the  vicinity  of  impact  is  solely  due  to 
incident  energy  if  the  geometric  boundaries  are 
sufficiently  away  from  the  impact  point.  The  rela¬ 
tive  importance  of  reflected  energy  (geometric 
boundary  conditions)  in  initiating  damage  may  be 
resolved  by  comparing  the  time  taken  by  the  dom¬ 
inating  frequency  component  of  the  initial  pulse  to 
reach  the  impact  region  after  reflection  within  the 
impact  duration  itself.  It  should  be  noted  that  in¬ 
terference  between  initial  pulse  and  the  reflected 
pulse  can  be  constructive  or  destructive.  We  will 
restrict  the  discussion  to  non-interference  cases. 

We  expect  three  dimensional  stress  field  in 
the  vicinity  of  the  impact  center  (a  region  with  di¬ 
mensions  of  the  order  of  maximum  contact  length) 
and  two  dimensional  (plate  or  shell  type)  stress 
field  in  the  far-region.  The  thickness  dimension 
plays  an  important  role.  The  Near-Impact-Stress- 
Field  (NISF)  is  considerably  different  in  thick  and 
thin  laminated  composites.  The  NISF  is  also  al¬ 
tered  by  the  boundary  conditions  on  distal  surface 
(free  surface  or  surface  resting  on  rigid  support 
etc.)  which  is  the  nearest  surface  from  impact 
point.  The  transverse  normal  stress  in  NISF  has  to 
satisfy  compressive  normal  traction  on  the  contact 
surface  and  zero  traction  on  the  rest  of  the  proxi¬ 
mal  surface  and  the  distal  surface;  this  surface  is 
free  in  most  of  the  practical  problems.  In  general, 
the  impact  duration  is  far  greater  than  the  time  tak¬ 
en  by  a  longitudinal  wave  to  travel  through  the 
thickness  and,  therefore,  the  through  the  thickness 
distribution  of  the  transverse  normal  stress  will  not 
be  altered  by  the  reflection  from  the  distal  surface. 
This  means  that  the  problem  is  quasi-static  in  the 
thickness  direction  which  also  justifies  use  of  stat¬ 
ic  indentation  laws.  These  two  attributes  of  the 
NISF  (boundary  conditions  and  quasi-static  na¬ 


ture)  restricts  transverse  normal  stress  from  be¬ 
coming  significantly  tensile. 

Crack  nucleate  in  region  of  high  stress 
concentration.  These  stress  produced  cracks  orient 
themselves  perpendicular  to  the  maximum  tensile 
principal  stress  direction  in  the  brittle  isotropic 
material  by  a  smooth  impactor  is  experimentally 
verified  [7].  In  case  of  anisotropic  solids,  it  is  nec¬ 
essary  to  take  into  account  the  orientation  depen¬ 
dence  of  both  elastic  constants  and  surface  energy. 
Clearly,  in  composites  with  strong  directional 
cleavage  tendencies,  the  anisotropy  in  surface  en¬ 
ergy  will  govern.  In  general,  though,  we  expect 
some  compromise  between  the  tendencies  for 
cracks  to  be  orthogonal  to  maximum  tensile  stress 
and  weaker  cleavage  planes.  Each  ply  of  the  lami¬ 
nated  composite  plate  may  be  considered  as  trans¬ 
versely  isotropic,  with  plane  of  isotropy  being 
orthogonal  to  fiber  orientation  of  the  ply.  There¬ 
fore,  the  crack  will  nucleate  orthogonal  to  the  ma¬ 
jor  component  of  tension  in  the  isotropic  plane  of 
the  ply.  Initial  impact  damage  can  be  explained 
with  the  help  of  this  hypothesis. 

Representative  numerical  and  experimen¬ 
tal  results  will  be  included  in  the  full  length  paper 
to  substantiate  the  arguments  and  elaborate  on  the 
various  aspects  of  the  low  velocity  impact  present¬ 
ed  above. 

This  research  is  partially  funded  by  NATO 
grant  No.  CRG  920  125. 
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An  analytical  study  is  conducted  to  pre¬ 
dict  the  nonlinear  transient  response  of  thin  imper¬ 
fect  laminated  curved  panels  subjected  to  impact 
loads.  Although  a  large  number  of  studies  deal¬ 
ing  with  impact  response  of  flat  laminated  plates 
has  been  conducted  in  the  past,  relatively  little 
work  seems  to  be  available  on  impact  of  laminated 
curved  panels.  The  imperfect  panels  are  modeled 
using  a  doubly-curved  thin  plate/shell  element  capa¬ 
ble  of  modeling  arbitrary  geometric  imperfections. 
The  element  was  previously  presented  by  Kapa- 
nia  aad  Yang  (1986)  and  was  later  extended  to 
study  the  large-displacement  behavior  of  laminated 
anisotropic  plates  and  shells  (Saigal  etnl.,  1986)  and 
the  buckling,  postbuckling  aad  nnnKw«r  vibrations 
of  imperfect  plates  (Kapania  and  Yang,  1987). 

The  element  is  quadrilateral  in  shape  and  has 
four  nodes,  one  at  each  comer.  Each  node  possesses 
12  DOFs:  u1,  dux/d£,  dv.lJdi),  and  sim¬ 

ilar  quantities  for  u3  and  «*,  where  u1,  u3  and  u3 
are  the  displacement  components  in  Cartesian  coor¬ 
dinates  xl,  x3  and  x3,  respectively,  and  (  and  rj  are 
the  curvilinear  orthogonal  coordinates  on  the  mid¬ 
dle  surface  of  the  shell. 

T.iiua,  nit  nonlmear  transient  responses  are 
obtained  using  the  modal  superpoaition  method 
(MSM)  and  a  reduction  method  baaed  on  Bits  vec¬ 
tors  (R-R),  both  in  conjunction  with  direct  integra¬ 
tion  schemes  (Wflaon-0  at  Newmark  methods)  (Ka¬ 
pania  and  Byun,  1993).  The  effect  of  the  number 
of  modee  used  is  also  discussed.  A  modified  contact 
law  is  incorporated  to  evaluate  the  impact  loads  due 
to  a  projectile.  Different  non-dimensional  shell  radii 
[r/h]  are  used  in  order  to  study  the  effect  of  the  non- 
dimensional  shell  radii  [r/h]  on  the  impact  response 
of  laminated  curved  panels.  The  effect  of  geometric 
imperfections  on  linear  end  nonlinear  transient  re¬ 
sponse  under  sudden  impact  loads  is  also  analysed. 
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The  laminated  composite  curved  plate  consid¬ 
ered  is  a  20- layered  [0*/45#/0*/  —  45*/0*]a»  lami¬ 
nate.  The  plate  dimensions  a  and  0  are  15.24  cm 
and  10.16  cm,  respectively.  The  thickness  of  the 
plate,  denoted  by  ho,  is  0.269  cm.  The  curved 
plate  is  assumed  to  be  impacted  at  the  center 
by  a  spherical  steel  impactor  of  diameter  1.27  cm 
with  an  initial  velocity  («o)  of  3  m/s  aa  shown  in 
Figure  1.  The  mass  density  of  the  impactor  is 
7.96  x  10~,Ns2cm~‘i.  The  boundary  conditions 
are  simply  supported  and  immovable  in  the  planes 
along  the  plate’s  edges.  The  material  properties  of 
a  graphite/epaxy  lamina  are  assumed  to  be:  Ej  s 
120  GPa,  Ei  =  7.9  GPa,  ui3  =  0.30,  Gxa  *  * 

Gis  =  5.5  GPa,  and  p  =  1.58  x  10~‘ 


Fig.  1.  A  laminated  curved  plate  subjected  to  im¬ 
pact  load  at  the 

The  finite  element  mesh  used  in  this  analysis  is 
an  8x8  finite  element  meah  (for  the  whole  plate). 
The  contact  force  histories  and  central  deflections 
of  the  curved  plate  were  obtained  by  using  both  the 
modal  superposition  (MSM)  aad  the  Bits  reduction 
(R-R)  methods  for  three  different  shell  radii.  The 
nondimensbnalaed  shell  radius  defined  aa  r0/ho 
waa  used.  The  contact  force  histories  and  central 
deflections  of  the  curved  plate  were  obtained  using 
a  time  step  of  2 pa. 

From  the  comparison  between  the  two  methods 
it  was  found  that  both  show  good  agreement  in  the 
contact  force  aad  displacement  histories  of  a  lam¬ 
inated  curved  plate  subjected  to  impact  loads.  In 


the  present  iterative  procedure  for  impact  response 
analysis,  the  time  increment  needed  to  be  very  small 
in  order  to  get  accurate  contact  forces.  For  lin¬ 
ear  analysis,  the  Rayleigh- Rits  reduction  method 
based  on  Rits  vectors  required  less  CPU  time  than 
the  mode  superposition  method  although  the  results 
were  very  similar. 

Both  msthods  required  a  large  number  of  ba¬ 
sis  vectors  to  account  for  the  impact  load  which  has 
high-frequency  characteristics.  Furthermore,  for  the 
nonlinear  analysis,  the  mode  superposition  method 
required  less  CPU  time  as  shown  in  Thble  1.  It 
should  be  noted  that  since  most  of  the  time  is  spent 
on  the  calculation  of  the  internal  nodal  force  vector 
the  use  of  reduction  methods  in  performing  non¬ 
linear  impact  response  analysis  may  not  result  in 
significant  reduction  in  CPU  time.  The  need  for 
development  of  an  efficient  way  of  handling  the  cal¬ 
culation  of  the  internal  load  vector  is  apparent.  Also 
note  that  baas  updating  should  be  performed  prop¬ 
erly  in  order  to  obtain  accurate  results  for  nonlinear 
impact  response  studies. 

Table  1  Comparison  between  the  CPU  time  re¬ 
quired  for  the  modal  superposition  method  (MSM) 
and  the  Bits  reduction  method  (R-R)  (nonlinear 
analysis). 


With  the  reduction  of  the  oon-dimensiooal  shell 
radius  [r/h]  the  curved  plate  appeared  to  be  stifier. 
The  aoo-dimwnaional  shell  radius  [r/A]  for  values 
less  than  1000  caused  <-h»itg—  jg  the  im¬ 

pact  response  of  the  laminated  curved  plate.  Fig¬ 
ure  2  shows  a  comparison  between  the  central  de¬ 
flection  histories  for  three  different  non-dimensional 
shell  radii  [r/h].  The  introduction  of  geometric  im¬ 


perfections  led  to  similar  conclusions,  that  is  a  stiffer 
structure  as  shown  in  Figure  3. 

M 
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Fig.  2.  Comparison  between  the  central  deflection 
histories  for  three  different  non-dimensional  shell 
radii  (r/h)  (plate  thickness:  0.269cm). 


Fig.  3.  Effect  of  simple  geometric  imperfections  on 
central  deflection  for  non-dimensional  shell  radius 
(r/h)ss  1000  (plate  thickness:  0.1345cm). 

REFERENCES 

Kapaaia,  R.  K.,  and  Byun,  C.t  “Reduction 
Methods  Bawd  on  Eigenvectors  and  Rits  Vectors 
for  Nonlinear  Transient  Analysis,*  Cwyrisfossl 
Mechanic*,  An  International  Jonmol,  VoL  11,  No. 
1, 1993,  pp.  65-82. 

Kapaaia,  R.  K.,  and  Yang,  T.  Y.,  “Formulation 
of  an  Imperfect  Quadrilateral  Doubly-Curved  Shell 
Element  for  Post-Budding  Analysis.”  AIAA  /.,  VoL 
24,  No.  2, 1986,  pp.  310-311. 

Kapaaia,  R.  K.,  and  Yang,  T.  Y.,  “Buckling, 
PnsthnckHng,  and  Nonlinear  Vibrations  of  Imper¬ 
fect  Plates,"  AIAA  /.  ,  VoL  25,  No.  10,  1987,  pp. 
1338-1346. 

Saiga!,  S.,  Kapaaia,  R.  K.,  and  Yang,  T.  Y.t 
“Geometrically  Nonlinear  Finite  Element  Analysis 
of  Imperfect  Laminated  Shells,”  Jomn sef  of  Compu¬ 
te  Materials,  VoL  20,  March  1986,  pp.  197-214. 


Control  Designs  for  on: art  Structures  and 
Their  Experimental  Validation 


Farshad  Khorrami 
email:  khorrami@pueel.poly.edu 
Control/Robotics  Research  Laboratory 
School  of  Electrical  Engineering  &  Computer  Science 
Polytechnic  University,  Six  Metroiech  Center 
Brooklyn,  NY  11201 


I.  Introduction 


with  a  piezoelectric  type  accelerometer  at  its  tip  and 
several  pieces  of  piezoceramics  along  the  links. 


Due  to  earth-based  and  space-based  applications, 
and  with  the  promised  advent  of  light-weight  high- 
strength  composite  materials,  much  attention  nas  been 
given  to  modeling  and  control  of  large  flexible  struc¬ 
tures  and  weapon  pointing  systems.  Another  class  of 
systems,  for  which  vibration  suppression  is  of  great 
importance,  is  high  speed  positioning  devices  such  as 
probes  for  electronic  circuit  boards,  robotic  manipula¬ 
tors,  and  machine  tools.  Utilization  of  smart  or  intelli¬ 
gent  materials  in  control  of  systems  exhibiting  flexibil¬ 
ity  is  receiving  increased  attention  recently.  Research 
and  development  of  completely  integrated  structures 
with  embedded  actuators,  sensors,  signal  processing, 
and  control  systems  are  of  prime  interest.  Different 
electroactive  materials  such  as  piezoceramics,  shape 
memory  alloys,  electrorheological  fluids,  polymer  bio- 
materials,  magnetorestrictive  materials  are  being  con¬ 
sidered  for  active  control  purposes. 

Several  control  design  strategies  for  flexible  mutli- 
body  systems,  large  flexible  structures,  and  smart 
structures  are  briefly  outlined  in  this  paper.  To  com¬ 
plement  our  theoretical  studies  of  these  systems,  sev¬ 
eral  experimental  setups  have  been  and  are  being  de¬ 
veloped  at  CRRL.  The  emphasis  of  the  experimental 
research  in  CRRL  is  on  control  of  flexible  mechanical 
structures,  smart  structures,  and  robotic  systems.  In 
the  next  section,  two  of  the  experimental  setups  devel¬ 
oped  at  the  CRRL  are  described.  Control  objectives 
for  these  experiments  range  from  nonlinear  controllers 
to  fuzzy  logic  and  neural  network  control  systems  for 
active  vibration  suppression  to  slewing  ana  pointing , 
and  combinations  thereof.  Furthermore,  system  ana 
parameter  identification  for  control  is  another  primary 
objective  in  the  development  and  utilization  of  these 
experiments. 

IL  Experimental  Setups 
A.  A  Two- Link  Flexible  Manipulator 

A  two-link  robot  arm  with  replaceable  links  has  been 
developed  at  the  CRRL.  Different  configuration  (i.e., 
rigidity)  can  be  studied  by  varying  the  lengths  and  the 
thickness  of  the  arms.  Figure  1  shows  a  picture  of  the 
experimental  setup.  The  actuator  for  the  first  link  is  a 
direct  drive  DC  motor  and  the  second  link  is  actuated 
by  a  geared  DC  motor  through  anti-backlash  gears. 
The  second  joint  has  been  designed  to  be  light-weight 
however,  the  joint  and  the  tip  are  supported  on  an  air 
bearing  on  a  granite  table.  Each  arm  is  instrumented 


Figure  1:  The  two  link  flexible  arm  at  CRRL. 


B.  A  Slewing  Truss  Structure  with  Active  Members 

A  slewing  truss  structure  has  been  developed  (Fig¬ 
ure  2).  The  struts  and  the  nodes  for  the  structure  are 
the  same  as  the  ones  used  in  NASA’s  and  JPL’s  se¬ 
tups.  One  of  the  nodes  is  attached  to  a  direct  drive 
motor  through  a  hub  assembly.  An  active  member  is 
being  installed  in  the  truss  setup  by  replacing  one  of 
the  struts  by  a  piezo  type  linear  actuator.  This  will 
provide  an  active  member  for  vibration  suppression 
along  the  truss.  We  are  also  instrumenting  the  truss 
by  various  sensors. 


Figure  2:  The  slewing  truss  structure  with  active  mem¬ 
bers. 


Figure  3:  Open  and  closed-loop  responses. 


III.  Control  Design  Methodologies 


A.  Adaptive  Nonlinear  Control 

Flexible  articulated  structures  are  nonlinear  dis¬ 
tributed  parameter  systems.  It  is  well  known  that  the 
vibrational  frequencies  of  such  systems  are  geomet¬ 
ric  configuration  dependent.  Therefore,  linear  based 
control  designs  tend  to  be  conservative  to  accommo¬ 
date  for  such  variations.  We  have  developed  nonlinear 
based  control  designs  to  significantly  reduce  the  fre¬ 
quency  variations  due  to  geometric  configuration.  This 
scheme  is  based  upon  asymptotic  expansions  on  the 
nonlinear  distributed  parameter  dynamics.  Further¬ 
more,  nonlinear  adaptive  controllers  have  been  con¬ 
sidered  to  guarantee  robust  performance  in  presence 
of  parametric  perturbations. 

B.  Decentralized  Control 

By  definition  of  a  smart  structure,  sensors,  actu¬ 
ators,  signal  processing,  and  control  algorithms  are 
embedded  in  a  smart  structure.  We  will  refer  to  a 
collection  of  sensors  and  actuators  and  the  associated 
signal  processing  and  control  algorithm  as  a  "smart 
patch.”  In  envisioned  smart  structures,  many  smart 
patches  may  be  utilised  to  perform  different  tasks  and 
therefore  these  systems  will  exhibit  a  decentralized  in¬ 
formation  structure.  Therefore,  it  is  highly  desirable 
to  have  decentralized  control  strategies  meaning  that 
each  smart  patch  in  the  system  will  have  its  own  de¬ 
cision  making  rules  and  algorithms.  This  is  highly  de¬ 
sirable  for  the  following  reasons:  1)  No  need  for  in¬ 
formation  transfer  among  "smart  patches,”  2)  Fault 
tolerant,  i.e.,  failure  of  one  "smart  patch”  does  not 
cause  the  failure  of  the  overall  system,  and  3)  Ease  of 
real-time  computational  load  on  each  subsystem. 

We  are  currently  pursuing  decentralized  control 
strategies  for  the  systems  under  consideration.  A  typ¬ 
ical  response  of  the  advocated  controllers  to  suppress 
vibrations  of  a  slewing  beam  is  given  in  Figure  3.  The 
beam  has  surface-mounted  piezoceramics  as  sensors 
and  actuators. 


Figure  4:  Architecture  of  the  fuzzy  adaptive  controller. 


C.  Fuzzy  Logic  Adaptive  Control 

Fuzzy  controllers  are  based  on  three  well  known 
stages:  1)  the  fuzzification  stage,  2)  the  rule  base  stage, 
and  3)  the  defuzzification  stage.  First,  the  partition¬ 
ing  of  the  fuzzy  input  space  is  carried  out  and  the 
corresponding  membership  functions  of  the  fuzzy  sets 
are  assigned.  Thereafter,  the  system  parameters  in¬ 
volved  must  be  identified  and  their  membership  val¬ 
ues  in  the  fuzzy  input  sets  evaluated  (fuzzification). 
We  have  developed  a  technique  utilizing  a  fuzzy  rule 
base  to  switch  among  several  controllers  designed  for 
different  operating  regimes  and/or  different  parameter 
sets.  The  technique  utilizes  a  real-time  system  identi¬ 
fication  algorithm  during  the  fuzzification  step.  Any 
parameter  estimation  method  may  be  employed  to  ob¬ 
tain  the  unknown  parameters  or  variables.  The  overall 
structure  of  the  scheme  is  depicted  in  Figure  4. 

D.  Neural  Network  Control 

Artificial  Neural  Networks  (ANNs)  have  shown 
great  potential  for  adaptive  control  applications  due  to 
their  capability  to  learn  the  unmodeled  system  char¬ 
acteristics  through  nonlinear  mappings.  This  feature 
is  desirable  in  the  smart  structures  due  to  the  inherent 
nonlinearities  of  the  system.  We  are  currently  utiliz¬ 
ing  a  multi-layer  feedforward  neural  network  with  a 
Trajectory  Pattern  Method  (TPM)  for  identification 
and  trajectory  tracking  of  nonlinear  systems  (espe¬ 
cially  non-minimum  phase  systems).  The  TPM  uti¬ 
lizes  specific  time  functions  for  generation  of  feedfor¬ 
ward  signals.  The  neural  network  is  advocated  to  re¬ 
duce  the  trajectory  error  by  optimizing  the  appropriate 
parameters  such  as  base  harmonics  in  the  TPM.  The 
training  algorithms  used  for  the  networks  are  the  back 
propagation  technique  and  the  genetic  algorithms. 
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Extended  Abstract 


Since  "Desert  Storm"  in  1990,  interest  in 
military  high  technology  has  came  to  an  all  time 
high,  especially  in  contemporary  weapon  systems 
that  use  electronics  and  composite  materials. 

The  potential  for  significant  weight  re¬ 
ductions  in  weapon  system  through  the  use  of 
advanced  composites  was  first  reported  on  a  broad 
scale  in  the  Air  Force  project  "Forecast"  con¬ 
ducted  in  1963.  This  observation  was  based  on 
the  then  recent  development  of  high-modulus, 
high-strength,  low  density  fiber  and  the  superior 
mechanical  properties  that  could  be  developed 
when  these  fibers  were  converted  into  composites 
laminates. 

In  military  aircraft,  composites  offer  im¬ 
proved  performances.  A  number  of  aircraft  de¬ 
sign  and  construction  programs  are  in  progress 
involving  major  use  of  advanced  composites.  For 
example,  the  European  Fighter  Aircraft  (EFA)  is 
designed  using  composite  wings  and  forward  fuse¬ 
lage.  The  use  of  composite  materials  will  result  in 
an  approximately  35%  weight  reduction  as  com¬ 
pared  to  that  of  conventional  materials.  Another 
example  is  the  construction  of  the  AV-8B  MK.  5 
Harrier  II  airframe,  where  the  use  of  composite 
materials  give  rise  to  a  25%  weight  saving  in  the 
airframe  compared  to  an  all-metal  construction. 

Recent  efforts  to  increase  ballistic  pro¬ 
tection  and  reduce  the  weight  necessary  for  this 


protection  have  led  to  a  new  focus  in  composite 
materials  among  the  armor  community.  A  vital 
aspect  of  armor  design  is  an  understanding  of  the 
properties  which  control  anti-ballistic  perform¬ 
ance.  Unfortunately,  for  these  new  materials,  the 
fundamental  research  necessary  for  this  under¬ 
standing  has  only  recently  received  attention. 
This  is  mainly  a  consequence  of  the  fact  that  the 
processes  occurring  during  the  projectile  pene¬ 
tration  are  at  extremely  high  strain  rates. 

As  armor-vehicle  designers  are  forced  to 
reduce  weight  and  increase  protection,  fibers  are 
being  combined  with  ceramics,  metals  and  plas¬ 
tics  to  create  more  effective  armors.  For  example 
the  MIA1  Abrams  MBT  Chobham  armor,  and  the 
M2/M3  Bradely  LAV  (Light  Armor  Vehicle)  in¬ 
terior  spaced  laminate  armor,  are  constructed  of 
these  composite  materials.  The  principle  advan¬ 
tage  of  fibers  is  their  high-strength-to-weight  ra¬ 
tio.  Other  properties  such  as  tensile  modulus, 
elongation,  energy  absorption  and  sound  speed 
have  proven  to  be  important  in  characterizing  im¬ 
pact.  Impact,  however,  is  a  complex  problem  and 
one  should  not  rely  on  any  single  property  as  a 
predictor,  often  combinations  of  properties  are  the 
key. 

Other  products  made  of  composite  ma¬ 
terials  in  the  military  include  the  KE-EP  light¬ 
weight  helmet  and  the  KE-EP  protective  vest. 


However,  very  likely  the  most  important  applica¬ 
tion  can  be  in  the  next  generation  'Rail  Gun” 
weapon  system,  of  which  the  "rail”  part  that  plays 
a  critical  role  could  be  made  of  composite. 

The  increasing  use  of  composites  in  the 
military  field ,  as  well  as  in  the  aerospace,  marine, 
and  automotive  fields,  has  spurred  great  interest 
in  developing  numerical  tools  to  optimize  designs 
made  of  composite  materials  Laminated  compos¬ 
ites,  by  their  nature,  are  heterogeneous  and  ani¬ 
sotropic.  This  makes  the  failure  characteristics  of 
the  composite  materials  more  complicated  than 
that  of  the  isotropic  materials.  Only  in  the  most 
recent  years,  the  finite  element  technique  together 
with  the  high  capacity  and  high  speed  of  modern 
computers  have  made  it  possible  to  analyze  fiber 
reinforced  composites  in  detail. 

Some  well-known  commercial  finite  ele¬ 
ment  codes  are  capable,  to  a  limited  degree,  of 
analyzing  composite  structures.  However,  few 
programs  attempt  to  perform  a  complete  thermo¬ 
mechanical  failure  analysis  of  composite  struc¬ 
tures.  Recently,  under  a  grant  from  the  NASA 
Goddard  Space  Flight  Center,  a  general  purpose 
thermo-mechanical  finite  element  program,  called 
COMPASS  (COMPosite  Analysis  of  Structural 
Systems),  was  developed  for  the  analysis  of  com¬ 
posite  structures.  Three  main  solution  sequences 
are  offered  by  the  code,  namely:  1)  non-linear- 
static  stress  analysis,  2)  delamination  growth 
analysis,  3)  progressive  failure  analysis. 

Focus  of  this  work  is  cm  the  progressive 


failure  analysis  which  is  the  third  case  available 
within  COMPASS.  All  elements  assume  linear 
stress-strain-temperature  relations.  Failure  of  a 
composite  element  is  determined  by  using 
Hashin's  failure  criterion,  while  the  failure  of  an 
isotropic  elements  is  determined  by  using  the  von 
Mises  failure  criterion.  The  solution  process 
involves  increasing  the  original  load  until  an 
element  fails.  The  element  that  failed  is  classified 
as  being  isotropic  or  composite.  Furthermore, 
composite  element  failure  is  identified  as  being  in 
fiber  or  matrix  mode.  The  stiffness,  the  heat 
conductivity,  and  the  thermal  expansion  matrices 
of  that  element  are  then  altered  according  to  the 
mode  of  failure.  As  a  consequence,  the  stresses 
and  the  temperature  in  the  entire  structure  are 
redistributed.  The  structure  is  then  analyzed  with 
the  set  of  matrices  to  determine  if  any  further  fai¬ 
lure  is  induced.  If  further  failure  does  not  occur, 
the  load  is  increased  to  the  level  needed  to  cause 
the  next  element  to  fail.  This  incremental  process 
continues  until  catastrophic  failure  occurs.  The 
ultimate  load  carrying  capacity  of  the  structure  is 
then  determined. 

The  constitutive  relations  of  anisotropic 
thermo-visco-elastic  solid,  the  failure  criteria  and 
post-failure  material  behavior,  and  the  finite  ele¬ 
ment  formulation  of  composite  structures  made  of 
thermo-mechanical  solids  are  presented.  Also,  for 
illustrative  purpose,  the  procedures  and  the  nu¬ 
merical  results  of  the  progressive  failure  analysis 
of  a  composite  plate  are  discussed  in  detail. 
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Summary 

An  extensive  experimental  program  to 
investigate  the  damage  and  energy 
exchange  which  accompany  the 
ballistic  impact  and  penetration  by 
rigid  spherical  steel  projectiles  of 
carbon  fiber  epoxy  composite  plates 
above  their  penetration  threshold 
has  been  conducted  .  Re w 
significant  findings  concerning 
these  impact  phenomena  are 
described.  Preliminary  applied 
mechanics  analyses  that  serve  for 
their  explication  or  correlation  are 
derived. 


Description  of  the  Experimental 
Program 

The  experimental  program 
consisted  of  firing  one-half  inch 
diameter  rigid  spherical  steel 
projectiles  at  carbon  fiber  epoxy 
matrix  composite  plates  which  were 
clamped  between  steel  plates 
containning  a  3  inch  diameter 
circular  hole  with  the  design  impact 
point  located  at  its  center.  The 
composite  target  plates  were  32 
plies  thick  layed  up  in  balanced 
quasi-isotropic  stacking 
arrangements  consisting  of  0,  45, 
and  90  degree  fiber  orientations  and 
also  in  balanced  spiral  staircase 
layups  in  which  the  fibre  direction 
increased  in  increments  of  90,  45, 
22.5,  or  11.25  degrees  from  ply  to 
ply.  Ten  inch  diameter  rigid 
plastic  cylinders  with  axes 
coincident  with  the  shotline  were 
lined  with  adhesive  paper  and  placed 
flush  with  the  front  and  rear  target 
support  plates  in  order  to  catch 
spall  fragments  eaunating  from  the 
target  plate.  Circular  end  plates 
containing  a  small  central  hole  to 
permit  uninterrupted  passage  of  the 
projectile  were  also  covered  with 
adhesive  paper  to  catch  fragments 


travelling  with  velocity  vectors 
contained  within  a  small  cone  of  the 
axis.  Projectile  velocity  was 
measured  ahead  of  the  target  plate 
by  means  of  a  set  of  fine  trip  wires 
of  known  separation  distance  by 
noting  the  instants  of  break  of 
electrical  continuity  and  behind  the 
target  by  two  sets  of  induction 
coils  of  known  location  along  the 
pipe  axis  which  were  wrapped  around 
the  outer  cylindrical  surface  of  the 
plastic  pipe.  Thus  the  energy 
imparted  by  the  projectile  to  the 
target  material  could  be  computed. 
Somm  experiments  were  performed 
without  the  adhesive  paper  and  the 
loose  fragments  collected  in  the 
plastic  pipe  were  swept  together  and 
weighed  and  compared  with  the  weight 
loam  of  the  target  plates.  A  mingle 
set  of  plates  consisting  of  plates 
which  had  been  penetrated  at 
different  velocities  were  stained 
with  Gold  Chloride,  pyrolysed  in  an 
air  oven,  deplied  with  a  scalpel, 
mounted  on  a  rigid-  backing  and 
subjected  to  measurement  of  the 
areas  of  the  penetration  hole  and  of 
the  envelope  of  surrounding 
gold-stained  delaaunated  area.  This 
was  done  with  the  aid  of  automatic 
image  analysis  equipment.  It  was 
thus  possible  to  obtain  the  shapes 
and  distribution  of  hole  and 
delaminated  areas  on  each  ply  of 
panels  penetrated  at  different 
velocities  over  the  range  of  200  to 
6000  feet  per  second.  An  additional 
novel  measurement  was  made  of  the 
total  surface  and  internal  fracture 
area  which  communicates  with  the 
surface  by  employing  a  surface  area 
measurement  apparatus  based  on  the 
gas  adsorbtion  technique  with 
Nitrogen  as  the  adsorbent  gas. 


Principal  Findings 
1 .  Energy  exchange  between 


projectile  and  target. 

The  .fraction  of  the  original  kinetic 
energy  associated  with  the 
projectile  after  conpietion  of  the 
Impact  event  at  firet  increaaea  with 
increase  in  the  Initial  impact 
velocity  and  then  appears  to 
approach  a  constant  aaymptotic 
value.  Logically,  thia  aaymptotic 
value  must  vary  with  the  thickneaa 
of  the  target  plate,  being  unity  for 
a  plate  of  vanishing  thickneaa  and 
zero  for  one  of  semi-infinite 
thickneaa . 

2.  Type,  morphology  and  diatribution 
of  the  daange. 

The  obaerved  damage  conaiated  of  a 
generally  irregularly  ahaped 
perforation  hole,  fragmented 
material  which  originated  from  the 
region  of  the  hole,  and  a 
delaminated  region  surrounding  the 
hole  located  between  each  pair  of 
pliea  in  which  there  was  a  change  in 
fiber  orientation. 

3.  Extent  and  variation  of  the 
total  delasunated  surface  with 
impact  velocity. 

The  curve  depicting  the  variation  of 
the  total  delaminated  area  with 
incident  velocity  of  the  projectile 
possesses  a  bathtub-like  shape.  The 
area  enclosed  by  the  outer  envelope 
of  the  delaminated  area  at  first 
decreases  with  increasing  projectile 
velocity  from  its  value  at  the 
penetration  threshold,  and  then, 
after  paaaing  through  a  minimum 
value,  increaaea  steadily  with 
further  increase  in  velocity. 


Relevance,  Utility  and  Value  of  the 
Research  Findings  to  the  interests 
of  the  Armed  Services 

The  initial  objectives  of  the 
reported  research  are  tied  to  the 
increased  presence  of  manmade  or 
composite  materials  of  the  general 
kind  investigated  in  thia  study  in 
the  construction  of  military 
vehicles  and  weapon  platforms  as  a 
result  of  now  well  recognised 
structural  efficiency  advantages 
such  as  strength  to  weight  ratio  and 
potential  tailorability  of  strength 
in  required  directions  to 
accommodate  the  directionality  of 
imposed  loads.  Decisions  to  employ 
new  materials  in  the  construction  of 
military  materiel  entail  a 


responsibility  to  understand  the 
associated  repairability  issue  that 
stems  from  their  susceptibility  to 
realistic  types  and  magnitudes  of 
damage  from  impact  by  ballistic 
projectiles,  and  the  effect  on 
vulnerability  as  a  result  of  their 
influence  on  the  performance  of  the 
engaged  projectile  and  of  the  degree 
of  strength  retention  in  the 
presence  of  realistic  types  and 
magnitudes  of  damage.  The 
investigation  thus  far  has  uncovered 
three  much  sought  after  items  of 
information  concerning  the 
interaction  of  spherical  projectiles 
with  composite  plates.  In 
particular,  the  first  finding  that 
the  degree  of  damage  is  bounded  from 
above  by  the  dastage  incurred  at  the 
perforation  threshold  is  considered 
a  useful  piece  of  knowledge  for 
maintenance  personnel.  A  second 
relevant  finding  is  the 
understanding  gained  from  the 
present  investigation  that  enables 
one  to  infer  the  form  of  the 
momentum  interchange  with 
projectilies  in  general  and  which 
consequently  allows  correlation  of 
the  thresholds  for  various  types  of 
functioning.  Thirdly,  the  finding 
concerning  the  asymptotic  value 
which  is  rapidly  approached  by  the 
fraction  of  initial  kinetic  energy 
retained  by  the  projectile  after 
penetrating  a  composite  panel  is  an 
important  piece  of  information  to 
aid  in  assessment  of  potential  for 
additional  damage  by  a  spherical 
projectile  which  has  penetrated  a 
multiplicity  of  panels  to  structure 
located  farther  along  the 
trajectory ,  or  in  the  design  of  a 
single  plate  capable  of  yielding  a 
specified  degree  of  attenuation  of 
the  .incident  energy  of  an  incident 
steel  sphere. 
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THE  DYNAMIC  RESPONSE  OF  THICK  ANISOTROPIC  COMPOSITE  PLATES 
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V*  develop  Che  analysis  and  nuaerical 
calculations  for  Che  response  of  thick 
anisoCropic  elastic  plate  strips  subjected  to 
transient  loadings.  Problems  associated  with 
the  vibration  of  mechanical  systems  have  long 
been  of  interest  in  structural  applications. 
Most  of  the  available  work  deals  with  various 
aspects  of  plates  constructed  of  only 
isotropic  materials.  Even  here,  mostly 
approximate  results  are  available  which  are 
based  on  classical  bending  theory  of  plates. 
Typically,  the  place  is  assumed  to  be  "Chin" 
so  that  variations  of  displacements  and 
stresses  through  the  plate* a  thickness  are 
neglected  and  that  the  plate  is  of  Infinite 
extent,  as  presented  by  Chow  (1971)  for  an 
orthotropic  laminated  plate. 

Comparatively  speaking,  however,  .only  a 
limited  number  of  results  have  appeared  on 
free  vibration  of  thick  bounded  plates, 
especially  thick  anisotropic  plates.  Much 
less  work  is  available  on  the  response  of  such 
structures  to  transient  loading  conditions. 

Specifically,  our  resulting  structural 
component,  made -up  of  multloriented 
unidirectional  fiber  composite  laminates,  is 
an  anisotropic  plate  having  monoclinic 
symmetry  Special  combinations  of 


orientations  of  such  laminates  may  result  in 
higher  than  monoclinic  symmetry.  However,  no 
combinations  of  such  laminates  would  lead  to 
less  chan  monoclinic  symmetry.  Thus,  the 
problem  we  deal  with  encompasses  all  classes 
of  symmetry  greater  than  trlclinic  symmetry. 

The  analysis  of  the  plate  scrip  of  Fig.  1 
is  conducted  in  two  steps. 


In  the  first  step,  the  characteristic  free 
vibrational  modes  of  the  strip  are  derived.. 
Linear  transformations  are  utilized  to  derive 
the  sixth  degree  characteristic  equation  and 
the  secular  equation  for  this  case  in  closed 
form.  Ve  subsequently  isolate  the 
mathematical  conditions  for  symmetric  and 
antisymmetric  wave  mode  propagation  in 
completely  separate  terms,  as  developed  by 
Nayfeh  and  Chlmenti  (1988) 

D11C  1coth(7ll)  -  Dl3C3coth(7*3) 

♦  DjjGjCoth(7ij)  -  0,  (la) 


011C1c«nh(TA1)  -  Dl3C3t«nh(7A3) 

+  °15C5  t«nh(7A5)  -  0.  (lb) 

W#  us«  the  two  secular  aquations  (la,b), 

together  with  the  charateristic  equation,  to 
obtain  the  proper  combinations  of  A  and  7  for 

valid  solutions.  In  the  second  step,  the 
applied  loads  are  expanded  in  teras  of  these 
normal  modal  functions,  and  the  response  of 
the  plate  is  obtained  by  linear  superposition 
of  the  appropriate  components.  Hater ial 
systems  of  higher  symmetry,  such  as 
orthotropic,  transversely  isotropic,  cubic, 
and  isotropic  are  contained  implicitly  in  our 
analysis. 

We  write  the  formal  solutions  (Jones 


(1971))  for  the  displacements 
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(2a, b) 


In  this  figure,  we  present  a  plot  of  natural 
frequency  versus  the  nondlmenslonal  quantity 
nd/L  (discrete  wave  number)  for  a  plate  with  4 
-  45*.  Lines  labeled  A^  and  correspond  to 

the  1th  antisymmetric  and  symmetric  natural 
modes,  respectively.  A  surprising  feature  of 
this  plot  is  the  crossing  of  the  second  and 
third  antisymmetric  modes  by  the  third 
symmetric  mode. 

Finally,  we  develop  the  transient  response 
of  the  plate  when  subjected  to  an  impulsive 
line  loading  applied  along  the  x2-axls. 


A* 
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Introduction 


Notched  Strength  Evaluation 


Fiber-reinforced  composites  have  been  used  in 
many  applications  of  Army  service  to  replace  exotic 
alloys  due  to  their  lightweight,  corrosion-resistance, 
and  high  strength.  While  the  strength  of  notched 
structures  under  various  loadings  has  always  been  of 
concern  to  the  Army,  satisfactory  evaluation 
approach  has  not  been  available  due  to  the  limitation 
of  existing  theories.  The  discovery  of  hole  sixe  effect 
made  it  dear  that  the  daseical  design  based  on  the 
stress  concentration  factors  may  no  longer  be 
acceptable  for  notched  composites.  Mainly  for  this 
reason,  great  efforts  have  been  made  since  the  early 
1970’s  to  evaluate  the  notched  strength  of  laminated 
composites  with  cutouts.  As  a  result,  several  models 
have  been  developed  eo  for  including  the  linear 
elastic  fracture  mechanics  model  (LEFM),  finite 
stress  distribution  models,  damage  zone  model,  and 
progressive  failure  model.  The  finite  stress 
distribution  model  drew  more  attention  than  other 
models.  Whitney  and  Nuismer  (1974)  proposed 
"point  stress  criterion  (PSC )"  and  "average  stress 
criterion  CAJSC)"  in  the  study  of  orthotropic  laminates 
containing  a  circular  opening  or  a  crack  under  tensile 
loads  by  assuming  that  the  strength  of  a  notched 
plate  was  controlled  by  the  normal  stress  adjacent  to 
the  circular  hole.  These  models  have  been  extended 
by  many  researchers  since  then. 

Almost  all  of  the  previous  studies  in  this  area  were 
limited  to  the  in-plane  loadings.  Although  composite 
laminates  with  cutouts  under  bending  have  many 
applications,  very  few  studies  have  been  reported. 
The  main  reason  is  the  difficulties  in  stress  analysis. 
In  this  study,  a  simplified  approach  baa  been  used 
which  combines  the  experiments  and  finite  element 
analysis  to  evaluate  the  notched  strength  of 
laminates  under  bending.  This  notched  strength 
depends  on  many  factors  such  as  the  geometry  of  the 
opening,  width  and  thickness  of  the  plate,  material 
properties,  fiber  angles,  and  stacking  sequence.  This 
paper  is  the  first  part  of  a  series  of  investigation  and 
focuses  on  the  laminates  with  elliptical  holes  under 
bending.  Further  studies  will  indude  laminates  with 
various  holes  under  differ  'at  loadings  and  the  affect 
of  stacking  sequence. 


For  an  isotropic  laminate  under  in-plane  tension, 
the  normal  stress  is  uniformly  distributed  across  the 
laminate  thickness.  But  when  the  laminate  is 
subjected  to  a  bending  load,  the  stress  varies  in  the 
thickness  direction.  This  difference  also  exists  in  the 
laminated  plates.  For  a  cross-ply  laminated  plate 
under  the  in-plane  tensile  load,  ah  of  the  0*  plies  are 
in  the  same  stress  level,  and  ah  the  90*  plies  are  also 
in  the  same  stress  level.  If  the  stress  in  one  0s  ply 
(or  90*  ply)  reaches  the  ultimate  strength,  all  the  0* 
plies  (or  90*  plies)  will  fail  simultaneously.  But  for 
the  cross -ply  laminates  under  bending,  each  0*  ply  or 
90*  ply  ia  in  the  different  stress  level.  The  stress  at 
one  point  or  the  average  stress  over  a  length  may  not 
represent  the  stress  level  or  etnas  concentration  of 
the  laminate.  Therefore,  these  criteria  may  not  be 
directly  applied  to  the  case  of  banding.  It  is  mom 
reasonable  to  indude  the  stress  variation  across  the 
laminate  thickness  in  the  failure  models.  By  this 
consideration,  the  following  failure  models  wan 
presented  to  evaluate  the  notched  strength  of 
laminates  under  bending. 


Modified  Point  Stress  Failure  Model  -  This  model 
extends  Whitney  and  Nuismer’s  point  stress  model 
by  defining  the  "point  stress"  as  the  average  stress  of 
the  plies  with  same  fiber  angle  across  the  laminate 
thickness.  For  the  laminates  other  than  on-axis 
unidirectional  or  cross-ply,  the  maximum  stresses 
usually  will  not  occur  along  a  radial  line  parallel  to 
the  applied  bending  moment  but  at  some  angles. 
Thus  the  modified  point  stress  criterion  should  be 
applied  along  a  path  which  endoses  the  elliptical 
opening,  aa  shown  in  Fig.  1. 
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when  o.  ia  the  stress  of  the  ply  with  fiber  angle  a, 
h,  is  the  total  thickness  of  the  pHes  with  fiber  angle 
a,  h  is  the  total  thickness  of  the  laminate,  a,  is  the 
ultimate  strength  of  the  material  along  or 
perpendicular  to  the  fiber  direction,  and  st  is  a 
characteristic  parameter  to  evaluate  the  failure 
strength  of  the  laminate. 


Modified  Average  Stews*  Failure  Modal  -  In  the 
ease  of  in-plane  loading*,  the  average  ■  trees  criterion 
includes  the  contribution  of  the  stress  distribution 
along  the  radial  direction  in  the  vicinity  of  the  hole 
edge.  For  notched  laminates  under  bending,  the 
contribution  of  stress  distribution  along  the  thickness 
direction  should  also  be  considered.  Mathematically, 
this  modified  model  can  be  written  as 

— r-  /  /  o«C*.  r,  ijfirdt.  -  o.  (2) 

Jt  0 

where  R  is  the  radial  distance  from  the  edge  to  the 
center  of  the  hole,  6  is  the  location  of  the  section 
where  maximum  stress  oocurs,  and  a,  is  another 
characteristic  parameter  to  predict  the  failure 
strength  of  the  laminates  under  bending. 


Fig.  1  A  Laminate  with  sue  Elliptical  Hole 
under  Bending 

As  discussed  earlier  in  the  introduction,  a 
theoretical/analytical  solution  to  the  problem  of 
notched  laminate  under  bending  would  be  extremely 
complicated  and  very  difficult  to  achieve.  Finite 
element  method  aeems  to  be  the  best  alternative  to 
serve  the  purpose.  This  has  been  realized  by  many 
investigators  when  studying  tensile  failure  of 
notched  laminates.  In  this  current  investigation,  a 
simplified  approach  was  utilized  which  combines  the 
sample  experiments  and  finite  element  stress 
analysis.  A  brief  description  of  this  procedure  is 
given  as  follows: 

•  Conduct  4-point  bending  test  and  record  the  Toad 
vs.  deflection"  curve. 

•  Obtain  'first  ply  failure  load"  from  Toad  vs. 
deflection"  curve  for  each  sample. 

•  Calculate  the  strength  reduction  factors  for  plates 
with  different  holes. 

•  Apply  modified  PSC  or  ASC  to  the  plate  without 
hole  to  determine  the  failure  load  of  the  plate. 

•  Determine  failure  loads  for  plates  with  different 
holes  using  strength  reduction  factors. 

•  Conduct  finite  element  stress  analysis  by 
applying  the  failure  loads  to  the  notched 
composite  laminates. 

•  Apply  modified  PSC  and  ASC  to  each  case  to 
determine  the  parameters  ax  and  a* 


Experimental  Study 

A  four-point  bending  test  was  performed  to 
determine  the  failure  load  of  laminated  composites 
containing  an  elliptical  hole  under  bending.  A  dial 
gauge  was  attached  to  the  bottom  surface  of  the 
sample  to  measure  the  deflection.  Scotchply 
laminated  samples  containing  an  elliptical  hole  were 
tested.  Each  sample  was  a  9"x  8"x  1/8"  [((VBOyO], 
laminated  plate.  Five  different  elliptical  hole  sizes 
were  used.  To  determine  the  failure  loads  for  each 
sample,  Toad  (force)  vs.  defection"  curve  was  plotted 
during  the  experiment  and  was  later  be  used  to 
calculate  the  strength  reduction  factors. 

Stress  Analysis 

The  finite  element  method  was  utilized  in  the 
stress  analysis  of  composite  laminates  with  different 
holes  under  bending.  The  software  COSMOS/M, 
which  is  capable  of  handling  laminated  composites, 
was  used.  A  finite  element  code  was  also  developed 
by  the  authors  to  compare  the  numerical  results. 
The  results  from  COSMOS/M  and  the  developed  code 
are  consistent.  The  stress  distributions  for  each  case 
were  obtained. 

Results  and  Discussions 

With  stress  distribution  provided  by  finite  element 
analysis  and  strength  reduction  factors  from 
experiments,  the  characteristic  parameters  ax  and  a, 
were  obtained  by  applying  modified  PSC  and 
modified  ASC.  The  parameters  ax  and  a,  evaluate 
the  notched  strength  of  laminate  under  bending  in 
two  ways.  The  parameter  ax  focuses  on  the  intensity 
of  stress  at  one  point  in  the  vicinity  of  the  opening, 
while  a,  considers  file  stress  distribution  in  the 
radial  direction  near  the  hole.  It  was  found  that  for 
an  elliptical  opening,  the  larger  the  aspect  ratio,  the 
smaller  the  strength  reduction  factors.  It  was  also 
found  that  both  a:  and  a,  are  reduced  as  the  aspect 
ratio  increases. 
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ABSTRACT 

The  salient  features  of  smart  structures  are  order  structural  systems  is  examined  that  uses 


the  generation  of  internal  forces  for  suppressing 
vibrations,  mode  shape  control  and  the  spatially 
distributed  nature  of  the  actuators  and  sensors. 
In  the  conventional  control  design  methodology, 
the  spatial  distribution  of  sensors  and  actuators 
is  not  taken  into  consideration.  Recently,  a 
number  of  researchers  have  proposed  designing 
distributed  control  schemes  for  infinite-order 
structural  systems  using  distributed  sensors  and 
actuators.  The  main  advantages  of  distributed 
control  are  (i)  the  ability  to  control  all  modes  of 
the  structure  thus  eliminating  spill-over 
problems,  (ii)  trade-offs  between  gain  margin 
and  bandwidth  of  the  controllers,  and  (iii)  the 
use  of  displacement,  rotation  and  strain 
measurements  as  feedback  signals  in  the 
infinite-order  systems.  The  design  and 
implementation  of  optimal  control  strategies  for 
infinite-order  smart  structures  present 
challenging  problems.  One  is  the  computational 
effort  required  to  design  optimal  controllers  for 
infinite-order  structural  models,  and  the  other  is 
the  development  of  optimally  shaped  distributed 
sensors  and  actuators.  A  procedure  for  sensor 
shape  optimization  and  the  design  of  controllers 
for  simple  cantilever  beam  test  article  has  been 
developed  and  discussed  in  this  paper. 

The  smart  materials  polyvinylidene 
fluoride  (PVDF)  and  shape  memory  alloys  can 
be  utilized  in  the  development  of  customized 
distributed  sensors.  These  sensors  can  be  shaped 
to  control  the  desired  response  of  the  structural 
system.  For  the  implementation  of  optimal 
control  strategies,  we  need  to  develop  sensor 
shape  optimization  techniques.  In  this  paper,  a 
method  for  the  design  of  controllers  for  infinite- 


finite  dimensional  approximations  to  desired 
curvature  and  curvature  rate  kernels.  These 
kernels  are  used  for  generating  the  shapes  of 
PVDF  distributed  sensors.  The  output  of  these 
spatially  distributed  sensors  will  be  utilized 
directly  as  the  control  signal  for  suppressing 
the  structural  vibrations. 

A  method  for  the  calculation  of 
distributed  sensor  shape  functions  was 
presented  by  Miller  and  van  Schoor  in  1991. 
A  finite  element  model  (FEM)  of  the  system 
was  used  to  obtain  a  state  space  model  of  the 
structural  system.  With  the  use  of  the  FEM 
interpolation  function,  shape  functions 
representing  a  continuous  approximation  to 
the  discrete  linear  quadratic  regulator  (LQR) 
solution  were  determined.  The  shape 
functions  were  then  transformed  into  the 
appropriate  shape  functions  required  for  the 
implementation  with  PVDF  distributed 
sensors.  This  method  is  computationally 
intensive  and  requires  covering  the  entire 
structure  with  two  distributed  sensors  for  the 
implementation  of  the  optimal  conUol 
solution.  Also,  only  the  LQR  control  design 
solution  was  demonstrated  and  the  shape 
functions  were  not  verified  with  simulations. 

In  order  to  verify  the  design  method 
discussed  above,  simulations  using  the 
distributed  sensors  were  necessary.  Initial 
simulations  of  the  LQR  controller  using  the 
state  space  model  of  the  structural  system 
were  performed.  The  results  of  this  simulation 
were  considered  to  be  the  standard  by  which 
the  performance  of  the  system  with  the 


distributed  sensors  was  evaluated.  A  second 
simulation  of  the  system  was  performed  to 
evaluate  the  distributed  sensors.  This  was 
done  by  using  the  FEM  interpolation  function 
to  calculate  the  distributed  parameters  of  the 
structure  and  the  resulting  response  of  the 
spatially  continuous  distributed  sensors.  The 
comparison  of  the  two  simulations  confirmed 
the  shape  function  calculation  of  the 
distributed  sensors. 

The  calculation  of  the  sensor  shape 
functions  is  a  mathematically  and 
computationally  intensive  process.  In  this 
procedure  it  is  required  to  calculate 
displacement  and  displacement  rate  kernels. 
These  kernels  are  then  transformed  into  the 
desired  kernels  by  the  appropriate 
mathematical  operations.  For  the  PVDF 
sensor  implementation,  this  requires  a  double 
integration  of  the  displacement  kernel  to 
obtain  the  curvature  kernel  and  a  double 
integration  of  the  displacement  rate  kernel  to 
obtain  the  curvature  rate  kernel.  A  method  is 
presented  for  the  direct  calculation  of  the 
curvature  and  curvature  rate  kernels.  This 
reduces  the  computation  time  required  for  the 
distributed  sensor  shape  optimization. 

Completely  covering  the  members  of  a 
structure  with  two  distributed  sensors  may  be 
possible  in  a  laboratory  setting,  but  it  may  not 
be  practical  for  a  large  structure  with  many 
structural  elements  and  physical  components. 
A  solution  to  this  problem  was  found  by  using 
a  relatively  small  part  of  one  of  the  distributed 
sensors  along  with  an  observer.  The  output  of 
the  resulting  distributed  sensor  provides  a 
portion  of  the  desired  control  signal.  The 
observer  may  be  driven  by  the  output  of  the 
implemented  sensor  providing  the  states  of  the 
state  space  model.  The  states  may  then  be  used 
to  generate  the  missing  part  of  the  control 
signal  that  would  be  provided  by  the  full 


sensor.  This  method  increases  the  amount  of 
computation  that  must  be  done  by  the  control 
system,  but  it  does  allow  for  the 
implementation  of  the  desired  control  with  a 
distributed  sensor  covering  a  part  of  the 
structure.  The  usefulness  of  a  single  relatively 
small  customized  PVDF  film  sensor  along 
with  an  observer  to  provide  full  state  feedback 
information  for  the  entire  structure  is 
examined.  The  complexity  and  performance  of 
the  proposed  control  scheme  is  compared  with 
the  original  distributed  sensor  implementation. 

Since  smart  structures  exist  in  a  variety 
of  sizes  and  shapes  and  in  various  applications, 
it  is  desirable  to  have  many  different  control 
design  algorithms  available  to  the  control 
design  engineer.  For  many  of  the  potential 
applications,  controller  design  methods  other 
than  the  LQR  may  provide  superior 
performance.  As  a  first  step,  a  pole  placement 
controller  was  designed  to  determine  if  shape 
functions  for  this  controller  could  be 
calculated  for  the  distributed  sensors.  The 
shape  functions  obtained  were  satisfactory 
indicating  that  state  space  controllers  other 
than  the  LQR  may  be  implemented  with  the 
distributed  sensors. 

In  summary,  we  have  developed  an 
alternate  method  for  computing  curvature  and 
curvature  rate  kernels  in  the  development  of 
optimized  sensor  shape  functions  and 
proposed  a  procedure  for  utilizing  a  partial 
distributed  sensor  along  with  an  observer  for 
the  implementation  of  controllers.  We  have 
also  examined  various  control  methodologies 
for  designing  distributed  controllers.  The 
usefulness  of  a  single  relatively  small 
customized  sensor  along  with  an  observer  to 
provide  full  state  feedback  information  is 
demonstrated  with  simulation  studies.  We  are 
in  the  process  of  implementing  these  sensors 
and  control  algorithms  on  a  simple  cantilever 
beam  tes*  ^icle. 
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ABSTRACT 

A  mathematical  model  is  presented 
for  the  dynamic  response  of  contractile 
fiber  bundles  embedded  in  or  around  elastic 
springs  that  are  either  linear  helical 
compression  springs  or  hyperelastic  springs 
such  as  rubber-like  materials.  The  fiber 
bundle  is  assumed  to  consist  of  a  parallel 
array  of  contractile  fibers  made  form  either 
contractile  polymeric  muscles  such  as 
polyacrylic  acid  plus  sodium  acrylate 
cross-linked  with  bisacrylamide  (PAAM) 
or  polyacrylonitrile  (PAN)  fibers.  The 
proposed  model  considers  the  electrically  or 
pH-induced  contraction  of  the  fibers  which 
may  include  resistive  heating  of  the  fiber  in 
case  of  shape-memory  alio  "he  theory 
then  branches  out  to  wr  different 
description  of  such  process  for  ionic 
polymeric  gel  fiber  bundles  and 
shape— memory  alloy  fiber  bundles. 

INTRODUCTION 

Intelligent  material  systems  and 
structures  have  become  important  to  all 
defense  and  civilian  sectors  of  the  society 
(see  Ahmad,  Crowson,  Rogers  and  Aizawa 
[1]).  Accordingly,  based  on  such  materials, 
structures  and  their  integration  with 
appropriate  sensors  and  actuators,  novel 
applications,  useful  for  all  defense  and 
civilian  programs,  have  emerged.  Numerous 
examples  may  be  cited  as  can  be  seen  from 
key  reference  articles  such  as  [1]  and 


Crowson  and  Anderson  [2].  In  the  present 
paper  a  mathematical  model  is  presented 
for  the  dynamic  response  of  contractile 
fiber  bundles  embedded  in  or  around  elastic 
springs  that  are  either  linear  helical 
compression  springs  or  hyperelastic  springs 
6uch  as  rubber— like  materials.  The 
proposed  theory  presents  two  different 
description  of  such  processes  for  ionic 
polymeric  gel  fiber  bundles  and 
shape-memory  alloy  fiber  bundles. 

For  the  case  of  ionic  polymeric  gel  fibers 
the  model  consists  of  an  encapsulated 
hermetically  sealed,  helical  compression 
spring-loaded  cylindrical  linear  actuators 
containing  a  counterionic  solution  or 
electrolyte  such  as  water+acetone,  a 
cylindrical  helical  compression  spring  and  a 
collection  of  polymeric  gel  fibers 
(polyelectrolytes)  such  as  polyvinyl  alcohol 
(PVA)  polyacrylic  acid  (PAA)  or 
polyacrylamide.  Furthermore,  the  helical 
spring  not  only  acts  as  a  compression 
spring  between  the  two  hermetically  sealed 
circular  end-caps  but  contains  snugly  the 
polymeric  gel  fiber  bundle  and  also  acts  as 
the  cathode  (anode)  electrode  while  the 
two  actuator  end-caps  act  as  the  other 
cathode  (anode)  electrodes.  In  this  fashion, 
a  DC  electric  field  of  a  few  volts  per 
centimeter  can  cause  the  polymer  gel  fiber 
bundle  to  contract  (expand).  This  causes 
the  compression  spring  to  contract  and  pull 
the  two  end-caps  closer  to  each  other 


against  the  elastic  resistance  of  the  helical 
spring.  By  reversing  the  action  by  means  of 
reversing  the  electric  field  polarities  the  gel 
is  allowed  to  expand  while  the  compression 
spring  is  also  expanding  and  helping  the 
linear  expansion  or  the  actuator  since  the 
polymeric  gel  muscle  expands  due  to  the 
induced  alkalinity  along  the  helical  spring 
body.  Thus,  electrical  control  of  the 
expansion  and  the  contraction  of  the  linear 
actuator  will  be  possible.  A  mathematical 
model  is  presented  based  on  the  proposed 
composite  structure  that  takes  into  account 
all  pertinent  variables  such  as  the  pH  of 
the  gel  fiber  bundle,  the  pH  of  the 
surrounding  medium,  the  hyperelastic 
parameters  of  the  fiber  bundle,  the 
electrical  variables  of  the  gel,  the  electric 
field  strength,' the  pH  field  strength  and  all 
pertinent  dimensions  followed  by  some 
numerical  and  experimental  simulations 
and  data. 

For  the  second  model,  we  consider  the  fiber 
bundle  of  SMA  to  be  either  circumscribed 
inside  a  helical  compression  spring  with  flat 
heads  or  in  parallel  with  a  number  helical 
compression  springs,  end-capped  by  two 
parallel  circular  plates  with  embedded 
electrodes  to  which  the  ends  of  the  SMA 
fibers  are  secured.  Thus,  the  fibers  can  be 
electrically  heated  and  subsequently 
contracted  to  compress  the  helical 
compression  spring  back  and  forth.  Design 
details  are  first  described.  In  essence  the 
dynamic  behavior  of  the  actuator  depends 
on  the  interaction  between  the  current 
supplied  to  the  wires  and  the  heat  transfer 
from  the  wires.  Further,  a  mathematical 
model  is  presented  to  simulate  the 
electro-thermo-mechanics  of  motion  of 
such  actuators.  The  proposed  model  takes 
into  account  all  pertinent  variables  such  as 
the  strain  c,  the  temperature  of  the  fibers 
T(t)  as  a  function  of  time  t,  the  ambient 
temperature  Tq,  the  martensite  fraction  £, 

the  helical  compression  spring  constant  k 
and  the  overall  heat  transfer  coefficient  h. 
Numerical  simulations  are  then  carried  out 
and  the  results  are  compared  with 
experimental  observations  of  a  number  of 
fabricated  systems. 
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Figure  1—  Composite  Fiber-Bundle 
Structure 


BUCKLING  AND  COLLAPSE  OF  SUDDENLY  LOADED  STRUCTURES 
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Problems  which  deal  with  stability  (or 
instability)  of  motion  have  concerned  researchers 
for  many  years  in  many  fields  of  engineering.  In 
structural  mechanics,  dynamic  stability  or 
instability  has  received  considerable  attention  in 
die  past  40  years.  Several  studies  have  been 
conducted  by  various  investigators  on  structural 
systems  which  are  dynamically  loaded  (sudden 
loads,  periodic  loads  or  nonperiodic  time- 
dependent  loads).  In  these  studies,  several 
attempts  have  been  made  to  define  critical 
conditions  and  to  develop  methodologies  for 
estimating  critical  conditions.  In  reviewing  these 
efforts  one  easily  concludes  that,  in  structural 
mechanics,  "dynamic  stability"  encompasses 
many  classes  of  problems,  and  many  physical 
phenomena.  It  is  not  surprising  then,  that  the 
term  finds  several  uses  and  interpretations, 
among  Structural  n^rhanioian^ 

Problems  of  parametric  resonance, 
follower  force  problems,  and  fluid-structure 
interaction  problems  are  all  listed  under  tire 
heading  of  "dynamic  stability",  among  others. 
See  the  review  article  of  Herrmann  (1967). 

Finally,  a  large  class  of  structural 
problems  that  has  received  considerable  attention 


and  does  qualify  as  a  category  of  dynamic 
stability  is  that  of  impulsively  loaded 
configurations  and  configurations  that  are 
suddenly  loaded  with  loads  of  constant  magnitude 
and  infinite  duration.  These  configurations  under 
static  loading  are  subject  to  either  limit  point 
instability  or  bifurcational  instability  with  an 
unstable  postbuckling  branch  (violent  buckling). 
The  two  types  of  loads  may  be  thought  of  as 
mathematical  idealizations  of  blast  loads  of  (1) 
large  decay  rates  and  small  decay  times  and  (2) 
small  decay  rates  and  large  decay  times, 
respectively.  For  these  loads,  the  concept  of 
dynamic  stability  is  related  to  the  observation 
that  for  sufficiently  small  values  of  die  loading, 
the  system  simply  oscillates  about  the  near  static 
equilibrium  point  and  the  corresponding 
amplitudes  of  oscillation  are  sufficiently  small. 
If  die  loading  is  increased,  same  systems  will 
experience  large-amplitude  oscillations  or,  in 
general,  a  divergent  type  of  motion.  For  this 
phenomenal  to  happen,  the  configuration  must 
possess  two  or  more  static  equilibrium  positions, 
and  escaping  motion  occurs  by  having 
trajectories  that  can  pass  near  an  unstable  static 
equilibrium  point  Consequently,  the 


methodologies  developed  by  the  various 
investigators  are  for  structural  configurations 
that  exhibit  snap-through  buckling  when  loaded 
quasi-statically.  See  die  book  by  Simitses  (1989) 
and  the  cited  references. 

In  addition,  there  is  considerable  work 
reported  in  the  literature  that  deals  with  suddenly 
loaded  columns  and  plates.  The  beginning  of  the 
study  of  these  problems  can  be  traced  to  the 
investigation  of  Koning  and  Taub  (1933)  who 
considered  the  response  of  an  imperfect  ("half- 
sine"  initial  shape),  simply  supported  column 
subjected  to  a  sudden  axial  load  of  specified 
duration.  Since  these  geometries  under  static 
loading  exhibit  bifurcational  buckling  wife  stable 
postbuclding  strength  (smooth  buckling)  their 
behavior  is  different  from  feat  of  geometries 
which  exhibit  violent  static  budding.  Therefore, 
fee  concepts,  criteria  and  estimates  of  critical 
conditions  are  different.  In  all  cases  of 
geometries  which  are  subjected  to  sudden  loads 
fee  ultimate  interest  is  whether  or  not  dynamic 
buckling  will  lead  to  collapse.  This  is  a  very 
important  design  consideration.  For  more 
references  see  Simitses  (1987). 

The  presentation  will  deal  wife  concepts, 
criteria  and  estimates  of  critical  conditions  for 
this  class  of  problems.  Moreover,  examples  of 
systems  which  are  subject  to  this  type  of  loading 
and  therefore  subject  to  dynamic  buckling  will 
also  be  presented.  These  include  shallow  arches 
and  spherical  caps,  frames,  cylindrical  shells,  and 


many  others.  These  configurations  are 
representative  of  domes,  submarine  hulls,  aircraft 
fuselages  and  jet  engine  casings.  In  addition,  fee 
presentation  will  address  such  important  effects 
as  those  of  duration  time,  static  pre-loading  and 
damping.  The  emphasis  will  be  placed  on 
identifying  problems  in  search  of  a  solution  and 
future  research  directions. 

In  this  task,  special  consideration  will  be 
given  to  structural  configurations  made  out  of 
fiber  reinforced  composites,  such  as  laminated 
plates  and  shells. 
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US  Army  Research  Laboratory 
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For  fiber  reinforced  composite  materials, 
the  matrix  dominated  properties  (i.e.  shear  and 
transverse  tensile)  increase  with  strain  rate. 
Fibers  such  as  graphite  and  glass  fibers  are  in 
general  not  rate-sensitive  materials;  however, 
the  fiber  dominated  properties  may  increase  due 
to  the  effects  of  strengthening  matrix.  This  is 
especially  true  for  compressive  strength  and 
modulus  since  a  rigid  and  strong  matrix  can 
provide  lateral  supports  increasing  the  buckling 
strength  of  fibers.  Dynamic  behavior  of 
materials  is  typically  obtained  from  standard 
testing  machines  up  to  strain  rates  of  1-10  /sec 
and  from  Hopkinson  Bar  experiments  at  strain 
rates  greater  than  100  /sec  [1,2,3].  This 
generally  leaves  a  gap  in  strain  rates  effects 
in  the  10-100  /sec  band,  precisely  the  zone  of 
interest  for  launching  of  projectiles.  In 
addition  most  current  testing  is  performed  at  a 
constant  strain  rata  up  to  failure  of  the 
specimen.  While  this  procedure  does  give  a  good 
picture  of  the  increase  of  yield  and  ultimate 
stress  with  increasing  strain  rate  for  rate 
sensitive  materials,  it  is  unclear  as  to  whether 
this  increase  in  strength  at  failure  can  be 
effectively  used  to  prevent  failure  of  materials 
subjected  to  short,  high  rate,  impulsive  loads 
under  interior  ballistic  cycles. 

This  investigation  proposed  a  experimental 
set-up  and  procedures  including  an  Air  Gun 
system,  design  of  test  fixtures  and  data 
acquisition.  The  experiments  will  be  able  to 
simulate  the  loading  condition  of  tank  gun  and 
artillery  projectiles  during  launching.  The  Air 
Gun  system  consists  of  a  series  of  various 
diameter  tubes  (2, 3, 4, and  7  inches)  of  great 
lengths.  Figure  1  shows  an  overall  layout  of 
the  Air  Gun  system.  The  breech  area  of  the  gun 
is  located  on  the  right  hand  side  where  a 
projectile  is  loaded.  The  gun  barrel  is  then 
evacuated  to  1  Torr  and  the  projectile  is 
released.  The  bird  is  accelerated  by  atmospheric 
pressure  (or  slightly  augmented)  and  travels  for 
several  hundred  feet  (up  to  300  feet  for  the  7 
inch  gun)  to  the  catch-box  area  in  the  terminal 
velocity. 

After  muzzle  exit,  the  bird  impacts  a 
mitigator  which  is  fabricated  from  honeycomb 
material  and  can  be  geometrically  tailored  for 
specific  loading  rates.  The  momentum  exchange 
mass  (MEM)  is  a  piece  of  heavy  steel  to  absorb 
the  residual  energy  from  impact.  Figure  2 
illustrates  the  mechanism  of  impact  tests.  By 
careful  engineering  of  theses  features,  it  is 
generally  possible  to  create  realistic  interior 


ballistic  loading  rates,  pulse  shapes  and 
durations  during  the  deceleration  of  the  bird. 
The  specimen  is  held  by  a  fixture  between 
mitigator  and  MEM  and  an  extra  mass  may  be  added 
in  front  of  the  fixture  to  adjust  the  force  in 
the  specimen.  The  speed  of  projectile  and  force 
acting  on  each  component  of  the  system  can  be 
predicted  from  modeling  of  dynamic  response 
during  impact. 

Fixtures  are  designed  and  fabricated  to 
hold  the  specimen  in  correct  alignment  during 
impact  and  provide  uniform  stress  transfer  into 
the  specimen.  The  gage  length  is  adjustable  and 
test  results  can  be  compared  to  those  from 
standard  static  test  fixtures  such  as  IXTRI  or 
Celanese  configurations  by  ASTM  Standard  D  3410 
[4,5].  A  0.5  inch  gage  length  is  used  for  this 
specific  investigation.  The  composite  test 
coupons  is  about  1.00  inch  wide  and  0.09  inch 
thick  with  a  layup  construction  of  [ (0/45/- 
45/0)«] .  Im7/8551-7  composite  prepreg  was  used 
for  fabrication  .  Figure  3  shows  a  typical 
failure  of  specimen  after  an  impact  test. 

Measurements  of  displacements  and 
deformations  were  performed  by  a  streak  cameras . 
A  strip  pattern  was  printed  on  the  surface  of 
specimen  and  a  high  speed  camera  was  aimed  at 
the  strip  constantly  during  impact.  Figure  4 
shows  a  streak  film  taken  by  the  camera  with 
speed  of  120  ft/sec.  The  total  impact  duration 
of  this  particular  test  is  about  5  mil  second. 
The  parallel  white/black  strip  in  the  film  is 
the  image  of  the  strip  pattern.  It  is  noted  that 
the  movement  of  strip  is  a  rigid  body  motion  of 
the  specimen  while  the  contraction  of  strip  is 
compressive  deformation  due  to  impact.  The 
streak  film  was  then  digitized  and  analyzed  by 
using  a  high  resolution  PDS  1010A  micro- 
densitometer  and  PDS  2300C  data  acquisition 
system.  Accordingly,  the  deformation  state  in 
terms  of  time  was  obtained  and  the  strain  due  to 
impact  could  be  calculated.  Figure  5  shows  the 
strain  as  a  function  of  time.  A  very  sharp  drop 
occurs  at  end  of  impact  (about  4  msec)  where 
material  failed.  The  ultimate  compressive  strain 
is  about  1%. 

An  accelerometer  was  placed  on  the  MEM 
surface  to  record  the  acceleration  in  terms  of 
impact  time.  The  force  in  the  specimen  is 
assumed  to  be  equal  to  the  force  acting  on  the 
MEM.  Therefore,  stress  in  the  specimen  can  be 
derived  and  is  equal  to  the  product  of  the 
acceleration  and  the  mass  of  MEM.  The  force  as  a 
function  of  impact  time  is  illustrated  in  Figure 


6.  The  zigzag  pattern  in  the  curve  is  due  to  a 
stress  wave  travelling  in  the  MEM.  The  maximum 
stress  in  the  specimen  can  be  calculated  to  be 
220,  165  and  122  Jcsi  at  points  A,  B  and  c, 
respectively.  The  static  ultimate  strength  was 
measured  to  be  115  Jcsi  by  using  ITTRI  test 
method  with  same  gage  length. 

Strain  rate  effects  on  compressive 
properties  of  composite  are  clearly  observed. 
In  fact,  the  tested  specimen  shows  the  failure 
is  initiated  by  transverse  delamination  then 
followed  by  fiber  buckling  (Figure  3) .  since  the 
polymer  matrix  is  highly  rate  dependent 
material,  the  dynamic  strength  increases  as 
strain  and  loading  rates  increase. 
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Figure  3.  Specimen  after  impact 


Figure  4.  Streak  film 
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Figure  5.  compressive  strain  in  the  specimen 
during  impact 


Figure  1.  Air  Gun  System 
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Figure  6.  Force  resulted  from  impact 


Technology  Programs  in  Composite  Structures 
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The  U.S.  Army  Missile  Command 
has  a  number  of  projects  in  the 
areas  of  advanced  composite 
structures.  Each  project  has 
the  ultimate  goal  of  improving 
weaponry  for  current  and  future 
military  systems.  Most  of  the 
projects  described  below  span 
the  FY93-95  time  frame. 

One  ongoing  project  at 
MICOM  is  3-inch  bottle  testing 
used  to  determine  ultimate 
strengths  used  in  design  of 
composite  structures .  The 
method  used  for  fabricating 
mandrels  is  tedious  and  time 
consuming.  This  program  will 
develop  a  more  efficient  method 
of  fabricating  bottles  in 
quantities  consistent  with 
program  needs. 

Proof  testing  of  composite 
structures  is  currently  the 
method  used  at  MICOM  to 
determine  the  quality  of  rocket 
motorcases  and  launch  tubes. 
Taguchi  methods  will  be 
utilized  to  determine  the 
effects  of  the  rate  of 
pressurization,  the  rate  of 
proof  load  to  burst  pressure, 
and  the  fiber/ res in 
combinations  on  composite  test 
vessels . 

Incorporation  of  optical 
fibers  into  a  composite 
structure  for  strength, 
pressure,  and  temperature 


measurement  is  currently 
underway  at  MICOM.  Optical 
fibers  will  be  incorporated 
into  the  3 -inch  pressure 
vessels  to  determine 
practicality.  The  specimens 
will  then  be  hydroburst  tested 
to  determine  if  structural 
integrity  is  compromised  due  to 
the  optical  fiber.  Analysis  on 
the  microstructure  will  be 
performed  to  determine  effects 
on  the  matrix. 

An  effort  is  underway  in 
the  area  of  service  life 
predictions.  A  literature 
search  will  be  performed  to 
obtain  data  on  long  term 
exposure  of  composites.  A 
series  of  similar  test 
specimens  will  be  fabricated 
and  exposed  to  similar 
conditions  at  accelerated 
times.  The  specimens  will  be 
tested  in  the  same  manner  as 
those  seeing  long  term  exposure 
and  a  correlation  between  the 
two  tests  will  be  developed. 

The  determination  of 
critical  fiber  length  in  resin 
transfer  molding  is  a  program 
underway  at  MICOM.  Specimens 
will  be  made  and  tested  in 
tensile  and  flexure  modes.  The 
test  specimens  will  be  of 
simple  geometry  and  consist  of 
various  chopped  fiber  lengths 
mixed  with  resin  and  injected 


into  a  mold.  Data  from  the 
tensile  and  flex  tests  will  be 
plotted  (load  vs.  fiber  length) 
thus  determining  a  "critical 
fiber  length"  for  an 
anticipated  load  condition. 

Another  program  underway 
at  MICOM  is  thermoplastic 
coatings  for  use  on  composite 
and  metal  substrates  as  a 
substitute  for  CARC  paint.  A 
series  of  both  metal  and 
composite  substrates  will  be 
flame  sprayed  using  selected 
thermoplastics.  The  coating  and 
material  combinations 
undergoing  successful 
application  will  be  tested  for 
coating  integrity  to 
environmental  exposure, 
adhesion  and  flexibility.  ASTM 
tests  will  be  used. 

A  high  use  temperature 
composite  structures  program  is 
underway  at  MICOM.  Polyimide 
resins  and  phenolic  triazine 
are  the  primary  candidates  that 
have  been  manufactured  to  date. 
The  MICOM  Advanced  Kinetic 
Energy  Missile  System  composite 
booster  motorcase  and 
centerbody  sections  have  been 
manufactured  by  four 
contractors  and  are  currently 
being  tested  for  bending  and 
ultimate  strength  at  MICOM. 

A  technology  effort  is 
underway  at  MICOM  on  low  impact 
damage  of  graphite/epoxy 
composite  structures .  Taguchi 
methods  are  being  utilized  to 
determine  design  enhancements 
which  can  be  used  to  mitigate 
impact  damage  and  increase  the 
survivability  of  a  variety  of 
composite  structures  such  as 
motorcases  and  launch  tubes. 
Additional  layers,  either 
filament  wound  or  braided  will 
be  added  to  a  baseline 
structure.  The  structures  will 


then  be  impacted  at  different 
levels  of  energy  and  hydroburst 
to  determine  ultimate  strength. 

A  composites  adhesive 
program  is  underway  that 
characterizes  adhesives  via 
rapid  heating  rates.  Double 
lap  shear  joint  test  specimens 
will  be  fabricated  with 
incorporated  thermocouples  near 
the  adhesive  joints.  Tests 
will  be  performed  at  both  high 
heating  and  loading  rates  to 
simulate  environements  that 
adhesives  encounter  when  used 
in  missile  systems. 

Each  one  of  the  above 
described  programs  is  being 
performed  jointly  by 
Government,  universities  and 
industry  to  encompass  all  the 
areas  of  expertise  to  assist  in 
solving  problems  encountered  in 
the  completion  of  the  tasks. 
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ABSTRACT 


Damage  mechanics  is  used  to 
investigate  the  various  damage 
mechanisms  in  loaded  metal  matrix 
composite  plates.  These  structural 
components  are  widely  used  in  the 
Army's  vehicles  and  installations.  They 
are  constantly  subjected  to  different 
types  of  loads  which  give  rise  to  the 
initiation  and  growth  of  various  damage 
mechanisms.  These  mechanisms  range 
in  scope  from  matrix-related  to  fiber- 
related.  In  addition,  certain  damage 
mechanisms  arise  due  to  the  effect  of  the 
interaction  between  the  matrix  and  fibers. 
In  general,  these  include  void  growth  and 
coalescence,  microcrack  initiation,  crack 
propagation,  fiber  fracture,  debonding 
and  delamination. 

The  damage  model  is  based  on  the 
concept  of  effective  stress  introduced  by 
Kachanov  (1958)  and  the  recent  work  of 
Murakami  (1988)  .  A  new  damage 
tensorial  variable  is  introduced  to  model 
crack  initiation  and  propagation.  The 
definition  of  this  tensor  is  based  on 
experimental  measurements  of  the  crack 
density  at  mutually  orthogonal  cross- 
sections  of  the  damaged  component.  The 
damage  model  is  coupled  with  existing 
material  model  to  formulate  an  integrated 
damage  characterization  constitutive 
model  for  metal  matrix  composite  plates. 
The  resulting  coupled  material  model  is 
capable  of  analyzing  both  elastic  and 
elasto-plastic  metal  matrix  composites. 


However,  in  all  cases,  the  brittle  fibers 
are  assumed  wO  be  elastic. 

The  measured  crack  densities  are 
normalized  with  respect  to  the  crack 
densities  at  the  failure  load.  The  results 
are  shown  in  terms  of  stress-strain 
curves  for  the  damaged  components.  In 
addition,  damage  curves  are  shown 
representing  plots  of  the  damage  variable 
components  versus  the  strain.  Certain 
conclusions  and  recommendations  are 
made  based  on  the  results.  In  particular, 
the  critical  state  of  damage  is  located 
and  identified. 

Two  approaches  are  considered  in  the 
formulation  of  the  damage  model  with 
respect  to  the  damage  of  the  matrix  and 
fibers  as  shown  in  Figure  1.  The  first 
approach  is  termed  "overall"  in  the  sense 
that  one  damage  variable  is  used  to 
characterize  damage  in  both  the  matrix 
and  fibers  (Kattan  and  Voyiadjis,  1993). 
The  second  approach  is  termed  "local"  in 
the  sense  that  two  damage  variables  are 
used  to  model  the  damage  mechanisms 
in  the  matrix  and  fibers;  one  for  each 
constituent  (Voyiadjis  and  Kattan,  1993). 
Each  approach  has  its  own  advantages 
and  disadvantages.  In  the  first  approach, 
utilizing  one  damage  variable  has  the 
benefit  of  capturing  all  damage  effects  in 
the  system.  However,  this  approach  lack 
a  detailed  description  of  each  damage 
mechanism  as  it  lumps  them  all  into  one 


Overall  and  Local  Approaches. 


variable.  On  the  other  hand,  the 
local  approach  has  the 
advantage  of  considering  "local" 
or  "constituent"  damage 
mechanisms  as  it  uses  an 
independent  variable  for  each. 
However,  its  weak  point  lies  in 
the 

characterization  of  interfacial 
damage;  that  is  it  cannot 
describe  properly  the  damage 
arising  due  to  the  interaction  of 
the  matrix  and  fibers,  like 
decohesion  and  debonding. 


In  the  derivation  of  the 
constitutive  model,  it  is  shown 
that  for  the  elasto-plastic  model, 
an  anisotropic  yield  function  is  derived  for 
the  damaged  composite  system.  This 
results  from  the  consideration  a  von 
Mises  type  yield  function  for  the 
undamaged  matrix  material  with  elastic 
fibers.  In  a  similar  fashion,  it  is  also 
shown  how  the  flow  rule  for  the  plastic 
strain  is  non-associated  for  the  damaged 
composite  system,  where  an  associated 
flow  rule  is  assumed  for  the  undamaged 
matrix  material. 

In  order  to  validate  the  theoretical 
predictions,  a  number  of  uniaxial  tensile 
loading  experiments  are  performed  on 
titannium  aluminide  SiC-reinforced  metal 
matrix  composite  specimens.  Two 
different  laminate  layups  are  considered; 
(0/90),  and  (±45),  .  Only  one  of  each 
orientation  was  loaded  to  failure.  The 
remaining  specimens  were  loaded  at  up 
to  certain  percentages  of  the  failure  load. 
The  damaged  specimens  are  then 
analyzed  using  image  analyses 
techniques.  Crack  densities  are 
measured  for  all  the  damaged 
sepcimens.  The  damage  tensor  is 
computed  based  on  these  measurements 
and  the  results  are  compared  with  the 
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ABSTRACT 


Beams  are  well  recognized  as  a  basic  structural 
dement  Engineering  students  are  introduced  to  classical 
beam  theory  in  conjunction  with  sophomore  strangth-of- 
materials  which  traditionally  deals  with  homogeneous 
metals.  Increased  use  of  composite  materials  in 
engineering  applications  has  induced  the  need  to  develop 
analytical  models  applicable  to  laminated,  anisotropic 

Initial  development  of  laminated  beam  theory  can  be 
found  in  elementary  stiength-of-materials  where  beams 
constructed  of  different  materials  are  considered.  In  this 
context  laminated  beams  are  analyzed  within  the 
framework  of  classical  beam  theory.  In  particular,  each 
layer  is  assumed  to  be  raider  a  state  of  onedimensional 
stress  (Le.  axial  normal  stress  and  transverse  shear  stress 
are  the  only  nan-vanishing  stresses).  For  materials  with 
similar  Poisson  ratios,  this  assumption  is  valid. 
However,  for  highly  anisotropic  materials  with  a  severe 
mismatch  of  Poisson  ratios  from  ply-to-ply  a  serious 
transverse  strain  incompatibility  results.  This  issue  has 
been  initially  addressed  by  considering  bidirectional 
laminates  under  uniaxial  tension  loading  (Dietz,  1949). 
Classical  strength-of-materials  framework  has  aim  been 
utilized  for  analyzing  laminated,  anisotropic  beams  by 
Vinson  and  Sierakowski  (1986)  and  move  recently  by 
Vasiliev  (1993). 

Other  approaches  to  laminated  beam  theories  have 
enforced  transverse  strain  compatibility.  However, 
confusion  exists  in  these  onedimensional  approaches 
between  cylindrical  bending  and  a  beam  theory.  The 
difference  between  die  two  models  is  analogous  to  the 
difference  between  plane  stress  and  pie  e  strain  in 
classical  theory  of  elasticity  and  is  clearly  distinguished 
by  Whitney  (1987).  Cylindrical  boding  has  been 
confused  with  beam  theory.  A  valid  i«mm«hirf  beam 
theory  should  reduce  to  classical  beam  theory  for  the  case 
of  a  homogeneous  material.  This  does  not  occur  with 
cylindrical  bending. 

Width  effects  produce  additional  difficulties  with 
laminated  composite  beams  due  to  the  large  width-to- 
thkkness  ratio  of  many  applications,  including  current 
flexure  test  specimens.  In  particular,  laminated  beams 
are  closer  to  a  narrow  plate  than  to  a  classical  beam. 
Shear  coupling  produced  with  die  inclusion  of  angle-ply 
layers  in  the  laminate  make  one-dimensional 
aasumptions  leas  valid  (Kedward,  1972). 

fothecurrent  paper  the  effect  of  kngth-to-width  ratio 
and  width-to- thickness  ratio  on  the  response  of 
anisotropic  laminated  beams  is  investigated.  The 
parameter  under  consideration  is  the  fundamental  free 


vibration  frequency.  Exact  solutions  are  obtained  for  the 
case  of  specially  orthotropic  thin  plates  with  two 
opposite  edges  simply-nipported  and  the  two  adjacent 
edges  free.  Minimum  leogth-to-width  ratios  compatible 
with  ooe-dimenaional  beam  assumptions  are  established 
for  having  a  high  major  Poisson's  ratio  raider 

bending.  In  order  to  consider  angle-ply  laminates  within 
die  framework  of  specially  orthotropic  plates,  stacking 
geometries  with  very  weak  bending-twisting  coupling  are 
investigated. 

The  effect  of  cross-section  width  is  investigated  by 
considering  long  narrow  (dates  with  various  widtb-to- 
thickness  ratios.  In  this  case  it  is  necessary  to  modify 
die  thin  plate  theory  to  include  transverse  shear 
deformation  in  the  width  direction. 

Beam  Theory  versus  Cylindrical  Bending 

Consider  a  laminated  (date  having  infinite  length  in 
the  y-dfrection  with  uniform  boundary  conditions  along 
die  edges  (see  Figrae  IX  If  any  surface  or  inplane  loads 
are  independent  of  y,  then  the  displacements  and  die 
corresponding  force  and  moment  resultants  are 
independent  of  y.  This  theory  is  not  compatible  with 
classical  beam  theory  as  force  and  moment  resnltants  are 
induced  relative  to  die  y-directian. 


Figure  1.  Cylindrical  bending. 


A  classical  beam  is  characterized  by  a  plate  having  a 
high  width -to-tinckness  ratio  with  the  long  sides  free  as 
illustrated  in  Figure  2.  In  this  case  the  force  and  moment 
resnltants  relative  to  the  y-dfrection  vanish.  The  inplane 
stresses  within  die  plies,  however,  are  two-dimensional 
in  nature. 
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Free- Vibration  Analysis 

Using  classical  laminated  plate  theory  based  on  the 
Kircbhoff  assumptions,  an  exact  solution  can  be  obtained 
for  a  specially  ortbotropic  laminated  plate  with  opposite 
sides  simply-supported  and  the  adjacent  edges  free. 
Numerical  results  are  obtained  for  [±45°]gs  and 
[±45°/0°2]js  laminates  with  ply  properties  which  are 
typical  of  contemporary  graphite/epoxy  laminates. 


U-L 


Figure  2.  Classical  beam. 


Both  of  these  laminates  contain  bending-twisting 
coupling  which  induces  foe  laminate  to  behave  as  an 
anisotropic  plate  rather  than  a  specially  orthotropic  plate. 
However,  the  coupling  is  very  weak  for  these  laminate 
stacking  geometries  and  foe  fundamental  vibration 
frequency  can  be  accurately  obtained  by  treating  foe  plate 
as  ortbotropic.  The  solution  is  obtained  by  seperation  of 
variables.  An  iteration  scheme  in  conjunction  with  a 
transcendental  equation  is  required  in  order  to  obtain  foe 
fundamental  vibration  frequency. 

Numerical  results  are  shown  in  Figure  3  for  the 
[±45°]w  with  10*  denoting  a  normalized  fundamental 
frequency. 


In  order  to  ascertain  foe  effect  of  widfo-to-foickness 
ratio  on  beam  response,  foe  Kirchhoff  theory  is  modified 
to  include  transverse  shear  deformation  in  foe  width 
direction.  Fundamental  frequencies  are  obtained  as  a 
function  of  widfo-to-thickness  rado6  for  plates  with  large 
values  of  R.  Length-to- thickness  ratios  are  small  so  that 
shear  deformation  is  not  required  in  foe  length  direction. 
Numerical  results  show  that  the  fundamental  frequencies 
are  rather  insensitive  to  widfo-to-thickness  ratios.  Thus, 
beam  behavior  is  controlled  by  R. 

References 

A.  G.  H.  Dietz,  Engineering  Laminates,  John 
Wiley,  1949. 

Keith  T.  Kedward,  "On  the  Short  Beam  Test 
Method,"  Fibre  Science  and  Technology,  Vol.  5, 1972, 
pp.  85-95. 

Valery  V.  Vasiliev,  Mechanics  of  Composite 
Structures ,  English  Edition  Editor  Robert  M.  Jones, 
Taylor  and  Francis,  Washington,  D.  C.,  1993. 

J.  R.  Vinson  and  R.  L.  Sicrakowski,  The  Behavior 
of  Structures  Composed  of  Composite  Materials, 
Martians  Nijhoff  Publishers,  Boston,  1986. 

J.  M.  Whitney,  "One-Dimensional  Analysis  of 
I  rnmutwi  Plates",  Chapter  4,  Structural  Analysis  of 
Laminated  Anisotropic  Plates ,  Technomic  Publishing 
Co_  Lancaster,  PA,  1987. 


Figure  3.  Normalized  fundamental  frequency. 


We  note  that  cylindrical  bending  (CB)  results  are 
obtained  for  length-to- width  ratios  (R)  lees  than  0.1, 
while  beam  theory  is  approached  for  R  >  8.  Similar 
reeulte  are  obtained  for  dm  [±45o/0o2l3S  laminate. 
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Viscoelastic  behaviour  is  much  more 
significant  in  polymeric  materials  than  in 
other  types  of  materials.  In  polymer  based 
composites,  the  viscoelastic  response 
manifests  strongly  in  matrix  dominated 
properties.  Service  conditions,  such  as 
stresses,  elevated  temperature  and  moisture 
can  accelerate  the  viscoelastic  processes  and 
thus  further  degrade  the  stiffness  and  strength 
of  the  composite.  However,  the  viscoelastic 
properties  are  not  always  harmful.  For 
instance,  the  damping  properties  of  the 
materials  can  be  utilized  in  controlling  the 
vibration  of  the  structures.  Fiber  reinforced 
polymer  composites,  incorporating  the  high 
stiffness  with  the  high  damping  properties, 
are  excellent  candidates  for  structures 
subjected  to  dynamic  loading. 

The  analysis  and  modeling  of 
viscoelastic  response  of  composite  structures 
are  based  on  the  viscoelastic  Correspondence 
Principle.  The  viscoelastic  properties  of  the 
composite  need  to  be  determined 
experimentally  and  then  modeled  by  a 
viscoelastic  constitutive  equation.  The 
viscoelastic  properties  of  the  composite  may 
also  be  predicted  by  micromechanics 
approaches  from  individual  constituents .  The 
latter  allows  a  broad  range  selection  of 
material  combination  in  design. 

The  viscoelastic  constitutive  modeling 
as  well  as  the  extension  of  micromechanics 
models  to  the  time  domain  are  based  on 
experimental  characterization  of  real  resin  and 
composite  systems.  This  paper  presents 
experimental  results  on  the  viscoelastic 
properties  of  PEEK  resin  and  PEEK 
composite  APC-2.  The  viscoelastic  response 
of  the  composite  is  compared  with  that  of  the 
pure  resin.  The  constitutive  modeling  for  the 


linear  and  nonlinear  viscoelastic  responses  is 
discussed. 

The  viscoelastic  properties  of  PEEK 
resin  arid  APC-2  composite  were 
characterized  by  creep  and  creep- recovery 
tests  following  an  accelerated  testing 
procedure  based  on  the  time-temperature  and 
superposition  principle  (TTSP)  and  time- 
stress  superposition  principle  (TSSP)  [1]. 
Additional  information  was  obtained  by 
dynamic  mechanical  measurements. 

As  compared  with  the  pure  resin, 
APC-2  composite  shows  a  higher  damping 

over  the  P'  relaxation  region  and  a  lower 
temperature  for  the  glass  transition.  The 
discrepancy  is  likely  caused  by  the  interphase 
layer  in  the  composite.  Thus  the  two-phase 
micro-mechanical  model  may  not  be  adequate 
to  predict  the  viscoelastic  properties. 

The  recoverable  creep  response,  which 
is  required  for  viscoelastic  modeling,  can  be 
determined  through  proper  mechanical 
models  rather  than  mechanically  conditioning 
the  specimens.  PEEK  resin  displays  both 
instantaneous  plastic  deformation  and  time- 
dependent  flow.  A  mechanical  model  is 
proposed  to  represent  such  a  behaviour. 

The  Schapery's  nonlinear  viscoelastic 
equation  is  mostly  used  in  viscoelastic 
constitutive  modeling.  For  the  case  of 
uniaxial  loading  it  is  given  by  [2]: 


e(t)  =  g0S0o+  J  giAS(t|MiO 
V?  =  “MKO  = 


(1) 

where  go,  gi,  g2  and  a©  are  the  nonlinear 
parameters  with  clear  physical  meaning  [1], 

AS(tj>)  is  the  kernel  function  representing  the 
linear  viscoelastic  response,  and  t  is  time. 

The  kernel  function  may  be  described 
by  a  modified  power  law  containing  four 
parameters,  the  so  called  general  power  law. 
For  the  creep  compliance  S(t),  it  has  the 
following  form: 


dr 

3o 


V  =  t|>(t') 


-/$ 


S(t)  =  Son 


S«-So 

(l+SL)n 


where  to  is  the  characteristic  time  parameter. 
So  and  Soo  are  the  initial  compliance  and  the 
compliance  at  infinite  time,  respectively,  and 
n  is  a  material  parameter.  It  has  been  shown 
that  the  general  power  law  can  describe  the 
viscoelastic  responses  over  a  broad  time 
domain,  such  as  fitting  the  master  curves. 
The  effect  of  temperature  can  be  readily 
considered  by  the  time-temperature  shift 
function  aT  and  trr,  which  are  known  from 
the  TTSP  procedure. 

For  die  case  of  uniaxial  loading,  under 

the  step  stress  input  o[H(t)-H(t-ti)],  the  creep 
response  predicted  by  the  Schapery's 
equation  with  a  kernel  of  general  power  law 


Soo-So 

e(t)  =  g0S0o+gig2O — -  t<ti 

(1+B&Z.)" 


As  can  be  seen,  it  is  impossible  to  determine 
gl,  g2  and  a o  from  creep  tests  only.  Lou  and 
Schapery  have  proposed  a  procedure  to 
determine  the  above  parameters  from  creep 
and  creep-recovery  tests  [2].  However,  the 
procedure  given  in  [2]  implies  that  the  creep- 
recovery  response  follows  the  Boltzmann 
superposition  principle.  PECK  resin  and  the 
composite  show  severe  nonlinearity  in  creep- 
recovery  response  and  thus  the  Schapery's 
procedure  cannot  be  used. 

It  has  been  shown  that  for  a 
graphite/epoxy  composite  gi  and  g2  vary 
with  stress  level  in  such  a  way  that  their 
product  is  nearly  unity  [13-  Some 
experimental  evidences  suggest  that  ap,  the 
nonlinear  parameter  for  the  dissipation 
element  in  Schapery's  model,  is  stress 
independent  and  thus  can  be  related  to  the 
time-temperature  shift  factor  ay. 
Consequendy,  the  stress  and  temperature 
dependent  nonlinearity  of  the  dissipation 


element  and  modulus  element  can  be 
decoupled  such  that  the  former  i  s 
characterized  by  ay  while  the  latter  is 
represented  by  parameters  ac  and  gig2-  If 
the  temperature  dependence  of  the  non¬ 
linearity  between  the  generalized  external 
force  and  internal  force  is  much  significant 
than  the  stress  dependence  of  the  modulus 
element,  gig2  becomes  the  shift  factor  by. 
Thus  Equation  3  becomes: 


£(t)=M2£( 

bp 
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thereby  the  nonlinear  parameters  go  and  a<, 
can  be  determined  from  creep  data  only. 

The  stress-dependent  viscoelastic 
responses  of  APC-2  composite  were 
measured  by  90’  and  15*  off-axis  creep  tests 
for  the  transverse  tensile  mode  and  shear 
mode,  respectively.  The  linear  parameters  in 
Eq.2  were  determined  by  a  graphical  TTSP 
procedure.  While  the  nonlinear  parameters 
go  and  a 0  were  determined  by  a  numerical 
TSSP  procedure.  The  temperature-time  shift 
factor  trr  follows  the  Arrhenius  relation.  The 
stress-time  shift  factor  a 0  appears  to  follow 
an  exponential  law,  which  can  be  traced  back 
to  the  energy  barrier  theory  for  the 
nonlinearity  of  viscosity.  With  the  equations 
and  parameters  presented,  the  linear  and 
nonlinear  viscoelastic  behaviours  of  APC-2 
composite  laminate  from  ambient  temperature 
up  to  the  glass  transition  temperature  can  be 
described  analytically. 
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